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1.0 INTRODUCTION 


This is the final report for the "Integrated Source and Channel 
Encoded Digital Communication System Design Study" (Contract No. NAS 9- 
15240). The contract is to evaluate the Ku-band system design, analyze 
S-band and Ku-band tracking techniques, and study payload communication 
techniques. 

As defined in the Statement of Work, the first task required 
the evaluation and assessment of the baseline Ku-band system design 
under development by Hughes Aircraft Corporation (HAC) and study the 
performance/compatibility problems associated with the SSO/TDRSS RF 
interface. More specifically, the particular Ku-band carrier, PN 
despreading, and symbol synchronization strategies, which are selected 
by HAC for implementation in the Ku-band transponder aboard the 
Orbiter, are to be assessed and evaluated from a systems performance .. 
viewpoint, verifying that system specifications are met. Furthermore, 
it was required that any critical areas impacting the detailed concep- 
tual and breadboard designs as they develop be identified, discussed 
and studied and appropriate recommendations made for parameter optimi- 
zation, as necessary. 

During the one-year period covered by the contract, Axiomatix 
continuously rendered support to HAC in all of the above areas, with 
the major thrust of its activity directed toward understanding, studying, 
and proposing potential solutions to the data asymmetry problem prevalent 
on the Mode 1 high data rate channel of the Ku-band return link. In 
this link, two selectable modes which provide three channels are avail- 
able. Channel 1 is used for operational data, whereas Channels 2 and 3 
may be assigned a wide variety of digital and analog signals. Addi- 
tional support was provided by Axiomatix to assess the anti-jam capa- 
bility of the SSO/TDRSS RF interface. The results of these activities 
are summarized in Section 2.0, with the detailed analyses, system per- 
formance comparisons, and tradeoffs included as Appendices A through D. 

As a part of the overall effort to assess the Ku-band system 
design and to evaluate possible system simplifications/modifications 
that may be achievable, the second task required a study of the design 
and implementation of tracking techniques which are suitable for incor- 
poration into the Orbiter Ku-band communication system. Emphasis was 
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placed on maximizing tracking accuracy and communication system flexi- 
bility while minimizing cost, weight, and system complexity of Orbiter 
and ground system hardware. 

Section 3.0 summarizes the results obtained for the performance 
of several implementations of biphase Costas loops when subjected to 
the unbalanced QPSK input. The subcarrier tracking loop performance 
is also summarized in Section 3.0. The detailed analysis of both 
carrier and subcarrier tracking loops is presented in Appendices E 
through L. 

Section 4.0 describes the payload communication study to assess 
the design and performance of the forward link and return link bent-pipe 
relay modes for attached and detached payloads. As part of this study, 
a design for a forward link bent-pipe is proposed which employs a 
residual carrier but which is tracked by the existing Costas loop. 
Detailed analysis of Costas loop tracking for signals with residual 
carrier is presented in Appendix L. Critical Ku-band system parameters 
have been identified and optimized such that cost and complexity have 
been minimized. Finally, Section 4,0 identifies the basic transmission 
capabilities that can be accommodated by both the forward link and 
return link bent-pipe relay modes. 
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2.0 KU-BAND SYSTEM DESIGN EVALUATION 

As defined in the Statement of Work, this task required Axiomatix 
to evaluate and assess the baseline Ku-band system design under develop- 
ment by Hughes Aircraft Corporation (HAC) and study the performance/com- 
patibility problems associated with the SSO/TDRSS RF interface. More 
specifically, the particular Ku-band carrier, PN despreading, and symbol 
synchronization strategies, which are selected by HAC for implementation 
in the Ku-band transponder aboard the Orbiter, are to be assessed and 
evaluated from a systems performance viewpoint, verifying that system 
specifications are met. Furthermore, Axiomatix was required to identify, 
discuss, and study any critical areas which impact the detailed concep- 
tual and breadboard designs as they develop and to make appropriate 
recommendations for parameter optimization, as necessary. 

During the one-year period covered by this contract, Axiomatix 
continuously rendered support to HAC in all of the above areas, with the 
major thrust of its activity directed toward understanding, studying, and 
proposing potential solutions to the data asymmetry problem prevalent on 
the Mode 1 high data rate channel of the Ku-band return link. Additional 
support was provided by Axiomatix to assess the anti-jam capability of 
the SSO/TDRSS RF interface. The results of these activities are summarized 
in what follows, with the detailed analyses, system performance compari- 
sons, and tradeoffs included as Appendices A through D. 

2.1 Impact of Data Asymmetry on Bit Error Probability 

Performance of the Ku-Band Return Link 

In the early months of the Task 1 effort, an area of concern 
developed over the data asymmetry* produced by the HAC design of the 
50 Mbps (100 Msps) rate 1/2 convol utional ly coded data channel of the 
Ku-band return link through the TDRSS. This data asymmetry was poten- 
tially capable of causing several dB degradation to the demodulated data 
at the TDRSS ground station. Therefore, a suitable asymmetry model had 
to be developed, the amount of degradation as a function of the amount 


The term "data asymmetry" as used here denotes the condition in 
which pulses of opposite polarity in the data stream unintentionally have 
unequal widths. 


of asymmetry had to be computed,, and the possible solutions to the 
problem had to be investigated. 

The high data rate link from the Shuttle Orbiter through the 
TDRSS to the ground takes NRZ symbols at 50 Mbps and encodes them with 
a rate 1/2, constraint length 7, convolutional code. The bit error rate 
performance of the convolutional decoder depends, among other things, on 
the symmetry of the modulation. Any asymmetry in the NRZ symbols entering 
the symbol synchronizer causes a misalignment in the symbol synchroniza- ■ 
tion clock which degrades the integrate-and-dump output and any soft or 
hard decisions derived from it for input to the decoder. For a specified 
degree of asymmetry (in terms of a. fraction of a symbol interval), the 
bit error rate degradation is dependent on the transition probability of 
the data. Clearly, if the data transmitted was either all ones or all 
minus ones, misalignment of the bit synchronization clock would have no 
degrading effect on th4' integrate-and-dump output since, for each symbol, 
this circuit would integrate up to its maximum value before being dumped. 
On the other hand, when the data is an alternating sequence, the worst 
case degradation results, since the transition which occurs at the end 
of each symbol in combination with the symbol synchronization clock 
misalignment prevents the integrate-and-dump output from reaching its 
maximum value. 

2.1.1 Asymmetry Models and Definitions 

To quantitatively determine the degrading effect of NRZ symbol 
asymmetry on error rate performance, one must develop a suitable asym- 
metry model which accurately describes the physical source from which 
the asymmetry originates. During the early definition phase of the 
study, much confusion reigned over the appropriate model and definition 
of data asymmetry to be used in meeting the performance specification. 
After expending considerable effort on resolving the differences among 
the various models and definitions proposed, it was determined that two 
possible models were appropriate and, provided that asymmetry was properly 
defined, either model produced the identical performance degradation due 
to this asymmetry. In the first model (the one adopted by Axiomatix), 
the assumption is made that +1 NRZ symbols are elongated by AT/2 (rela- 
tive to their nominal value of T sec) when a negative-going data 
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transition occurs and -1 symbols are shortened by the same amount when 
a positive-going data transition occurs.* Otherwise (when no transitions 
occur), the symbols maintain their nominal T sec width. Thus, AT repre- 
sents the relative difference in length between the elongated +1 and 
shortened -1 symbols. An example demonstrating this model is illustrated 
in Figure la. The second asymmetry model (used in [1,2]) makes the 
assumption that positive NRZ pulses are shortened whenever adjacent 
prises are negative. Thus, a given positive pulse preceded and succeeded 
by a negative pulse would be reduced in duration at both ends. Stated 
another way, a positive-going transition occurs early and a negative- 
going transition occurs late relative to the nominal transition time 
instants. Letting 6 rep.^esent the fractional (relative to the nominal 
bit duration T) increase in positive pulse length due to a single adjacent 
negative pulse, then for a given random data sequence, the shortest pulse 
v/ould have length T(1 - 26)., while the longest would have length T(l+26). 
Figure lb illustrates the application of the second asymmetry model to 
the same bit stream as that used in Figure la. 

Regardless of which asymmetry model is used, data asymmetry is 
defined as the difference in length between the shortest and longest 
pulses in the sequence divided by their sum. For the first model, this 
definition gives 


Asymmetry = n 


T(1+~) - T(l-f) 
T(l+f) + T(l-|) 


( 1 ) 


whereas for the second model, we get 


T(1 + 26) - T(1 - 26) 
T(1 + 26) + T(1 - 26)' 



( 2 ) 


In the absence of noise, the timing instants for the in-phase 
integrate-and-dump (i.e., the epoch of the symbol synchronization clock) 
are determined in Appendix A. For the first asymmetry model, these 
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Due to symmetry in the data itself, it is immateriaT whether the 
elongated pulse is of positive or negative polarity, and vice versa for 
the shortened pulse. 
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occur at t = T(n + -|-); n = 0,±l,±2, For the second model, it can easily 

be shown that, on the average, the symbol synchronizer will lock up at 
the nominal transition points of the equivalent symmetric data waveform, 
i.e., 0,±T,±2T,±3T, . . . . 

2.1.2 Error Probability Performance (Perfect Symbol Synchronization) 


Based on the foregoing definition of asymmetry and the accompany- 
ing clock misalignment for the two models, it is straightforward to 
show that, for random NRZ data, the average symbol error probability 
associated with hard decisions made on the in-phase integrate-and-dump 
output of the symbol synchronizer is given by 




JEs/Nq + i erfc 


Je^/NoO - n) 


+ jg erfc 


(3) 


where denotes the symbol energy, Nq the channel noise spectral den- 
sity, and 


. A 2 
erfc X = — 


f” U.2 

e"^ dt . 


/rr •' X 


(4) 


Table 1 contains the symbol energy-to-noise ratio degradations 
(in dB) for asymmetry values nxlOO of 3 %, 7 %, ] 0 %, 15 %, and 20?o, and 
Es/No=0, 0.75, and 1.5 dB. The values of E^/Nq selected correspond 

to bit energy-to-noise ratios E,/N„ = 3, 3.75, and 4.5 which correspond 

° ^ -3 -4 -5 

respectively to decoder bit error probabilities P[j“'10 > 10 and 10 . 

The degradations are the additional E /Np, required due to asymmetry to 

s u 0 

produce the same value of symbol error probability when n=0, i.e., P^ , 

where 

= lerfc (Ve^] . (5) 


It should be noted that the symbol energy-to-noise ratio degradations 
given in Table 1 assume no channel bandwidth limitation, i.e., ideal 
rectangular pulse shapes have been assumed for the NRZ data. Any 
rounding of the pulses caused by channel bandwidth limitation produces 
an additional symbol energy-to-noise ratio degradation over and above 
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that due to data asymmetry. However, the bandlimiting degradation and 
the data asymmetry do not add algebraically but, rather, combine in a 
way determined by the particular symbol synchronization implementation. 
More will be said about the combined effects of bandlimiting and NRZ 
data asymmetry shortly. 

Table 1, Performance Degradation Due to Data Asymmetry with Random Data 


Asymmetry {%) 

Es/Nq (dB) 

'b 

Degradation (dB) 

3 

0 

10~3 

0.135 


0.75 

10"4 

0.135 


1.50 

10"5 

0.135 

7 

0 

10'3 

0.333 


0.75 

lO-'^ 

0.337 


1.50 

10“5 

0.340 

10 

0 

10-3 

0.495 


0.75 

10-‘^ 

0.505 


1.50 

10-5 

0.517 

15 

0 

10-3 

0.799 


0.75 

10-"^ 

0.824 


1.50 

10-3 

0.854 

20 

0 

10-3 

1.149 


0.75 

10-"^ 

. 1.201 


1.50 

10-5 

1 .264 


2.1.3 Error Probability Performance (Imperfect Symbol Synchronization) 

In the previous section, we presented the effect of NRZ data asym- 
metry on the error probability performance of a convolutionally coded 
channel, namely, the 50 Mbps channel of the Ku-band return link. 

Inherent in those results v;as the assumption that the symbol synchro- 
nizer produced a perfect clock which locked up with a misalignment equal 
to half the asymmetry for asymmetry Model 1 and no misal ignment for 
Model 2. The effect of additional clock misal ignment (commonly referred 
to as symbol sync jitter) which is typical of practical symbol syn- 
chronizers, was the next item studied under this task. The results 
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of this study, which are presented in detail in Appendix B, are now 
summarized. 

For random NRZ data with equi probable symbols, the average proba- 
bility of error conditioned on the misalignment (XT) of the symbol syn- 
chronization clock relative to its nominal position is given by 

Pe(x) = ^ + Jjerfc (1 - „ - 2|»| ) 

+ + 2|x|)] +i^e>-fc j^/yw^O - 2n)] ; 

0 < |X| < i , (6) 






# 


# 




where either (1) or (2) still applies as the definition of asymmetry. 
Note that, when x=o, (6) reduces to (3). Figure 2 illustrates the 
SNR degradation at a symbol E_/Np, of 1.5 dB (corresponding to a decoded 
bit error probability of 10”^) as a function of X with percent asymmetry 
nxlOO as a parameter. It is observed that, in the neighborhood of the 
nominal symbol synchronization lockup point (x=0), the sensitivity of 
SNR degradation due to asymmetry is extremely small, even for large 
asymmetry values. This makes the results given in Section 2.1.2 rela- 
tively insensitive to the assumption of a hard-locked symbol synchroni- 
zation, even though this assumption yields the minimum SNR degradation 
due to asymmetry. 

Another approach (not discussed in Appendix B) to assessing SNR 
degradation due to both asymmetry and symbol synchronization misalign- 
ment is to assume a model probability density function p(x) (p.d.f.) 
for X and average the conditional error probability P^(x) over this 
p.d.f. to obtain the average error probability performance [2]. In 

this regard, we postulate a Tikhonov p.d.f. for p(x) which is entirely 

2 

characterized in terms of the variance a of the synchronization error. 

A 

Thus, for NRZ data, we have [3], 
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Figure 2. SNR Degradation Due to Data Asytninetry as a Function of Symbol Sync 

Clock Misalignment 



In order to average the conditional error probability over the p.d.f. 
in (7), we must, in addition to (6), characterize P^(a) over the inter- 
val n/2 <. |x| j< 1/2. Following the procedure given in Appendix B, it 
is straightforward to show that (also see [2]) 


P^(X) = erfc ^ erfc 


IjE^/Ng (l-n-2[Ai) 


+ 7T erfc 

O 


Je^/Nq (l-n+2|A.|) ; n/2 < | a| < 1/2 . (8) 


Thus, using (6), (7) and (8), the average error probability can be 
computed from 

fl/2 

Pc = P.(A)p(A)dA . (9) 

^- 1/2 ^ , 

Figures 3 and 4 illustrate P^ [as computed from (9)] versus E^/Nq in 
dB with as a parameter for asymmetry values of 3 % and 6%, respec- 
tively. (These numerical values are taken from [2].) Again at 
Es/Nq= 1.5 dB, the additional degradation due to symbol synchronization 
jitter is on the order of a few tenths of a dB, even for a timing error 
(as measured by its standard deviation a^) as large as 7 %. 

2.1.4 Error Probability Performance for Asymmetric 

Manchester Coded Symbols 

Although Manchester coding cannot be employed on the high rate 
channel because of bandwidth limitations on the Ku-band return link, 
such coding is used on the lower rate channels, and thus the effect 
of asymmetry on the performance of these channels is potentially of 
interest. 

When Manchester coding is employed, then relative to the NRZ 
sequence, the Manchester coded waveform has 3/2 as many transitions. 

Thus, since SNR degradation due to asymmetry is directly related to 
the average transition density of the data sequence, one would intu- 
itively expect that, for the same asymmetry, the Manchester coded case 
should yield a larger SNR degradation than the corresponding NRZ case. 
Whether or not and to what extent the above intuitive notion is indeed 
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true depends on how one defines percent asymmetry for the Manchester 
coded case. 

To demonstrate this point, consider the NRZ sequence of Figure 1 
and the corresponding asymmetric Manchester waveforms which, for the 
two previous asymmetry models, are illustrated in Figure 5. Here, 
corresponding to Model 1, A/2 denotes the fractional (relative to the 
half-symbol time T/2) elongation of the positive half pulse in the 
Manchester coded waveform, and for Model 2, 6 is accordingly defined 
relative to the same half-pulse width. Then, defining asymmetry as 
was done for the NRZ case but now relative to the half-pulse duration 
gives (for Model 1) 


^ Id +4/2) - 7(1 - A/2) 

Asymmetry = n = 4- f 

i(l+A/2) +^(l-A/2) 


and, for Model 2, 

1(1+26) - J(1 - 26) 

T) = j f — = 26 , 

|(1 + 26 ) + ^(1 - 26 ) 


( 10 ) 


( 11 ) 


which are identical with (1) and (2).* Once again, as in the NRZ case, 
the in-phase integrate-and-dunip output depends, in general, on the 
polarity of the symbol over which it is integrating and the preceding 
and succeeding symbols. Thus, evaluating this integrate-and-dump 
output for the eight possible three-symbol sequences and noting that, 
for Model 1, the nominal bit synchronization lock-up misalignment is 
now AT/8, we get an expression for the average probability of error 
corresponding to asymmetric Manchester coded random data, namely. 


erfc 


/yNgd-n) 


+ ^ erfc 


(1 - f) 


( 12 ) 


* 

Note that, while this definition allows for equal percent 
asymmetry when compared with NRZ, the actual time asymmetry (as mea- 
sured in seconds) is different for the two cases. 




# # # ^ 










16 


Using (12) and (3), Figure 6 illustrates the SNR degradation in dB 
versus percent asymmetry for Manchester and, for comparison, NRZ coded 
data at E^/Nq= 1.5 dB. Here note that the Manchester code yields a 
larger SNR degradation than NRZ for small asymmetry values, while the 
reverse is true for large percent asymmetries. 

If, on the other hand, we define data asymmetry relative to the 
symbol time T for both NRZ and Manchester data, then equal percent 
asymmetry implies equal amounts of asymmetry (in seconds) as measured 
by the actual time displacements of the waveform transitions. Thus, 
for Manchester coded data, we have that n=A/4 for Model 1 and ^= 6 
for Model 2. The corresponding expression for error probability is 
identical to (12) with n replaced by 2ri. Using this definition. 

Figure 7 plots SNR degradation in dB versus percent asymmetry for 
Es/Nq= 1.5 dB for Manchester coded data and the corresponding results 
for NRZ data obtained from (3). Note that now the Manchester coded 
case always yields a larger SNR degradation for a given percent asym- 
metry. The conclusion to be reached here is that, in making comparisons 
between asymmetric NRZ and Manchester coded systems, one must exercise 
care in selecting a definition which is appropriate to the particular 
application at hand. 

To generalize the results given in (12) to the case of imperfect 
symbol synchronization, one merely follows the procedure outlined in 
Appendix B for NRZ data, whereby (6) and (8) are obtained as extensions 
of (3). In particular, for a timing misalignment of AT/2 (consistent 
with our first definition of asymmetry relative to the half-pulse width 
and the fact that the time base for the symbol synchronizer is now T/2), 
the conditional error probability is given by (also see [2]) 


erfc 




+ ^ erfc 




Pe(») 


+ -g erfc 


+ IM) 


0 <. I A I 


erfc 


JE^/N„{1-2|a|) 


+ -g erfc 
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+ g- erfc 


^E^/Nod+f - dl) 
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Figure 6. SNR Degradation vs. Percent Asyninietry for Manchester and NRZ 

Coded Data 
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Substituting (7) and (12) into (9) results in the average error proba- 
bility performance for Manchester coded data with asymmetry. The cor- 
responding performance plots analogous to Figures 3 and 4 for NRZ data 
are given in Figures 8 and 9. Comparing the NRZ and Manchester data 
results, we observe that, in terms of error probability performance, 
Manchester data appears to be less sensitive to symbol synchronization 
error. Of course, it should be pointed out that this conclusion is 
highly dependent on the time base used as normalization for the asym- 
metry definition. In fact, if the alternate definition is used, wherein 
both NRZ and Manchester asymmetries are defined as fractions of the 
symbol time, then for fixed values of n and a , the error probability 

A 

performance of Manchester data will always be worse than that for 
NRZ [2]. 

2.1.5 The Effect of D.C. Restoration on Error Probability 

Performance in the Presence of Asymmetry 

The investigation of the effects of d.c. restoration* on communi- 
cation link performance was prompted by test results [4] conducted in 
the Electronic Systems Test Laboratory (ESTL) at JSC which indicated 
significantly less performance degradation than that predicted by the 
analytic results of the previous sections. In particular, it was found 
that d.c. restoration tends to reduce the degrading efretts of data 
asymmetry and thus it was necessary to incorporate d.c. restoration into 
the analytical model in order that the predicted performance better 
agree with the test results. 

The effect of d.c. restoration on data detection is most easily 
accounted for by artificially shifting the decision threshold (nominally 
at zero) against which the matched filter output is compared. The 
amount of this artificial shift in threshold depends upon the specific 
way in which d.c. restoration comes about. 

The simplest method of achieving d.c. restoration is to capacitive- 
couple the input signal to the symbol synchronizer. In this case, the 
artificial threshold shift equals the d.c-. component of the asymmetric 


D.C. restoration refers to the process by which the d.c. value 
of the asymmetric data waveform is forced to zero. 
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data waveform in front of the capacitor which, for random data with 
transition density D, is easily shown to be [4] 


\ ~ nD , 


(14) 


where «^E^/T is the data pulse amplitude in Figure 1. Computing the 
matched filter output for the eight possible three-symbol sequences 
made up of the present, preceding, and succeeding symbols and shifting 
these outputs by gives the result for the error probability per- 
formance of as 3 mimetric NRZ data with d.c. restoration by capacitive 
coupling, namely. 


= erfc 


UE./No(l-nD) 


+ Tc erfc 
1 6 


^E^/Nq (1 - 2n+ nD) 


+ ^ erfc 


jEs/N(jO-n+nD) 


+ ^ erfc 


JE^/No (1 + nD) 


. (15) 


For equiprobable data symbols (D=0.5), (15) simplifies to 


erfc 


r 




+ ■ erfc 




+ jg- erfc 




(16) 


Comparing (16) with (3), we observe that the effect of d.c. restoration 
is to compensate for the data asymmetry by shifting the effective 
decision threshold away from the shortened symbol.s. 

For a given value of asymmetry, the value of E^/Nq required to 
obtain [as computed from (16)] equal to 10"^ can be calculated. 
Comparing this value of E^/Nq With that obtained from (5) for the same 
Pg gives the SNR degradation for asymmetric NRZ data with d.c. restora- 
tion by capacitive coupling. Figure TO illustrates this SNR degradation 
versus asymmetry along with the comparable results obtained from (3) 
corresponding to no d.c. restoration (direct coupling). 

Another method of achieving d.c. restoration, which depends 
specifically on the symbol synchronizer implementation itself, is to 
require the matched filter output to have zero crossings at the center 
of each symbol period that starts with a data transition. In this case, 
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Figure 10. Performance Degradation for Capacitively Coupled Matched Filters; 

Random Data (D = 0.5) 



the effective shift in decision threshold relative to its nominal (zero) 
value is [4] 

\ = n ^€71 . (17) 

Comparing (17) with (14), we can immediately conclude that the error 
probability for this method of d.c. restoration is given by (15) with 
D = 1 , i .e. , 

08 ) 
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Again by determining those values of E^/Nq required to obtain P^=10 
for various values of asymmetry, one can compute the SNR degradation for 
d.c. restoration based on symbol timing. The results of these calcula- 
tions are illustrated in Figure 11, along with experimental test results 
taken in the ESTL for the sake of comparison. It is to be noted that 
the experimental results include the effect of bandl imiting, whereas 
the theory as predicted by (18) in no way accounts for this effect. 
Furthermore, the data detector used in the experimental setup is not 
a true matched filter as is assumed for the analytical model. Surpris- 
ingly enough, however, the analytic and experimental results show rea- 
sonably good agreement. In the next section, we consider the combined 
effects of bandlimiting and asymmetry on the performance of a filter- 
sample type data detector which is a more realistic model of the detector 
used in the ESTL experimental tests. 

2.1.6 The Combined Effects of Bandlimiting and Data Asymmetry 
on Error Probability Performance 

In the previous sections, the degradation of error probability 
performance due to data asymmetry alone was considered. Since, in 
reality, a practical receiver operates in the presence of both band- 
limiting and data asymmetry, it is important to study the combined 
effects of these two sources of degradation on its error probability 
performance. Assuming a filter-sample type detector (see Figure 12) 



SHR Degradation (dB) 






Figure 12. Band Limit and Sample Detector for Binary Data With Asymmetry 





and NRZ data. Appendix C investigates in detail the combined degrading 
effects of intersymbol interference and SNR reduction of the desired 
signal due to bandl imiting by an ideal low-pass filter and data asym- 
metry. The results of this analysis, which are first derived assuming 
no d.c. restoration and later modified to include this effect, are sum- 
marized below. 

Figure 13 illustrates the bit error probability as a function 

of energy- to-noise ratio with BT as a parameter and 10% asymmetry while 

Figure 14 illustrates as a function of asymmetry the required energy-to- 

-5 

noise ratio to achieve a bit error probability of 10 with BT still the 
parameter. Both of these figures assume no d.c. restoration. The cor- 
responding results which include d.c. restoration by capacitive coupling 
are illustrated, respectively, in Figures 15 and 16. It is concluded 
that, for small values of asymmetry (less than about 10%), the optimum 
filter bandwidth - bit time product remains equal to 0.9, and the corre- 
sponding amount of energy-to-noise ratio degradation at a fixed error 
rate is virtually independent of whether or not d.c. restoration is 
present. In particular. Tables 2 and 3 summarize typical results for 
the two cases when BT=0.9. The quantity aE^/Nq represents the additional 
E$/No required at a given value of data asymmetry relative to its value 
at zero asymmetry. 

On the other hand, for values of BT other than the optimum, d.c. 
restoration has a beneficial effect in reducing energy-to-noise ratio 
degradation due to data asymmetry. 

The other conclusion which can be reached by comparing the above 
numerical results with those in the previous section is that the sample 
detector is much less sensitive to data asymmetry than the integrate- - 
and-dump detector. This is not surprising when one realizes that the 
degradation due to data asymmetry for the integrate-and-dump detector 
comes about because of a reduction in the signal energy as a result of 
integrating only over a fraction of the total bit interval. By comparison 
the reduction in the peak of the filter response at the sampling instant 


* 

A previous study [5] of the degrading effect of ideal low-pass 
bandlimiting alone on the bit error probability performance of a 
filter-sample detector revealed that BT=0.9 was optimum from the 
standpoint of minimum degradation. 
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Figure 16. Energy- to-Noise Ratio as a Function of Data Asymmetry for Fixed Bit Error Probability and 
Filter Bandwidth - Bit Time Product as a Parameter; D.C. Restoration by Capacitive Coupling 





Table 2. No D.C. Restoration 


= 10 ' 


E^/Nq (dB) aE^/Nq (dB) 


9.97 

9.98 
10.00 
10.06 
10.13 
10.31 


Pj = 10-® 

n (%) 

E^/Ng (dB) 

dEj/Ng (dB) 

0 

11.28 

0 

2.5 

11,29 

0.01 

5.0 

11.32 

0.04 

7.5 

11.37 

0.09 

10.0 

11.43 

0.15 

15.0 

11.63 

0.35 


Table 3. D.C. Restoration by Capacitive Coupling 

Pe " p£ = 10"^ 

E^/Nq (dB) aE^/Nq (dB) n {%) E^/Nq (dB) . 

9.97 0 0 11,28 

9-99 0.02 2.5 11.30 

10.02 0.05 5.0 11,32 

10.12 0.15 10.0 11.40 

10.25 0.28 15.0 11.54 

10.44 0.47 ?n.n n 7 ^ 


= 10"^ 

n {%) 

E3/Ng (dB) 

AE^/Nq (dB) 

0 

11.28 

0 

2.5 

11.30 

0.02 

5.0 

11.32 

0.04 

10.0 

11.40 

0.12 

15.0 

11.54 

0.26 

20.0 

11.73 

0.45 





clue to data asymmetry should have a considerably lesser effect. In fact, 
in the absence of bandlimiting, the integrate-and-dump detector would 
still yield considerable E^/Nq degradation due to data asymmetry, whereas 
the sample detector would show none. Of course, the sample detector 
without bandlimiting would produce infinite noise power and, thus, on 
an absolute E^/Nq basis would be far worse than the integrate-and-dump. 

2.1.7 Performance of a Gated Integrate-and-Dump Filter 

for Detection of Asymmetric Data 

This section deals with still another implementation of a data 
detector for asymmetric data, namely, the gated integrate-and-dump 
filter (see Figure 17). The motivation for studying the performance 
of such a detector stems from several considerations. First, from an 
implementation point of view, the gated integrate-and-dump has the 
advantage that operation at high data rates can be accomplished with 
smaller circuit losses than the ideal integrate-and-dump since the 
constraint on its switching times at the symbol transition instants 
can now be considerably relaxed. Second, since. the input data stream 
possesses asymmetry, the ideal integrate-and-dump is no longer the opti- 
mum detector, and it is thus possible that an alternate (possibly simpler 
to implement) detector could yield superio.r performance. 

Appendix D discusses in detail the performance of the gated 
integrate-and-dump filter as a function of data asymmetry and gating 
interval and demonstrates the relation between these parameters which 
optimizes this performance. Only the case of NRZ data is treated; how- 
ever, as before, the results are obtained for both the case of no d.c. 
restoration and the case of d.c. restoration by capacitive coupling. 

These results are summarized as follows. 

In the absence of d.c. restoration, the average error probability 
performance of the gated integrate-and-dump is given by 



^ erfc p + i erfc [^E^/N^ 


+ ^ erfc 


/I - 2s /J 
0 <. e ^ n/2 


1 erfc pEj/N(,(l - 2e)‘] ; 


n/2 e <_ 1/2 


(19) 
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Figure 18 is a plot of Pr versus e with n as a parameter and 

E^/No“9.6 dB (corresponding to P^=10” when e=n=0). We observe 

from this figure that, for a given value of data asymmetry n, P^ is 

minimized by choosing e-n/2. Figure 19 is an illustration of the 

symbol energy-to-noise ratio (in dB) required to achieve an average 
-5 

error rate of 10 in the presence of data asymmetry, The curve labeled 
e=0 corresponds to the performance of the ideal integrate-and-dump 
(see also Figure 3 of Appendix A), The remaining curves indicate a 
constant for values of n:l2e in accordance with the second equa- 

tion of (19) followed by an increase in E^/Nq with n as required by 
the first equation of (19). Note that each of these curves cross the 
e=0 curve at some value of n» say r;Q» which means that, for n>riQ» 
the gated integrate-and-dump outperforms the ideal integrate-and-dump 
in the sense of requiring less E^/Nq for a given average probability 
of error. The dashed curve in Figure 19 represents the performance 
corresponding to selecting e = n/2 at each value of n and is thus the 
best achievable with the gated integrate-and-dump. 

When d.c. restoration by capacitive coupling is present, then 
the analogous result to (19) is 



n/2 < e < 1/2 . (20) 


Figure 20 again illustrates P^ versus g With n as a parameter 
and E-^/Ng= 9.6 dB, where P^ is now computed from (20). We observe from 
this figure that, for a given value of data asymmetry, there exists a 
value of e which minimizes Pj^; however, unlike Figure 18, this value 
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Figure 18. Average F.rror Probability versus Gate Interval at Symbol Edge with 
Data Asymmetry as a Parameter - No D.C. Restoration 


symbol rperay-to-Noise Ratio versus Data 
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of e, namely, is not equal to n/2. Figure 21 is the analogous 

figure to Figure 19 when d.c. restoration is present. Again, the dashed 
curve corresponds to which represents the best achievable per- 

formance using the gated integrate-and-dump as a data detector. Com- 
paring Figure 21 with Figure 19, we observe the considerable reduction 
in SNR degradation due to data asymmetry when d.c. restoration is employed. 
This improvmenet is analogous to that achieved when other types of data 
detectors are used. 

In conclusion, the use of a gated integrate-and-dump filter for 
detection of asymmetric NRZ data can, depending upon the amount of data 
asymmetry present, produce significant improvement in SNR degradation 
due to asymmetry relative to that of an ideal integrate-and-dump filter. 


2.2 Antijam Techniques for Shuttle Communication 

There are three basic techniques for antijam communication. 

These techniques are direct sequence pseudonoise (DSPN), frequency 
hopping (FH), and time hopping (TH), and any hybrid of these techniques. 


2.2.1 Direct Sequence Pseudonoise Spreading 

The most straightforward way to widen the spectrum of a data 
signal is to multiply (modulate) it by a wideband signal. . Such a 
spreading signal must have correlation properties that aid in acqui- 
sition and tracking. One of the best signals fitting these specifica- 
tions is the binary pseudonoise (PN) signal. This type of spread 
spectrum is used on the S-band and Ku-band Shuttle forward links to 
provide the low flux density over a 4 kHz bandwidth. 

The antijam performance of a DSPN system is easily computed. 
Assume that the additive interfering signal has noise power N within 
the PN bandwidth. No matter what its initial bandwidth, when the 
received signal is demodulated, this interference is spread over the 
PN bandwidth. Note that the ratio of required energy per bit to 
single-sided noise spectral density is given by 



0/req 


NT, 


( 21 ) 




Figure 21. Symbol Energy-to-Noise Ratio. versus Data Asymmetry with Gate 
Interval at Symbol Edge as a Parameter - D.C. Restoration by 

Capacitive Coupling 
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where S is the signal power, N is the noise power in the PN bandwidth, 
T|^ is the bit duration, and is the PN chip duration. 

For a tone jammer of the carrier frequency with power J and 
white noise power in the PN bandwidth equal to 

N, 


Therefore, 



SL 


N = + J . 
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ST, 
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Let the available Ej^/Nq = STj^/Nq be denoted (Eb/NQ)^; then 



J/Tc 
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1 


S \T 


<V«0)req ■ ■ 


(24) 


Define the circuit margin in the presence of white noise only equal 
to M; that is. 


Rewriting (24), 


or 


M = 






req 



= M - 1 


T, 


(M-1) 


(25) 
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Note that 1/Tj^ is equal to the bit rate R|^ and 1/T^ is equal to the 

(yNo),(j/s) 

(fToT • 


PN chip rate R^. Therefore 


( 28 ) 


As an example of the use of (28), consider the S-band forward link. 
Let the required E|j/Nq for a bit error rate of 10"^ be 



rate 1/3, K=7 
convolutional 
coding required 
for Pj^ = 10-4 


1.5 dB + 

1.0 

+ 1.0 

loss in 

loss in PN 

TORS 

bit sync 

despreading 

loss 


= 7 dB . 


(29) 


Let the available (Ejj/Nq)^ for Rjj=32 kbps be 



and let J/S be 20 dB; then 

(J/S) (E. /N.)^ 

= (32 kbps)— -ppJy-yj- = 78 Mbps . (31) 

Note that, since the energy is normalized to the bit rate rather than 
the coded symbol rate (3 x 32 = 96 ksps), the bit rate is used in (28). 
Also note that (28) assumes a PN code length of infinite period but, 
in actual practice, the PN code rate must be greater than Tj^/T^. 
Therefore, in this case, the code length must be 4095 or greater. 

This example illustrates the disadvantage of using DSPN to achieve 
antijam communication for moderate data rates and small circuit margins 
in white noise. The main disadvantage is that extremely high chip 
rates are required and, as a consequence, the synchronization time 
becomes unacceptably long. 

2.2.2 Frequency Hopping 

Another method of anti jam modulation is frequency hopping (FH). 
Using this technique, the carrier is hopped over a large number of 
frequencies (i.e.. Targe bandwidth). The carrier frequency selection 
process is controlled by a PN sequence which could be nonlinear to 


:d 


43 








3 


3 


D 




deter determining the selection code by the jammer. Since the frequency 
hopping rate is of the order of the bit rate and much less than the 
total hopping bandwidth, the acquisition of a FH system is much easier 
than a PN system with similar bandwidth. 

The antijam performance of a FH system depends on the type of 
jamming. For example, if the jammer placed equal power at each of 
the possible carrier frequencies, then the FH processing gain over a 
non-frequency hopping system is 



number of carrier frequencies 



(32) 


where B is the total available bandwidth and R is the larger of the 
frequency hop rate or the bit rate. If, instead of jamming each pos- 
sible carrier frequency, the jammer could jam 0.2 of the carrier 
frequencies with enough power to cause an error (probability 1/2), 
then the bit error rate would be 10%. Thus, to protect against par- 
tial band jamming, an error correcting code must be employed; otherwise, 
an error would be made every time a carrier frequency from the jammed 
band was selected. 

It is important to note that only noncoherent frequency hopping 
is being considered, since large hopping bandwidths make it impractical 
to maintain phase coherence across the total bandwidth. Therefore, 
each hop has an independent phase but the phase is constant over the 
hop. Because of the noncoherent frequency hopping, unless a very slow 
frequency hopping rate is used which would allow frequency following 
by the jammer, it is impractical to lock up a carrier tracking loop at 
each frequency hop. Therefore, in the Shuttle application, a modifica- 
tion to the carrier recovery must be made to the demodulator. Probably 
the simplest modification is to employ differential PSK (DPSK) and 
derive the carrier phase by comparing the successive bits. 

To calculate the anti jam performance of a FH system, effective 
signal-to-noise ratio is defined as 



0 


9 


(33) 


where B is the hopping bandwidth, is the bit rate, J is the interfer- 
ence power, and S is the signal power. Note that, in using 
it is assumed that the thermal noise is insignificant to the jammer 
power (i.e., 10 dB less than J). For DPSK modulation with soft deci- 
sion Viterbi decoding (rate 1/2 or rate 1/3), the required 
is shown in Table 4. 

For FH-PN/DPSK, the data modulation is first spread with a PN 
code as in DSPN and then the carrier is hopped. Using (33), the 
required bandwidth can be determined for the example used in the DSPN 
case. For rate 1/3 Viterbi decoding, = 7.8 dB at Pj^=10“^. 

Thus, 



with = kbps for the S-band and (J/S) = 20 dB, then 

B = 19.6 MHz . (35) 

For the Ku-band forward' link, is equal to 216 kbps. With no 
coding on the Ku-band forward link, = 42.6 dB for P|^=10”^. 

To achieve (J/S) = 20 dB, then 

B = 393 GHz , (36) 

which is clearly impractical. However, if rate 1/3 Viterbi decoding 
is used, then (E^fj/f'^g^eff ~ B=130 MHz. 

Note again that (33) assumes a reasonably high circuit margin 
in thermal noise. If this cannot be achieved, then a finer anti jam 
performance analysis can be made. But even if B was required to be 
greater than 130 MHz, the synchronization requirements would not be 
increased over the modest requirements of the FH system in the example 
because the synchronization time is independent of the hopping bandwidth 

2.2.3 Time Hopping 

Whereas FH channelizes via different carrier frequencies within 
a given bandwidth, time hopping (TH) channelizes via time slots within 
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a given time frame. A burst transmission is employed within the slot 
M ^ selected by a PN code or nonlinear sequence. Provided the TH trans- 

mitter can operate at a high peak power but with low duty factor, the 
TH system has the same antijam performance as the FH system. The TH 
system does not appear practical for the Shuttle system unless the 
Q transmitter is converted to a high peak power pulsed operation. 

As a conclusion to the antijam techniques for the Shuttle 
communications, it appears that a frequency hopping system could 

\ provide the best antijam performance for the Shuttle S-band and Ku-band 

I 

I Q forward links with some modification to the demodulator in deriving 

I carrier phase, 



3.0 TRACKING TECHNIQUES ANALYSIS 

As a part of the overall effort to assess the Ku-band system 
design and to evaluate possible system simplifications/modifications 
that may be achievable, this task requires Axiomatix to study the design 
and implementation of tracking techniques which are suitable for incor- 
poration into the Orbiter Ku-band communication system. Emphasis is 
placed on maximizing tracking accuracy and communications system flexi- 
bility while minimizing cost, weight, and system complexity of Orbiter 
and ground system hardware. 

In a host of previous studies [6-13], much consideration has been 
given to implementation of the three-channel Orbiter modulator, whose 
purpose is to generate a signal for simultaneous transmission of three 
channels of information on the Ku-band return link. One of the possible 
structures for Mode 1 is referred to as the Three-Channel Quadrature 
Multiplex Modulator [9], wherein the high data rate channel (50 Mbps) 
is biphase modulated on an in-phase carrier and the two lower rate chan- 
nels (192 kbps and up to 2 Mbps), after being biphase modulated onto 
quadrature squarewave subcarriers and summed, are amplitude modulated 
onto the quadrature carrier (see Figure 22). The receiver for the 
three-channel quadrature multiplex signal employs biphase Costas loops 
for carrier and subcarrier recovery. While the performance of the carrier 
recovery loop has previously been investigated [6, Appendix C], the com- 
panion performance analysis for the subcarrier loop has not. Thus, the 
first part of this task focuses on developing the theory necessary to 
carry out this analysis. In particular, we begin with a fundamental 
study of the tracking performance of biphase Costas loops when subjected 
to an unbalanced QPSK input. The plural usage in the phrase "biphase 
Costas loops" refers to the fact that we shall study several different 
Costas loop implementations, each of which has advantages over the others, 
depending on the application. Following these more general analyses, we 
shall apply the results specifically to the subcarrier tracking loop of 
the Ku-band return link. 

ORIGINAL PAGE IS 
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3.1 Tracking Performance of Unbalanced QPSK Demodulators 

Unbalanced quadriphase-shift-keying (QPSK) used on the Ku-band 
return link is an attractive means for transmitting two digital data 
streams which have different average powers. The two data streams are 
not constrained to have identical data rates nor must they have the 
same data format; e.g., one might be an NRZ sequence and the other a 
Manchester code. In fact, it is the difference in data rates which 
causes the unbalance of power when it is desired to have symbol energies 
and therefore error rates on the two channels within the same order of 
magnitude. 

3.1.1 Biphase Costas Loop With Passive Arm Filters 

Previous results [14,15] have indicated that, when the unbalanced 
power ratio is large, e.g., approximately 4:1 or greater, a biphase 
Costas loop is a more efficient demodulator than a fourth-power tracking 
loop. These results, however, accounted only for the filtering effect 
produced by the two passive arm filters of the loop on the equivalent 
additive noise perturbing the' loop. When the bandwidth of these filters 
is selected on the order of the data rate, as is typical of optimum 
Costas loop design for tracking purposes [16,17], the filtering degra- 
dations of the data modulations themselves and the cross-modulation noise 
produced by their multiplication in the loop often cannot be neglected. 

An analysis that incorporates these additional filtering effects 
into the performance characterization of a biphase Costas loop demodu- 
lator of unbalanced QPSK is presented in Appendix E. Many of the results 
obtained in this appendix are in the form of closed-form expressions 
which can easily be evaluated numerically for design and performance 
prediction purposes. The generality of the results enables them to be 
applied to a wide variety of applications, in particular, the performance 
of the subcarrier tracking loop for the three-channel Orbiter Ku-band 
return link, which will be discussed later on. For the moment, we shall 
briefly summarize the results given in Appendix E. 

It is common practice to characterize the tracking performance 
of a Costas loop by the variance of the loop phase error 2 (jp , often 
referred to as the tracking phase jitter. In the linear region of 
operation, an expression for this quantity is given by 
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where P = Pj/^qB^ is the total power (Pj) to noise spectral density (Nq) 
ratio in the loop bandwidth (Bj^) and is the so-called squaring loss , 
which represents the additional degradation in equivalent loop signal- 
to-noise ratio relative to a linear loop, i.e., one in which there 
exists no nolsex noise or signal x noise error signal components. Thus, 
from (37), we see that characterization of the loop squaring loss in 
terms of system parameters is sufficient for predicting the loop's 
tracking performance. 

As a numerical illustration of this characterization, consider 
the case where the low rate modulation m-j(t) is Manchester coded data, 
the high rate modulation m^(t) is NRZ data, and the arm filters are 
single-pole Butterv/orth (RC) filters. Then Figure 23 illustrates (for 
fixed fractional channel powers n-j and r\^ corresponding to a 4:1 power 
ratio on the two channels) the behavior of as a function of the ratio 
of two-sided arm filter noise bandwidth to the higher of the two data 
rates i?2 = ratio of data rates and Pi-T2^^0 Pentam- 

eters. Assuming Pj/Ng to be fixed, then the variation of squaring loss 
with PjT 2 /Nq directly reflects the effect of changing the high data rate 
i? 2 ' Furthermore, at low values of B^/i? 2 , we observe from Figure 23 
that additional interesting peaks and valleys of the squaring loss char- 
acteristic occur. These extrema represent tradeoffs between SxS distor- 
tion and crpss-modulation noise or S^N power, depending on which of the 
latter dominates the total noise. 

The corresponding numerical evaluation of the tracking jitter, 
for a fixed ratio of arm filter noise bandwidth to the loop noise band- 
width {B^./B|^) is shown in Figure 24. The minimum values of a^p for some 
ratios of P^P'^^s^nt best design points when the combined effect 

of NxN distortion and cross-modulation noise or S^N power is minimal. 
Assuming Pj/Ng to be fixed, the variation of a^p with PjT 2 /Ng is shown 
in Figure 25. As is intuitively true, the tracking jitter performance 
improves with the increase of PjT 2 /Ng. 
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Figure 24. Tracking Jitter Standard Deviation versus Ratio of Arm Filter Band 
width to High Data Rate; PtT 2/I^0 and /? 2 //?-j are parameters; m^(t) i 
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3.1.2 Biphase Costas Loop With Active Arm Filters 


In a previous Costas loop study for biphase modulation [16], it 
was demonstrated that considerable improvement in tracking performance 
could be obtained by employing active arm filters of the integrate-and- 
dump type as opposed to passive arm filters:, An investigation of whether 
a similar performance improvement can be obtained for an unbalanced QPSK 
modulation is presented in detail in Appendix F. The highlights of this 
investigation are summarized as follows. 

In Section 3.1.1, we observed that, for given values of the data 
rate ratio power ratio P-|/P 2 ’ total signal-to-noise ratio 

in the high data rate bandwidth PjT 2 /Nq, an optimum arm filter band- 


width or, equivalently, exists in the sense of maximizing 

(i.e., minimizing the squaring loss). Using that value of B^/i? 2 > namely, 

(B./i?o)„_., and defining the corresponding value of S. by , then 
1 c ope L opt 

the minimum improvement in tracking performance (or, equivalently, in 

squaring loss performance) obtained by employing active integrate-and- 

dump arm filters as opposed to active arm filters is given by 

1 = S./S, , where S. is used here to denote the squaring loss of 

L "-opt L 

the Costas loop with integrate-and-dump arm filters. 

Closed-form expressions for are derived in Appendix F for all 
combinations of NRZ and Manchester data formats for the two channels 
and both synchronized and unsynchronized symbol clocks. As an example, 
assuming single-pole (RC) arm filters as the basis of comparison for 
the Costas loop with passive arm filters and assuming unsynchronized 
symbol clocks for the Costas loop with integrate-and-dump filters, 

Figures 26 and 27 illustrate I (in dB) versus the channel power ratio 
P-j/P^ with the data rate ratio ^ parameter and values of total 

power-to-noise ratio PjT 2 /Nq typical of coded and uncoded systems. We 
observe from these figures that the improvement in squaring loss per- 
formance of using integrate-and-dump filters as opposed to single-pole 
arm filters is an increasing function of P 1 /P 2 depends heavily on 
the choice of data formats for m-j(t) and m 2 (t). 

Similarly, using many numerical illustrations, it is shown in 
Appendix F that, for a fixed ratio of data rates and total power-to- 
noise ratio in the higher data bandwidth, the squaring loss itself 


I (dB) 



Figure 26. Squaring Loss Improvement {in dB) Using Integrate-and-Dump Arm ’ 
Filters as Opposed to Single-Pole Passive Arm Filters; m,(t) is 
Manchester, m 2 (t) is NRZ I 
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increases as the ratio of powers in the two channels increases, and the 
rate at which this loss increases (tracking performance deteriorates) 
also depends heavily upon the data formats in each channel . Thus, it 
is concluded that, when the ratio of data rates is on the same order 
of magnitude as the inverse of the power ratio, i.e., approximately 
equal signal energies in the two channels, the biphase Costas loop can 
be used as an efficient demodulator of QPSK. On the other hand, if the 
energy in the two channels is very unbalanced, e.g., one channel is 
coded and one is uncoded, then it is still possible to efficiently use 
a biphase Costas loop for demodulation of unbalanced QPSK, provided that 
the higher data rate channel is Manchester coded. It is understood that 
the foregoing conclusions are quite general and are not intended to rule 
out specific design situations in which sufficient total power- to-noise 
ratio is available to tolerate large squaring losses. In an individual 
situation, one must resort to the specific numerical results given in 
the illustrations to determine the suitability of employing a biphase 
Costas loop for demodulation of unbalanced QPSK. 

3.1.3 Biphase Polari ty-Type Costas Loop With Passive Arm Filters 

In Sections 3.1.1 and 3.1.2, we addressed the problem of tracking 
an unbalanced QPSK signal with a conventional biphase Costas loop with 
analog input phase detectors and an analog third multiplier (the one 
that forms the loop error. signal ). Because of dc offsets associated 
with analog multipliers, it is common practice to hard-limit the 
in-phase* channel arm filter output and replace the analog third multi- 
plier with a chopper-type device (switched multiplier) which typically 
exhibits much less offset (see Figure 28). While it is also possible to 
replace the input in-phase and quadrature analog phase detectors with 
switched multipliers, the impact of doing so on the resulting tracking 

performance is mininsal since the arm filters will pass only the first 

2 

harmonic of these phase detector outputs. Thus, aside from the 8/ir 
power loss associated with the first harmonic of a square wave, the 


For unbalanced quadriphase, we shal larbitrarily refer to the 
in-phase channel as that corresponding to the point of data extraction 
for the higher power signal. 









performance of the loop would be identical to that given in Section 
3.1.1 for an analog third multiplier or that to be presented here for 
a switched third multiplier. For ease of terminology, we shall refer 
to a conventional biphase Costas loop with a switched third multiplier 
as a "biphase poTari ty-type Costas loop" or, even simpler, a "polarity- 
type Costas loop." 

Generally speaking, introduction of a limiter (hard or soft) 
into a system results in signal suppression, the amount of which is 
a function of the signal-to-noise ratio at the limiter input. This 
signal suppression, in turn, reduces the total loop gain and, as a 
consequence, the loop bandwidth. Another potential problem with the 
limiter under strong signal conditions is that it may increase the 
tendency of the loop to false lock. 

Appendix G discusses in detail the tracking behavior of the 
polarity-type Costas loop with unbalanced QPSK input and compares its 
performance with that of the conventional Costas loop discussed in 
Section 3.1.1 and Appendix E. In particular, for the case of single- 
pole Butterworth (RC) arm filters and a particular combination of NRZ 
and Manchester coded data on the two channels, the squaring loss 
(tracking jitter penalty relative to a linear loop) is evaluated and 
illustrated as a function of the ratio of arm filter bandwidth to 
higher data rate and total signal power-to-noise ratio in this higher 
data rate bandwidth. Also numerically illustrated is the corresponding 
mean-squared tracking jitter performance as a function of these same 
receiver parameters. A summary of these results is given as follows. 
Figure 29 illustrates the variations of versus with 

as a parameter. Superimposed on these curves (in dashed lines) are 
the corresponding results obtained from Figure 23 for the biphase 
Costas loop with passive arm filters. We observe from these numerical 
results that, for high signal-to-noise ratios, the hard-limited loop 
actually outperforms the conventional loop and, depending on the data 
rate ratio, the improvement (in terms of squaring loss) might be as 
high as 2.8 dB. Also, for a given signal-to-noise ratio and a given 
arm filter bandwidth to high data rate ratio, the squaring loss does 
not change significantly with data rates when the ratio of the data 




rates is high. This is particularly true for small values of arm filter 
bandwidth to high data rate ratio. A comparison with the dotted curves 
of Figure 29 reveals that the same is not true for the conventional 
Costas loop with passive arm filters. However, the fact that the 
polarity- type Costas loop produces an improvement in tracking perform- 
ance at high signal-to-noise ratios over the biphase Costas loop with 
passive arm filters is not surprising in view of similar results 
recently demonstrated for biphase modulation [18]. 

3.2 Performance of the Subcarrier Tracking Loop for the 

Three-Channel Orbiter Ku-Band Return Link 

Having examined the many ways in which a biphase Costas loop can 
track an unbalanced QPSK signal, we are now in a better position to 
analyze the tracking performance of the subcarrier loop for the three- 
channel Orbiter Ku-band return link. The primary difference between 
the subcarrier loop tracking behavior and that of the loops discussed 
in Section 3.1 is the fact that the in-phase and quadrature data modu- 
lations which are input to the subcarrier loop are biphase modulated 
onto quadrature square-wave subcarriers as opposed to sine-wave sub- 
carriers. A secondary difference is the fact that, since the output 
of the quadrature phase detector of the carrier tracking loop serves 
as the input to the subcarrier loop, the performance of the latter 
depends on the carrier tracking loop phase error. An analysis which 
takes both of these differences into account is given in Appendix H. 

For simplicity of numerical evaluation, it was found convenient to 
assume perfect carrier tracking. The additional degradation due to 
the phase tracking jitter of the carrier loop itself can be easily 
assessed from the results of a previous report [6, Appendix C] wherein 
the performance of a Costas loop for recovering the carrier from the 
three-channel quadrature multiplex signal was studied. 

Some specific results from Appendix H are summarized below. In 
addition to perfect carrier tracking which was already mentioned, the 
following assumptions were made. The subcarrier frequency is 8.5 MHz, 
the high rate modulation m-jCt) in Mode 1 is a rate 1/2, constraint 
length 7, convolutional code (NRZ format) with data rate R^~ 50 Mbps, 
and the arm filters in the subcarrier loop are one-pole Butterworth (RC) 
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Under these assumptions. Figures 30 and 31 illustrate the subcarrier 
tracking jitter behavior for a fixed ratio of arm filter noise bandwidth 
to loop noise bandwidth In particular, Figure 30 illustrates 

the case where m 2 (t) is NRZ data at 500 kbps and m^tt) is Manchester 
coded data at 192 kbps, while Figure 31 corresponds to the same param- 
eter values as Figure 30 with the exception that m 2 (t) is now also 
Manchester coded. In both figures, the rms tracking jitter, (in 
% radians), is plotted versus the ratio of two-sided arm filter noise 
bandwidth (B^) to the higher data rate (i? 2 ) with total power-to-noise 
ratio in the higher data rate bandwidth (PjT^VNq) as a parameter. It 
is observed that the changes in the subcarrier tracking jitter as a 
function of obvious when m 2 (t) and m^it) are both Man- 

chester codes than when m 2 (t) is NRZ and m^Ct) is Manchester. 

Upon establishing a subcarrier reference signal, the two lower 
rate modulations, m 2 (t) and m^Ct), can then be demodulated. The error 
probability performance associated with these two data demodulations 
is the subject of discussion in the next section. 

3.3 Error Probability Performance of Channels 2 and 3 

of the Three-Channel Ku-Band Return Link 

Evaluation of error probability performance of BPSK, QPSK, and 
offset (staggered) QPSK receivers has been extensively covered in the 
literature [19-24]. While the techniques used there are certainly 
applicable to demodulation of unbalanced QPSK as on the two tow,er rate 
channels of the three-channel signal, the complexity of the evaluation 
when the ratio of data rates in the two channels is large prompts one 
to look for a simpler calculation procedure. Indeed, such an approach 
is possible when the noisy reference loss* is small or, equivalently, 
the effective signal^^to-noise ratio in the tracking loop bandwidth is 
large, i.e., the loop operates in its so-called linear region. Making 
such an assumption for purposes of error probability performance evalu- 
ation is quite reasonable when one realizes that this very same assumption 


"Noisy reference loss", is defined here as the equivalent increase 
in signal power required to produce the same error probability as obtain- 
able in a perfectly synchronized system. 
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has already been implied in assessing the tracking performance of the 

subcarrier loop. The approach taken is to expand the error probability 

conditioned on the subcarrier loop phase error (p in a power (Maclaurin) 

series in <f> and then, keeping only the first few terms of this series, 

average this conditional error probability over the probability density 

function (p.d.f.) of (j>. By doing this, we obtain the additional error 

probability due to a noisy subcarrier reference as an additive term 

2 

directly proportional to the mean-squared phase jitter directly 
associated with the receiver's subcarrier tracking loop. In this 
regard, the results derived in Appendix F play an important role in 
assessing this error probability performance. Finally, similar argu- 
ments can be advanced to give closed-form results for the noisy refer- 
ence loss itself. 

Appendix I contains a detailed discussion of the general problem 
of assessing the error probability performance of unbalanced QPSK 
receivers under the above-mentioned assumptions. The key results from 
this discussion are summarized as follows. In terms of the total 
signal-to-noise ratio R^2 ~ higher data rate band- 

width and the transmitted modulation indices ri2 ahd defined by 

(P„+P,)T ■ P,. 

'*T2 ■ ’ '’l ' P2 + P3 ’ 

the error probability performance of Channels 2 and 3 become 


= i- erfc 
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( 39 ) 


where ratio of data rates and the normalized mean- 

squared crosstalk iri2^2 tabulated in Table 2 of Appendix F for various 
combinations of data formats in the two channels. The corresponding 
expressions for the noisy reference loss (in dB) itself are 
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(40) 


In applying the results of (39) and (40) to Channels 2 and 3. to 
the three-channel SSO Ku-band return link signal, we again assume that 
m 2 (t) is NRZ data with a maximum rate Mbps and m 2 (t) is a Manchester 

coded data stream at 192 kbps. We further assume that the power 
allocation is chosen so that, for the given data rate ratio Y-j- = ^ 2/^3 
= 10.42, the signal-to-noise ratios R 2’^3 1''^° channels are made 

equal, i.e., both channels operate at the same error rate. Thus, 

R 2 = R 3 implies ^ 2^2 ~ ^ 3^3 modulation indices n 2 »n 3 become 


P + P 
-2 ^3 


Yt + 1 


= 0.9124 




P + P 
^2 ^3 


Yj+ 1 


= 0.0876 


(41) 








Further, it is typical to design the Costas loop bandwidth on the order 
of i? 2/100 (or less) since most of the power is in the high rate channel 
which controls the performance of the tracking loop. Thus, assuming 
Bi/i? 2 = 0*01 ,* Figures 32 and 33 illustrate L 2 and of (40) versus 
Bi/i ?2 for error probabilities of 10 "^, 10 '^, and 10 "^, corresponding 
respectively to R 2 = R 3 = 8.4, 9.6, and 10.5 dB. Several conclusions 
may be drawn from these figures. First, the noisy reference loss on 
Channel 2 is considerably smaller than that of Channel 3. The principal 


Smaller values of B[^/i ?2 as would be typical in practical receiver 
design would yield insignificant losses in L 2 and L^- 
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Figure 33. Channel 3 Noisy Reference Loss (in dB) vs. the Ratio of Two-Sided 
Costas Loop Arm Filter Bandwidth to Channel 2 Data Rate; m 3 (t) is 
192 kbps Manchester coded data, iTi 2 (t) is 2 Mbps NRZ data. 




reason for this can be easily explained in terms of the result in (40), 

where i t is observed that the effective cross-modulation loss on Chan- 
~ 2 

nel 2, m ^2 » is divided by which in this case has a value equal to 
10.42. Secondly, for either channel , the noisy reference loss decreases 
with increasing error probability. This is intuitively satisfying when 
one realizes that the slope of the error probability versus signal -to- 
noise ratio curve becomes steeper as becomes smaller and thus, for 
a given a^, the parallel ideal and noisy synchronization error proba- 
bility curves become closer together. Finally, we observe that there 

exists an optimum arm filter bandwidth (for fixed R^) in the sense of 

2 ^ 

minimizing L., i = 2,3. Since only depends upon this bandwidth, it 

is clear that this bandwidth choice is identical to that which minimizes 
2 

or, equivalently, the loop squaring loss. Note that, if is 

decreased, then the noisy reference loss will also decrease, since the 
equivalent loop signal-to-noise ratio p increases. 

In summary then, it is concluded that the crosstalk degradation 
due to noisy subcarrier demodulation references is quite small (on the 
order of tenths of a dB or less, depending on the particular channel 
and the ratio of loop bandwidth to data rate in that channel). When 
the higher data rate channel is 1 Mbps Manchester coded data, then since 
both channels are now Manchester coded, the crosstalk loss would be even 
smaller. 

While the results of Appendix I have been directed principally 
toward the demodulation of unbalanced QPSK by a conventional (single- 
channel) Costas loop, the expressions for average error probability 
[see (39)] and noisy reference loss [see (40)] apply in a much broader 
sense. In particular, the two-channel type Costas loops discussed in 
[25,26] have a mean-squared phase jitter given by (37) where, however, 
is a much more complex function of the various system parameters 

such as data rates and channel power ratios. Nevertheless, once s, 

2 L 

(and thus ) is determined, (39) and (40) apply directly toward 

evaluation of the noisy synchronization reference effects of these 

loops on error probability performance. Further discussion of the 

performance of two-channel Costas loops along these and other lines 

is presented in the next section. Finally, other possible applications 

of the results derived in Appendix I pertain to demodulation of 
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unbalanced QPSK using a biphase Costas loop with switching type multi- 
pliers. Once again, (39) and (40) apply, provided that an expression 
for the squaring loss can be found (see, for example, Appendix G). 



As part of the ongoing effort to investigate tracking techniques 
for improving communications efficiency on board the SS0» Axiomatix 
undertook to reexamine previous work [25,26] on the problems of carrier 
synchronization of an unbalanced QPSK (UQPSK) signal format, starting 
with the well-known maximum a posteriori (MAP) estimation technique as 
motivation for deriving closed loop tracking configurations. Such 
carrier reconstruction implementations are suggested upon examining 
the gradient of the likelihood function whose solution is the MAP esti- 
mator of carrier phase. 

In general, the closed loop implementation which results from 
application of MAP estimation theory to the problem of carrier phase 
estimation (herein referred to as the MAP estimation loop ) is impractical, 
primarily because of the difficulty of implementing the hyperbolic 
tangent nonlinearity suggested by using the gradient of the likelihood 
function as an error control signal in the loop. To arrive at practical 
realizations, one ordinarily approximates this nonlinearity with simpler, 
more easily implementable functions. In the past, authors have employed 
the approximations 



tanh X = sgn x *, 

X large 

(42) 

and 

tanh X = X ; 

X Small . 

(43) 


Since the input to the nonlinearity is a monotonic function of signal- 
to-noise ratio (SNR), then the approximations of (42) and (43) corre- 
spond, respectively, to conditions of high and low SNR. 

In Appendix J, we reexamine some of these approximations and the 
optimum structures which result. Using a slightly more sophisticated 
approximation of the nonlinearity, new and interesting structures are 
presented, first for QPSK as an example and then for UQPSK, which over- 
come some of the deficiencies posed by MAP receiver structures for these 



same modulations suggested by authors in the past. The equivalence of 
the new QPSK structure {Figure 34) with the well-known quadriphase Costas 
loop is discussed in Appendix K. Furthermore, it is demonstrated that 
the new UQPSK structure (Figure 35) allows carrier reconstruction at 
all ratios of data rates and powers in the two channels, even in the 
limit of balanced QPSK. 

In addition to arriving at new structures for optimum tracking 
of balanced and unbalanced QPSK, the practicality of the two-channel 
type carrier reconstruction loop for UQPSK is discussed in detail in 
Appendix J, paying particular attention to its sensitivity to variations 
in channel gains. Rather than go into the details of the performance 
of the loop in Figure 35, we instead investigated in Appendix J the 
sensitivity to gain variations of the simpler loop (Figure 36) found 
in [25,26] and used these results qualitatively as being indicative of 
two-channel Costas-type configurations. Actually, to make life even 
simpler, we performed our sensitivity analysis on the equivalent loop 
to Figure 36, where active (integrate-and-dump) arm filters were used 
in place of the passive arm filters and, furthermore, the channel gains 
2 /pJ/Nq and 2 /P^/Nq were replaced by arbitrary channel gains /K^ and /i^. 

It is shown there that, whereas the channel gains /K^ and /i^ 
were "optimally" chosen in [25,26] from MAP estimation considerations, 
a different selection of these gains based upon directly optimizing the 
loop's tracking performance can yield as much as a 10 dB improvement in 
this performance. The significance of this statement is not so much the 
fact that the tracking performance can be enhanced by a better choice of 
gains, particularly since these gains are now theoretically signal-to- 
noise ratio dependent, but rather the high degree of sensitivity of the 
performance of the MAP estimation-type loop to variations in the channel 
gains themselves. 

To provide quantitative verification of these points, Figures 37 
through 39 illustrate the squaring loss corresponding to the two 
possible selections of channel gains, versus the power ratio Xp with 
the data rate ratio A.|. = n (n integer) fixed and the total power- to-noise 
ratio in the high data rate bandwidth c-j - P^T-j/Ng = -3 dB and 10 dB. The 
values of n selected for these plots range from balanceo (n = 1 ) to 
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Figure 39: Squaring Loss Versus Power Ratio for Unbalcinced QPSK; Two-Channel 

Costas Loop; n = Tp/T-, = 10 
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highly unbalanced (n=10) data rates. Several points are worthy of note 
from the results in these figures. First, from Figure 37, we observe 
that, for the balanced data rate case (n=l, arbitrary), both^ chan- 
nel gain selections produce the identical squaring loss. Another case 
where the two channel gain selections produce the same result is that 
for which the channels have equal energy, i.e., A^^l/n. Since, from 
[25,26], the "optimum" value of K 2 /K 1 as motivated by MAP estimation 
theory i s 



(44) 


which is independent of signal-to-noise ratio, then for the equal energy 

2 

case, we have ~ ^p ‘ 

Second, the difference between the true optimum squaring loss 
and the value given by selecting the channel gains as in [25,26] 
increases with increasing n . Also, for fixed n, these differences 
are much more significant at higher values of total power-to-noise 
ratio in the high data rate bandwidth than at the lower values. As 
an example, for n= 10 and = 10 dB, we see from Figure 39 that, for 
equal powers (Ap = l), the optimum squaring loss is 10 dB smaller than 
that predicted by MAP estimation channel gain selection. The same com- 
parison at = -3 dB only shows a 1 ,4 dB improvement. 

Third, the true optimum (minimum) squaring loss decreases with 
increasing Cp regardless of the values of n and X, except for the 

case of small X. and n=l. 

P 

Finally, it can be shown that Figure 37 is representative of the 
one-channel Costas loop performance for all n and Apj hence, these 
curves represent lower bounds on the squaring loss performance for the 
two-channel configuration. That is, for fixed and Ap, the two- 
channel Costas loop will exhibit a smaller squaring loss for any data 
rate ratio n than the value indicated in Figure 37. This performance 
comparison between one- and two-channel Costas loops is made in Figure 40 

where the ratio 5^ = ^L|2.channev'®L|l -channel P’°“ed 

versus Ap with c-j (in dB) as a parameter and Ap = l/n. Clearly, for 

A approaching unity, the improvement becomes infinite; however, we 
P 


9 






recall that, in this limiting case, neither the one-channel nor the 
two-channel loop is capable of tracking at all, i.e., goes to zero 
for both. 

In conclusion, we point out that the MAP estimation theory pro- 
vides good intuition for implementing closed loop tracking configurations 
for UQPSK modulation; however, care must be exercised in selecting the 
ratio of gains between the two channels of the resulting configuration. 

In this regard, two options are available, namely, (1) that gain ratio 
which is motivated by MAP estimation theory and is independent of 
signal-to-noise ratio and (2) that ratio which minimizes the loop 
squaring loss and is dependent on signal-to-noise ratio. The differ- 
ence in squaring loss between these two choices of gain ratios can be 
as much as 10 dB for certain values of data rate ratio and power ratio. 
When the signal energies in the channels are chosen equal, then both 
gain ratios also become equal and yield identical squaring loss perform- 
ance. Even in this case, the two-channel Costas loop of [25,26] is 
capable of tracking balanced quadriphase. By an extension of the power 
series approximation used for the hyperbolic tangent nonlinearity which 
arises from the MAP estimation approach, we have been able to demonstrate 
a carrier reconstruction loop for UQPSK (Figure 35) which should yield 
better performance than the above-mentioned two-channel loop as the 
modulation becomes more balanced; in particular, it acts like a quadri- 
phase Costas loop in the limiting case of balanced QPSK. A detailed 
analysis of the tracking performance of this loop will be considered 
in the future. 

3 . 5 Costas Loop Tracking Performance on the Shuttle S-Band 

Uplink in the Presence of Residual Carrier 

We conclude the tracking techniques analysis task with an inves- 
tigation of the effects of residual carrier on the Costas loop tracking 
performance of the Shuttle S-band uplink. 

Traditionally, a Costas loop is intended for use in receivers 
which must reconstruct a carrier reference from an input signal v/hose 
carrier component is totally suppressed, e.g., a biphase modulated 
carrier. On the Shuttle S-band uplink, however, during the time when 
both data and a ranging subcarrier are linearly modulated on the same 


81 


€ 




m 


m 


m 












carrier, the data modulation index is not tt/ 2. Thus,- in this instance, 
the Costas loop is called upon to accurately track a signal whose 
carrier component is not completely suppressed. 

A simple block diagram which, for the sake of analysis, charac- 
terizes this situation is illustrated in Figure 41. Included in this 
illustration are the means by which the data and the ranging subcarrier 
are extracted using the in-phase demodulation reference generated by . 
the loop. 

Several key questions arise relative to the performance of the 
Costas loop under these unorthodox conditions: 

1. Is the loop capable of successfully tracking the input 
independently of the value of the data modulation index? 

2. Is it possible to extract the ranging subcarrier as shown 
if the data modulation is removed? 

3. What is the additional threshold power-to-noise ratio 
required to operate the loop in the PM (ranging and data) mode as 
compared to the PSK (data only) mode? 

4. What tradeoffs exist between power in the ranging channel 
and loop threshold performance as a function of the data and ranging 
modulation indices? 

The answers to these and other questions related to the perform- 
ance of Costas loops in the presence of residual carrier are given in 
Appendix L, where a theory is established for such performance as a 
function of the parameters which characterize the system tracking 
behavior. The key results from Appendix L are summarized below. 

First, it is demonstrated that, in the residual carrier mode of 
operation, there exists a critical data modulation index below which 
the loop will not operate, regardless of the value of signal -to-noise 
ratio. An expression for this critical modulation index is given as 
follows. Consider the Costas demodulator illustrated in Figure 41 
whose input signal is of the form 

s(t,0) = sin [togt + Bj d(t) + sin (co^^t + 0^(t)j + e] , (45) 

where P is the total received power, Wq is the carrier radian frequency 
and 0 the corresponding input phase to be estimated, is the data 









modulation index with d(t) the data waveform, and is the ranging 
subcarrier modulation index with u the radian subcarrier frequency 

V w 

and 0^(t) the tone ranging modulation. Then, if 



<d^(t)> = 


S (f)lG(j2;,f)l^df 

•'—OO 


( 46 ) 


represents the power in the data modulation at the arm filter output 
with Sj^(f) denoting the power spectral density of the unfiltered data 
d(t), when 


cot 



(47) 


the loop S-curve vanishes and the loop will not lock at any loop signal- 
to-noise ratio. For values of less than the critical value satis- 
fying (47), the loop locks at ^ = ±(2n+l)Tr/2, n=0,l,2, For 

values of 3^ greater than this critical value, the loop locks at 
cp = ±mr, n=0,l,2,.... 

Table 5 tabulates the critical values of data modulation index 
for the various data symbol rates of interest on the S-band uplink, 
wherein the data is Manchester coded and the Costas loop arm filters 
are assumed to be single-pole Butterworth (RC) with 3 dB cutoff fre- 
quency f and two-sided noise bandwidth B. = Trf . 

, In* 


Table 5. Critical 

Values of Data 

Modulation 

Index on the S- 

-Band Uplink 

R 

s 



m 

^^d'crit 

32 kbps (Low 

Data Rate Uncoded) 

134 kHz 

13.155 

0.886 

0.8157 rad 

96 ksps (Lov^ 

Data Rate Coded) 

134 kHz 

4.385 

0.6636 

0.8872 rad 

72 kbps (High 

Data Rate Uncoded) 

308 kHz 

13.44 

0.8884 

0.815 rad 

216 ksps (High 
Data Rate Coded) 

308 kHz 

4.48 

0.670 

0.8848 rad 
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Second, the performance of the ranging channel is characterized 
both when data modulation is on and when it is off. For either case, 
the signal power in the ranging channel is given by 

= 2P cos^ (Sg) cos^(j£? (48) 

where, for large 

cos^ (p = 1 — ^ . (49) 

Then, when the data modulation is on, we note from (45) that, as the 

value of is increased beyond its critical value as determined from 

(47), the carrier becomes more suppressed and thus the Costas loop 

tends to track better. However, as 3^ is increased, cos 6^ decreases 

and thus [from (48)] the power in the ranging channel P^ tends to 

decrease because of this. Since the power in the ranging channel is 

2 

proportional to both cos and the accuracy of Costas loop tracking 
(through the factor cos^ cp), a tradeoff exists with regard to the 
selection of the data modulation index- 

When the data modulation is removed, i.e., then as pre- 

viously mentioned, the loop will now lock up around cp = ti/Z (as opposed 
to (p 0) . This is equivalent to saying that the in-phase and quadra- 
ture demodulation reference signals switch roles. Thus, since the 
ranging subcarrier component in the input signal s(t,e) is now demodu- 
lated by the quadrature carrier reference signal, then from (48), we 
have that (for large loop signal -to-noise ratio) P^=0. The conclusion 
to be reached then is that, in order to extract range information as 
in Figure 41, the data modulation cannot be turned off at the transmitter. 

The third key result concerns the comparison of threshold tracking 
performance in the Pfi versus the PSK mode. Tracking threshold is typi- 
cally specified as the minimum input signal level for which the mean 
time to loss of lock is greater than or equal to 10 seconds. This con- 
dition is uniquely specified by the loop signal-to-noise ratio or, 
equivalently, • Thus, two modes of operation, e.g., PM and PSK, 
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will have identical tracking thresholds if the total powers to-noise 
ratio P/Nq is adjusted in one relative to the other to produce the same 
1n both. 

The baseline design of the S-band network transponder under 
development by TRW calls for rad and 3^ = 1.1 rad. Performance 

tests on the engineering model, as reported in the network transponder 
CDR package [27], reveal the following tracking threshold signal levels 
in the PSK mode (through TORS): 


Table 6. Tracking Threshold Signal Levels in the PSK Mode 


Data Rate 

Temperature 

Signal Level 

Low Data Rate 

70°F 

-104.2 dBm 

(32 kbps uncoded or 

120°F 

-102.4 dBm 


t 

96 ksps coded) 

-20®F 

-104.2 dBm 


High Data Rate 

70°F 

-103.7 dBm 


(72 kbps uncoded or 

120°F . 

-101.4 dBm 

• 

216 ksps coded) 

-20°F 

(SSP will not hold 



lock -100 dBm) 


The input noise power spectral density is -151 dBm/Hz. Thus, 
the maximum measured P/N^ is 48.6 for the low rate and 49.6 for the 
high rate modes. Assuming these threshold values of P/Nq for PSK 
operation through the TDRS, Table 7 tabulates the values of various 
system parameters leading up to the calculation of P'/Ng (total power- 
to-noise ratio in the PM mode). Also indicated in the table is the 
dB increase in power- to-noise ratio required to go from the PSK to the 
PM mode through the TORS. 

Table 8 illustrates the power- to-noise ratio P^/Ng in the ranging 
channel [as computed from (48)] for the same parameters as in Table 7 
and operation in the PM mode. 

Finally, tracking threshold performance is examined as a func- 
tion of data and ranging modulation indices. This threshold, when 
operating in the PM mode, can be reduced (relative to that in the PSK 
mode) by reducing the ranging modulation index and/or increasing 
the data modulation index (i.e., suppressing the carrier more). 


Table 7. System Parameters for Calculation of P'/Nq 
P/Nq = 48.6 dB, = 1.0 rad, 3^, = 1-1 rad 


Data Rate 

D m 

PT/Nq (dB) 

Pi (dB) 

P7Nq (dB) 

aP/Nq (dB) 

32 kbps 

0.829 

3.549 

-4.642 

54.23 

5.63 

96 ksps 

0,5078 

-1.223 

-4.642 

55.01 

6.41 

72 kbps 

0.8326 

0.0267 

-8.257 

53.83 

5.23 

216 ksps 

0.5166 

-4.745 

-8.257 

54.50 

5.90 


Table 8. Power- to-Noise Ratio in the Ranging Channel 

P/Nq - 48.6 dB, = 1.0 rad, 3^ = 1-1 rad, 

= 200 Hz 


Data 

Rate 


(dB) 

32 

kbps 

43. 

21 

96 

ksps 

43. 

.97 

72 

kbps 

42. 

,78 

216 

ksps 

43. 

,41 
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Unfortunately, however, both of these changes also reduce the power- 
to-noise ratio in the ranging channel. To see these effects quanti- 
tatively, Figure 42 illustrates the tracking threshold power-to-noise 
ratio for PM operation as a function of ranging modulation index for 
fixed data modulation index, and Figure 43 illustrates the ranging 
channel power-to-noise ratio as a function of these same parameters. 
The value of P/N^ for PSK operation is again chosen equal to 48.6 dB 
and the results are given for both uncoded and coded low data rates. 

We observe from Figure 42 that, as 6^ approaches it/ 2 (fully suppressed 
carrier), the tracking threshold becomes virtually insensitive to data 
rate. A similar phenomenon is observed in Figure 43 for the ranging 
channel power-to-noise ratio. Using the information in these curves, 
along with the ranging channel bandwidth and the required signal-to- 
noise ratio in this bandwidth turned around to ground, one can select 
values of 6^ and to meet the additional requirements on increase 
in tracking threshold for PM versus PSK operation. 


Figure 42. Tracking Threshold Power-to-Noise Ratio for PM Mode vs. Ranging 
Modulation Index with Data Modulation Index as a Parameter 
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4.0 PAYLOAD COMMUNICATIONS STUDY 

The Ku-band system currently provides a return link wideband 
and narrowband bent-pipe. This study assessed the design and perform- 
ance of the return link bent-pipe relay mode for attached and detached 
payloads. As part of this study, a design for a forward link bent-pipe 
is also proposed, and its performance is evaluated for attached and 
detached payloads. Critical Ku-band system parameters have been 
identified and optimized such that cost and complexity have been 
minimized, and the basic data transmission capabilities that can be 
accommodated for payloads have been identified for both forward link 
and return link bent-pipe relay modes. 


Forward Link Bent-Pipe Relay Design 


The present Ku-band system design does not include provision for 
a forward link bent-pipe relay mode. This section proposes a possible 
modification to the Hughes Aircraft Company (HAC) Ku-band forward link 
design to accommodate forward link bent-pipe signals to attached or 
detached payloads. A simplified block diagram of the HAC forward link 
design is shown in Figure 44. The ''F receiver is in the Deployed 
Electronic Assembly (DEA). The input to the Electronic Assembly 1 
(EA-1) is at 647 MHz IF. The EA-1 provides IF ampl ification, automatic 
gain control (AGC), and automatic level control (ALC) for the remaining 
demodulation circuits in EA-1. The PN spread spectrum modulation is 
removed by the PN tracking loop while the Costas loop tracks the 
carrier to complete the demodulation of the data. The output of the 
EA-1 is buffered with an AGC. The input to the Signal Processor 
Assembly (SPA) is either switched to the Network Signal Processor 
(NSP) directly or bit and frame synchronized in the SPA. In the 
second mode, the demultiplexed 72 kbps operational data is sent to 
the NSP while the demultiplexed 128 kbps data is sent to the payload 
processor. In order to implement a forward link bent-pipe mode, an 
additional output is needed from the Costas loop and an additional 
switch is needed in the SPA. The switch in the SPA allows either 
normal data, bent-pipe data modulated on the carrier, or bent-pipe 
data modulated on a subcarrier to be output to the payload. 




















The outputs of the Costas loop for normal data and the bent-pipe 
data are presented in Figure 45. In the normal data demodulation, the 
Costas loop has been designed to reconstruct a carrier reference from 
the input signal whose signal component is totally suppressed (e.g., i 
a biphase modulated carrier). For the bent-pipe mode, the Costas loop 
must accurately track a signal whose carrier component is not comp’^jtely 
suppressed. As is the case for a Costas loop with a ranging waveform 
modulated on the subcarrier of the received signal, several key questions 
arise relative to the performance of the Costas loop when bent-pipe 
data is present. These questions are repeated here for convenience. 

1. Is the loop capable of successfuly tracking the input 
independent of the value of the data modulation index? 

2. Is it possible to extract the bent-pipe data modulated 
on a subcarrier if the normal data modulation is removed? 

3. What is the additional threshold signal-to-noise ratio 
required to operate the loop in the PM (i.e., residual carrier) mode 
as compared to the PSK mode? 

4. What tradeoffs exist between the power in the bent-pipe 
channel and loop threshold performance as a function of the normal 
data and the bent-pipe data modulation indices? 

Appendix L presents the detailed analysis to answer these questions for 
the general case of Costas loop tracking in the presence of a residual 
carrier. To determine the performance of the Costas loop with the 
bent-pipe data, consider the input signal of the form 

s(t,0) = v^P sin [o)Qt+ d(t) + 3 ^ sin (wg^t+ 0 j^p(t)) + e] , (50) 

where P is the total received power, Wq is the carrier radian frequency 
and 0 is the corresponding input phase to be estimated, 3^^ is the data 
modulation index with d(t) either the normal data waveform or the bent- 
pipe data v/aveform when modulated on the carrier, and 3^ is the bent- 
pipe subcarrier modulation index with u) the radian subcarrier fre- 
quency and the bent-pipe data waveform. Using the results 

presented in Section 3.5 and Appendix L, the performance of the sub- 
carrier bent-pipe channel is characterized both when 3^j is equal to zero 
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(i.e., data modulating the carrier directly is not present) and when 3^ 
is nonzero (i.e., data modulating the carrier directly is present) by 


the bent-pipe signal power Pj^p to the total power P. 


= 2 cos^ 


1 - 


pS. 


(51) 


where p is the signal-to-noise ratio in the loop bandwidth given by 

P 


p = 




(52) 


and where is the squaring loss in the Costas loop. The squaring 
loss is derived in Appendix L to be 




0o(6s)[cos (6<,) + D^KpSin ^ (bJ — 


(53) 


for NRZ data. The data distortion factor D for NRZ data with RC 

m 


filters as the lowpass arm filters in the Costas loop is given by 

1 


D„ = 1 


m 




1 - exp 


2B.^ 


R 


(54) 


The factor in (53) is found in Appendix L to be 

3 - (3+2B./i?g) exp (-2B.//?g) 


m D 


= 1 - 


4B|//?, 


(55) 


The term is the ratio of the two-sided noise bandwidth of the 


lowpass arm filters in the Costas loop the the data rate. The signal 
energy per bit-to-single-sided noise spectral density, R^, is given by 

n _ P 




(56) 
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Figure 46 presents the ratio of subcarrier bent-pipe signal power 
Pbp to total power P as a function of modulation indices. The curves 
of Figure 46 are computed from (51). It may be observed that the curves 
for = 5 dB and R^= 10 dB coincide, and hence the bent-pipe data 
demodulation performance is insensitive to R^. 

Another key performance parameter is the additional power required 
to operate the Costas loop in the PM mode. Let P' denote the signal 
power in the PM mode to achieve the same performance as PSK data with 
power P. Then the ratio of P'/P is the factor by which the signal power 
must be increased in order not to degrade the normal data bit error 
probability. Appendix L shows that 



where 

A = ^ 0 ^( 35 ) sin^ (e^) - cos^ (e^)] 


B = 


Jo^(35)[cos^ (^d^^^mS ^^d^^ 


m 


Dm Kn + 
m D 




4R 


d 



(W 

'"'d 




4R^ 





(58) 


Figures 47 and 48 present the ratio P'/P as a function of modu- 
lation indices for R^j = 5 dB and Rjj = T0 dB, respecti vely . From these 
figures, it may be sesn that, as increases (i.e,, approaches PSK), 
the ratio P'/P approc les 1. Also, is increased, there is only 

a small incree.f:?. in e ratio P'/P for large 3 ^ but, for 3 j = l.l> 
increasing from 5 dB to 10 dB increases the required P'/P by 0.4 dB 
The tracking performance of the Costas loop in the PM mode is 
also an important design parameter. The tracking phase jitter which 
degrades the power in the' data and subcarrier bent-pipe channels has 


variance 








where p and 5, are given by (52) and (53), respectively. For B =1.0 
and 3^ = U , the standard deviation of the tracking phase jitter 
in percent radians is presented in Figure 49 as a function of and 
the ratio of Costas loop arm filter to normal data rate. Note that 
increasing the ratio of greater than 2 drastically degrades the 

tracking performance. The optimum value of for the PSK mode is 

1.4 and this seems to be a good value for the PM mode as well. There- 
fore, the Costas loop need not be modified for bent-pipe operation 
but the additional output from the quadrature arm is required in EA-1 
if bent-pipe data is to be modulated on a subcarrier. Also, as Figures 
46 through 48 have shown, there will be some degradation to the normal 
data in the PM mode, depending on the modulation indices. 

4.2 Return Link Bent-Pipe Relay Design 

The return link bent-pipe mode for payloads allows transmission 
of data that is not in the standard NASA format. Thus, multiple formats 
and multiple modulators/demodulators can be used by the payloads. The 
design goal for the bent-pipe mode is to minimize the signal processing 
in the Orbiter for data that does not meet the standard NASA format. 
Therefore, no Orbiter control of command or telemetry for the bent-pipe 
data is proposed. Rather, the Orbiter will act as a relay which will 
either merely make a frequency translation at IF or perform RF demodu- 
lation and remodulation on a new carrier. ■ 

The block diagram for the return link bent-pipe is shown in 
Figure 50. Channel 1 is primarily an analog commercial TV channel 
with a bandwidth of 4.2 MHz. /nternately, however, Channel 1 may 
consist of one of the other data channels shown, namely, 4.5 MHz of 
analog data or up to 4 MHz of NRZ digital data. 

The remaining channels are to be modulated onto an 8.5 MHz square 
wave unbalanced QPSK subcarrier. The power split is 4:1. The 20% input 
is designated as Channel 3 and consists of a 192 kbps Bi-(f)-L data 
channel from the network signal processor. The 80% input, Channel 2, 
is generally from the payload system and consists of up to 2 Mbps of 
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Figure 50. Return Link Bent-Pipe Link Implementation 












data in the NRZ format. As shown in Figure 50, Channel 2 may also con- 
sist of information from recorder playbacks. 

Among the various candidate signal formats that are expected to 
be transmitted via nonstandard bent-pipe is a 16 kbps Bi-(ii-L data 
stream on a 1.024 MHz subcarrier. A block diagram of the transmitter 
is shown in Figure 51. Figure 52 shows a block diagram for a possible 
implementation of the ground receiver. At the receiver, the post- 
detection filter for the 8.5 MHz subcarrier has a bandwidth of 7.6 MHz. 
The LPF for the analog TV channel is 4.2 MHz. The 16 kbps data signal 
on the 1.024 MHz subcarrier was also tested in Channel 1 as shown in. 
Figure 52. The bandwidth of the prefilter for the 1.024 MHz subcarrier 
is set at 160 kHz. 

The FM deviations are set at the transmitter by adjusting the 
gains G 2 and G^. each individually with the other signal inputs 
removed, as follows: 


Gain 

Channel 

+ AF 

Si 

8.5 MHz subcarrier 

6 MHz 

S2 

TV - Channel 1 

11 MHz 


16 kbps on 1.024 MHz sub- 

5 MHz 

0 

carrier in Channel 1 



The total range of the instantaneous frequency deviation is 2 aF in each 
of the above cases. 

In the wideband bent-pipe mode, the bent-pipe payload data is 
transmitted in Channel 1. The concern in this mode is that, for low 
signal bandwidths (or low data rates), the S/N in the 4.5 MHz band- 
width at the transmitter will be so small that the FM discriminator 
on the ground will be below threshold. However, the S/N in the 4.5 MHz 
bandwidth at the transmitter will not affect the performance of the 
bent-pipe mode, but rather the S/N in the signal bandwidth is the 
parameter that determines the performance. This occurs because the 
noise at the transmitter is treated as signal by the FM discriminator 
so FM threshold will not occur unless the total received signal power 
(i.e., signal plus noise from the transmitter) is too small compared 
to the noise power at the receiver. Thus, the FM discriminator restores 
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Figure 51. Ku-Band Mode 2 Ff1 Transmitter for Return Link Bent-Pipe Mode 
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the signal -plus-noise waveform that was present at the transmitter. 
Following the FM discriminator, the post-detection signal demodulator 
must contain a filter matched to the transmitted signal if the signal 
is to be optimally demodulated. If the post-detection filter is 
matched to a signal with a 4.5 MHz bandwidth, then demodulating a 
16 kbps signal will not be optimum, and required S/N in the 4,5 MHz 
bandwidth will have to be increased drastically. Therefore, the per- 
formance of the wideband bent-pipe mode is a function of the ground 
post-detection demodulation filter being matched to the transmitted 
signal and not due to the FM threshold. 

In the narrowband bent-pipe mode, the bent-pipe payload data 
is transmitted in Channel 2. The main concern in this mode is the 
S/N in the 4,5 MHz bandwidth of the Payload Interrogator (PI) will be 
so small that the ground subcarrier tracking loop presented in 
Section 3.2 will lose lock and the operational data will also be 
lost. If the S/N in the PI 4.5 MHz bandwidth is small, then the 
effective data rate in Channel 2 will be 4.5 Mbps. Since the ground 
subcarrier tracking loop is designed for a maximum of 2 Mbps in 
Channel 2, the signal distortion due to lowpass filtering of 4.5 Mbps 
will increase, and the squaring loss of the Costas loop will increase 
also. Therefore, the ground subcarrier tracking loop performance will 
degrade and lose lock at a much higher received S/N in the 2 MHz band- 
width. In order to establish the required S/N in the PI 4.5 MHz band- 
width for the narrowband bent-pipe mode, a test of the HAC digital 
implementation of the subcarrier UQPSK modulator should be made for 
the narrowband bent-pipe signals at low S/N. Depending on the effective 
bandwidth of the HAC subcarrier modulator, the PI 4.5 MHz bandwidth 
may have little significance. 

A proposed way of combating low S/N in PI bandwidth is to devise 
a squelch circuit that would not send a signal to the Ku-band system 
from the PI unless the S/N in the PI bandwidth was greater than 7 dB. 

The design of the squelch circuit is a problem because, without knowl- 
edge of the signal structure, which is unlikely in the bent-pipe mode, 
it is difficult to determine the signal waveform from noise. It should 
be noted that, instead of a squelch, the ground operations could always 


command the Ku-band system out of the narrowband bent-pipe mode if the 
ground subcarrier tracking loop loses lock and the operational data 
is not being demodulated. 

To solve the low S/N problem in the narrowband bent-pipe mode, 
a test should be made with the HAC digital subcarrier UQPSK modulator, 
and an analysis of low S/N signals in Channel 2 should be developed to 
compare against the test results. 
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BIT ERROR RATE DEGRADATION OF THE KU-BAND RETURN LINK CHANNEL 
DUE TO NRZ DATA SYMBOL ASYMMETRY 

by 

Marvin K. Simon 


The high data rate link from the Shuttle Orbiter through the TDRSS 
to the ground takes NRZ symbols at 50 Mbps and encodes them with a rate 1/2, 
constraint length 7, convolutional code. The bit error rate performance 
of the convolutional decoder depends, among other things, on the symmetry 
of the modulation. Any asymmetry in the NRZ symbols entering the symbol 
synchronizer causes a misalignment in the synibol synchronization clock 
which degrades the integrate-and-dump output and any soft or hard decisions 
derived from it for input to the decoder. For a specified degree of asym- 
metry (in terms of a fraction of a symbol interval), the bit error rate 
degradation is dependent on the transition probability of the data.' 

Clearly, if the data transmitted was either all ones or all minus ones, 
then misalignment of the bit synchronization clock would have no degrading 
effect on the integrate-and-dump output since, for each symbol, this cir- 
cuit would integrate up to its maximum value before being dumped. On the 
other hand, when the data is an alternating sequence, then the worst case 
degradation results, since the transition which occurs at the end of each' 
symbol in combination with the symbol synchronization clock misalignment 
prevents the integrate-and-dump output from reaching its maximum value. 

To quantitatively determine the degrading effect of NRZ symbol 
asymmetry on error rate performance, consider first the alternating NRZ 
sequence illustrated in Figure la, where the +1 symbols are elongated by 
AT/2 (relative to their nominal value of T sec) and the -1 symbols are 
shortened by the same amount. Thus, AT represents the relative difference 
in length between the +1 and -1 symbols. The data asymmetry is defined 
as the ratio of the difference in length between the +1 and -1 symbols 
to the sum- of their lengths, i.e.. 


Asymmetry 



T(l+|) - T(1 -|) 
T(l+|) + T(l-|) 




In the absence of noise-, the timing instants for the in-phase 
integrate-and-dump (i .e. , the epoch of the symbol synchronization clock) 
are determined as follows. The mid-phase integrate-and-dump in the symbol 






m 








II 
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I « 


synchronizer integrates across the transitions in the data symbol stream 
and determines the magnitude of the symbol synchronization error signal. 

In the steady state, this error signal must have zero value, on the average 
Letting eT denote the misalignment of the symbol synchronization clock, 
then from Figure lb we see that, when integrating across a negative data 
transition, the output of the mid-phase integrate-and-dump is given by 


V_ = (]- + f -c)T - (i- |+e)T = (A- 2 e)T. ( 2 ) 


When integrating across a positive data transition, the same output is 
given by 

V+ = - (]• - e)T + +e)T = 2eT . . (3) 

Figure Ic illustrates the in-phase integrate-and-dump output for the two 
types of data transition (negative and positive). For the negative data 
transition, we notice that the in-phase integrate-and-dump output reaches 
its maximum value and thus no degradation results. For the positive data 
transition, this same output is degraded by the factor (1 -4e). Thus, 
since the symbol synchronization clock misalignment e is determined by 
equating the average of the product of the mid-phase integrate-and-dump 
output and the sign of the difference of two adjacent in-phase integrate- 
and-dump outputs (or, equivalently, the sign of the data transition) to 
zero, i.e., V^-V_ = 0. Thus, from (2) and (3), 


2 sT - (A - 2 e)T = 0 , 

( 4 ) 

A 

^ " 4 • 

( 5 ) 


Assuming then a clock misalignment as in (5) and letting E denote the 
symbol energy and Nq the channel noise spectral density, then the average 
symbol error probability associated with hard decisions made on the • 
in-phase integrate-and-dump outputs is given by 


= I erfc (v^f^N^) + I erfc [/E^ (1 - A)] , 


where 


X- A 2 

erfc X = 


/n^ 


dt , 


( 6 ) 

( 7 ) 


X 


and the “A" superscript on refers to the result for an alternating 
data symbol sequence. 

When the input data is random with equiprobable symbols, then the 
in-phase integrate-and-dump output depends, in general, on the polarity 
of the symbol over which it is integrating and that of the preceding and 
succeeding symbols, Thus, one must compute the in-phase integrate-and- 
dump output for each of the eight possible three-symbol sequences (see 
Figure 2) and their corresponding conditional error probabilities. Then 
averaging these conditional error probabilities over the equal probabili- 
ties of the eight equal ly-1 ikely three- symbol sequences gives the average 

D 

error probability P£‘, i.e.. 


J " 2[i erfc (.€ 75 ^ (1 -|)[ 


1 


- t erfc 


(1-4)) 


^ erfc i erfc 


+ ^ erfc 0 - 4 ) (8) 


where the "R" superscript on P^ refers to the result for a random data 
symbol sequence. 

Tables 1 and 2 contain the symbol energy- to-noise ratio degradations 
(in dB) for asymmetry values [see (1)] of 3, 7, 10, 15 and 20% (A = 0.06, 
0.14, 0.20, 0.30, and 0.40) and E^/Nq= 0, 0.75, and 1.5 dB. The values 
of E^/Nq selected correspond to bit energy- to-noise ratios Ej^/Nq = 3, 3.75, 
and 4.5 which, respectively, correspond to decoder bit error probabilities 
P|^=10''“’, lO”^, and 10~ . The degradations are obtained from (6) and (8) 
by computing the additional E /N« required due to asymmetry to produce 

s ^ 0 

the same value of symbol error probability when A=0, i.e,, P^ , where 

^ 1 erfc (/ET^) . (9) 

It should be noted that the symbol energy- to-noise ratio degradations given 
in these tables assume no channel bandwidth limitation, i.e., ideal rec- 
tangular pulse shapes have been assumed for the NRZ data. Any rounding 
of the pulses caused by channel bandwidth limitation produces an additional 
symbol energy- to-noise ratio degradation over and above that due to data 
asymmetry. However, the bandl iniiting degradation and the data asymmetry 
do not add algebraically but rather combine in a way determined by the 
particular symbol synchronization implementation. 
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Table 1. Alternating Data 


Asymmetry [%) 

E3/N(,(dB) 

Pb 

Degradation (dB) 

3 

0 

Eh 

0.275 




0.275 



H 

0.275 

7 

0 

10"3 

0.68 



10"^+ 

0.69 


T . 50 

10"5 

0.70 

10 

0 

Eh 

1.025 


0.75 

mm 

1.045 


1.50 


1.070 

15 

0 

10"3 

1.693 


0.75 

mm 

1.748 


1.50 


1.811 

20 

0 

10“3 

2.530 


0.75 

lO"*^ 

2.645 


1.50 

10"5 

2.778 


Table 2. Random 

Data 


Asymmetry (%) 

E^/Nq (dB) 

Pb 

Degradation (dB) 

3 

0 

10"3 

0.135 


0.75 

lO-'t ' 

0.135 


1.50 

10" 5 

0.135 

7 

0 

10"3 

0.333 



10“'+ 

0.337 


1.50 

10-5 

0.340 

10 

0 

10"3 

0.495 


0.75 

10-4 

0.505 


1,50 

10"5 

0.517 

15 

0 

10"3 

0.799 


0.75 

10-4 

0.824 


1.50 

10“ 5 

0.854 

20 

0 

10“3 

1.149 


0.75 

10“'+ 

1.201 


1.50 

10”5 

1.264 


# 
















For uncoded NRZ data, we find from (9) that, for P^ = 10~^, we 
require E^/Nq= 9.6 dB. The corresponding additional E^/Nq (SNR degrada- 
tion) due to asymmetry required to produce a 10”^ error probability is 
illustrated in Figure 3 as a function of percent asymmetry. In this, 
figure, we have let denote T(1 + A/2) and T_ denote T(1 - A/2) in accord- 
ance with the definition of asymmetry given in (1). Also illustrated in 
this figure is the SNR degradation due to asymmetry where the asymmetry 
is in accordance with Goddard's definition given by 


\ 

Asymmetry = — ? 


2T 


T(H-f) - T(1-|) 
2T(1 -j) 


2(1 - f) 


(10) 


or equivalently. 


A = 2 (Asymmetry) 

1 + Asymmetry 


( 11 ) 


Thus, for a given asymmetry, A as computed from (11) is substituted in (8) 
from which the necessary E^/Nq to yield P^^=10'^ is computed. Comparing 
this E^/Nq with 9.6 dB yields the SNR degradation plotted in Figure 3. 

In conclusion, we note that an analysis of the data symmetry prob- 
lem is also presented in [1,2] using, however, a slightly different model 
for the way in which the asymmetry comes about. The assumption made there 
is that positive NRZ pulses are shortened whenever adjacent pulses are 
negative. Thus, a given positive pulse preceded and succeeded by a nega- 
tive pulse would be reduced in duration at both ends. Letting 5 represent 
the fractional (relative to the nominal bit duration T) increase in posi- 
tive pulse length to a single adjacent negative pulse, then for a given 
random data sequence, the longest pulse would have Iqngth 1 + 26, while 
the shortest would have length 1-25. If, as before, asymmetry is 
defined as the difference in length between the shortest and longest 
pulses in the sequence divided by their sum, then using the model in 
[1 ,2], we get 


Asvmmetrv = ^(1 + 26) - T(1 - 25) ^ 
Asymmetry t(1 + 26) + T(1 - 26) • 


( 12 ) 


Also, using this same model, it can easily be shown that, on the average, 
the symbol synchronizer will lock up at the nominal transition points of 
the equivalent symmetric data waveform, i.e. , 0, T, 2T, 3T, . ... 


SNR Degradation (dB) 
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Finally, using the definition of (12) and the above statement 
regarding the relative timing of the symbol synchronization clock and the 
data, the result given in [1,2] for average error probability due to data 
asymmetry is identical to our Equation (8). Thus, properly applied, 
either symbol synchronization asymmetry model will yield the same results. 
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FURTHER CONSIDERATIONS ON THE EFFECTS OF DATA ASYMMETRY ON 
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S-BAND RETURN LINK CHANNELS 

by 

Marvin K. Simon 

In Appendix A and [1], we investigated the bit error rate degrada- 
tion of the Ku-band return link channel due to NRZ data symbol asymmetry. 
Inherent in these results was the assumption that the symbol synchronizer 
was of such a type (e.g., data transition tracking loop [2]) as to pro- 
duce a symbol sync clock which, on the average (e.g., no noise), locks 
up with a misalignment equal to half the asymmetry. While this assump- 
tion is believed to be typical of a wide variety of symbol synchronizers 
which employ an Integra te-and-dump filter as a data detector, one might 
nevertheless ask the question: How sensitive is this SNR degradation due 

to asymmetry to perturbations of the symbol sync clock misalignment about 
its above-assumed value? Having this information would also tend to make 
the results somewhat independent of symbol sync configuration. In this 
appendix, we derive an expression for symbol error probability in the pre- 
sence of data asymmetry and conditioned on an arbitrary clock misalignment 
whose value is allov/ed to range over the data asymmetry interval. SNR 
degradation for the coded NRZ case is then given as a function of this 
clock misalignment with percent data asymmetry as a parameter. 

Despite the fact that the high rate Ku-band return link is coded 
NRZ data, it is also desirable to have theoretical results for SNR degra- 
dation due to asymmetry corresponding to coded Manchester data. This 
would allow comparison with measurements previously made by TRW on the 
S-band return link. Thus, the second part of this report has as its end 

result a plot of SNR degradation in dB versus percent asymmetry for coded 

-5 

Manchester data at a bit error rate of 10 (symbol SNR equal to 1.5 dB). 
Also plotted on the curve for purposes of comparison will be the corre- 
sponding results (taken from Table 2 of [1]) for NRZ data. 

Probability of NRZ Symbol Error in the Presence of Data 
Asymmetry and Arbitrary Clock Misalignment 

For random NRZ data with equi probable symbols, the in-phase 
integrate-and-dump output depends, in general, on the polarity of the 
symbol over which it is integrating and that of the preceding and 


2 


succeeding symbols. Letting eT denote (as in [1]) the misalignment of 
the symbol sync clock, then the in-phase integrate-and-dump output for 
each of the eight possible three-symbol sequences is tabulated below. 

In-Phase 

Ihtegrate-and-Dump 
Symbol Sequence Output 


1 1 1 

1 1 -1 

1 -1 1 

1 -1 -1 

-11 1 

-1 1 -1 

- 1-1 1 

-1 -1 -1 


T 

T 

-T(l-A) 

-T(l - A + 2e) 
T 
T 

-T(l-2e) 

-T 


Thus, since each of these three-symbol sequences are equally likely, the 
average probability of error conditioned on e is given by 

Pe'^( 0 = ^ erfc + Jg erfc [/E^ 

+ ^erfc [vf^(l- A+2e)] +^ erfc (1-A)3 

( 1 ) 

where the "R" superscript on refers to the result for a random data 
sequence. First note that, for e = A/4, (1) reduces to (8) of [1]. 

Secondly, P^ (e) is a symmetric function of e around the point A/4, i.e., 
equal positive and negative variations of e around the nominal position 
A/4 produce equal SNR degradations above that previously given in [1]. 
Figure 1 illustrates the SNR degradation at a symbol E /N« of 1.5 dB 

S 

(corresponding to a decoded bit error probability of 10"^) as a function 
of £ for 0 X e ^ A/2 and percent asymmetry A/2x 100 as a parameter. We 
observe that, in the neighborhood of the nominal symbol sync lockup point, 
e = A/4, the sensitivity of SNR degradation due to asymmetry is extremely 
small, even for large asymmetry values on the order of 25%. Thus, we 
conclude that the results given in [1] are relatively insensitive to the 
assumption of a hard-locked symbol sync at the midasymmetry point, even 
though this assumption yields the minimum SNR degradation due to asymmetry. 


SNR Degradation (dB) 





Probability of Manchester Symbol Error in the 
Presence of Data Asymmetry 

When Manchester coding is employed, then relative to the NRZ 
sequence, the Manchester coded waveform has 3/2 as many transitions. 

Thus, since SNR degradation due to asymmetry is directly related to the 
average transition density of the data sequence, one would intuitively 
expect that, for a fixed amount of asyironetry (in seconds) , the Manchester 
coded case should yield a larger SNR degradation than the corresponding 
NRZ case. Indeed, this is true, as we shall soon demonstrate. On the 
other harid, depending on how one defines percent asymmetry for the Man- 
chester coded case (two possibilities are suggested), a plot of SNR 
degradation versus percent asymmetry might result in a different con- 
clusion than the above. This shall also be demonstrated by numerical 
example. 

Consider an NRZ sequence and the corresponding Manchester waveform 
with asymmetry illustrated in Figure 2. Here A/2 denotes the fractional 
(relative to the half symbol time T/2) elongation of the positive half 
pulse in the Manchester coded waveform. Once again, as in the NRZ case, 
the in-phase integrate-and-dump output depends, in general, on the polarity 
of the symbol over which it is integrating and the preceding and succeeding 
symbols. For the eight possible three-symbol sequences, the in-phase 
integrate-and-dump output is tabulated below, assuming a nominal bit 
sync lock-up misalignment of aT/8. 

In-Phase 

In teg rate -and -Dump 

Symbol Sequence Output 

111 ■ T(1 - A/2) 

1 1 -1 T(1 - A/4) 

-1 1 -1 T(l-A/4) 

-1 1 1 T(l-A/2) 

-1 -1-1 -T(l - A/2) 

-1-1 1 -T(l-A/2) 

1 -1 -1 -f(l - A/4) 

1 -1 1 -T(l - A/4) 

Thus, the average probability of error for random Manchester coded data 








= y erfc (1 - 4/2)] + i erfc 0 - a/4)] . (2) 

If we define percent asymmetry relative to the inverse of the data rate 
(i.e., T), then a given amount of asymmetry (in seconds) produces the 
same percent asymmetry for either NRZ or Manchester coded data. Thus, 
from Figure 2, we have, for Manchester coded data, that 

Asymmetry = a/4 . (3) 


Using this definition. Figure 3 plots SNR degradation in dB versus percent 
asymmetry for Manchester coded data and E^/Nq= 1.5. For comparison, the 
corresponding results for NRZ obtained in [1] are also illustrated. Note 
that the Manchester coded case always yields a larger SNR degradation for 
a given percent asymmetry. This conclusion is in accordance with our 
previous observations relative to the increased average transition density 
of Manchester coded data relative to NRZ. 

If asymmetry is defined as was done for NRZ but now is relative to 
the half-pulse duration, then 


Asymmetry 


(T/2)'^ - (T/2)~ ^ 2^1 +a/2) - |(1 - A/2) 
(T/2)^ + (T/2)‘ I(l+A/2) + 1(1 -A/2) 


A 

2 • 


(4) 


Note that, for this definition, a given percent asymmetry (i.e., A/2xl00) 
results in twice as much asymmetry (in seconds) for Manchester than it 
does for NRZ. Nevertheless, using this definition. Figure 4 illustrates 
SNR degradation in dB versus percent asymmetry for Manchester coded data 
at E^/Nq= 1.5 dB. The NRZ curve is again included for comparison and is 
identical with that given in Figure 3, Here note that, for small asym- 
metry values, the Manchester code yields a larger SNR degradation than 
NRZ while, for large percent asymmetries, the reverse is true. Thus, in 
drawing conclusions relative to the performance of Manchester versus NRZ 
coded data in the presence of data asymmetry, one must exercise care in 
applying their definitions of data asymmetry. 



Figure 3. SNR Degradation vs. Percent Asymmetry for Manchester and NRZ 

Coded Data 



SNR Degradation (dB) 


T 


Asymmetry = a/ 2 
£3/N(,= 1.5dB 


Manchester/ 


•Asymmetry {%) 

Figure 4. SNR Degradation vs. Percent Asymmetry for Manchester and NRZ 

Coded Data 
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I a„(t;0) = I - nT) , 

n=-oo n"-“ 


where is the polarity (±1) of the jith data bit, and p(t) is the pulse 
shape which, for NRZ data, is defined by 


p(t) = 


A; 0<.t<T 


0 ; otherwise . 


with A denoting its amplitude and T its nominal bit period. In the 
presence of data asymmetry, this same data stream is described by 


nT<t<(n+l)T ifd^=d^^^ = l 


a„{t;n) “ 


nT <.t <.(n+l+n)T if d^ = = 1 


■A; (n+n)T:< t < (n+l+n)T if d^ = d^^l=-l 


■A; (n+n)t ^t.^ (n+1 )T ~ "^^n+l " 


Since the ideal low-pass filter in Figure 1 has transfer function 


H,(f) = 


1 ; -B< f <.B 


0 ; el sewhere 


then the transfer function B^(f;n) of the response of this filter to the 


n^th data bit a^^ltin) is 


d AT [—tt 
n n \ Tr fT 


B„(f;n) 2 <r|b^(t;n)) = 


sin 

^ T* I e 


■B^fxB 


0; elsewhere 


where 




T(l+d^n); 1fd„ = -d„„ 


the corresponding time response b^(t;n) is given by 


b^(t;n) = B (f;n) e 

J_B '' 






2d„A fiBT „ F / 2T rii + 4(1 -d„) 

— cos x (l - |i + ^ 

* Jo L \ 'n 'n , 
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The total response of the filter to the doubly infinite pulse train can 
be expressed as 


y(t;n) = I b (t;n) = 
n=-« 


-Jo T 




“ 2d A firBI' ^ 
r n R/Sinx 


(^) cos X n - ^ + 


2T[n*f(l-d„0 


The first term in (8) is the desired signal whichi for BTq< 1, 
peaks at a value of t such that the argument of the cosine function 
equals zero. Thus, for ^^=1, this peak occurs at [see (6)] 


if dp=d, 


^(1 + n) ; if d^= -dy 


For dQ= -l, the peak occurs at 


, ^(1 +2n) '. if d^ = d^ 


fO + n) ; if dQ= -d^ 


Since, for random data, dg and d^ are each equally likely to take on 
values +1 and -1, then, on the average , the peak of the desired signal 
occurs at t = ^ (1 + n) - t^ , where t^ denotes the desired sampling time. 

The second term in (8) is the intersymbol interference due to 
bandlimiting the input signal. Thus, when sampled at t=t^, the total 
response becomes 


y(tg-,n) = 


2 dg A rTrBT 


( cos jx^(dg + d^) dx 


+ I 

n=-oo 

n/0 


~ 2d A firBT? . f / (l+nd ) pn 

=-00 ^ Jn ^ \ 


where 


= 1 + -^(d -d .,) 
2 ' n n+1 


( 12 ) 
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Alternately, after some trigonometric simplification, 

2 djs A / 1 1 \ 

y(t^;n) = {^Si [trBT(l +ndQ)] + IsjTrBTd - nd^)]} 


where 


+ l -“{j-Si 8T(2n + 1 -nd„^,)] + isi [,BT (2n-l - nd^)]} 

(13) 


Si(y) = 


Sin X 


dx 


(14) 


« 


The total filter response sampled at t^ is then hardlimited to 
determine the estimate d^ of the zero th bit polarity d^. The probability 
of error in making this decision is given by 


= lprob{j-(t^;n) + n,(t^) < 

+ i Prob |y(t^;n) + n,(t^) > 0|d^=.i| 


(15) 


where n-,{t) is the low-pass filtered noise process which is zero mean 

'2 

Gaussian with variance a =NqB, and the overbar denotes statistical 
averaging over the joint distribution of the doubly infinite data 
sequence {d^}; n / 0. As an approximation to (15), we shal l assume 
that the contribution to the total intersymbol interference caused 
by bits further away than N bit intervals from the bit being detected 
(dg) is negligible. Stated another way, we shall consider only the 
intersymbol interference effects of the N preceding and N subsequent 
bits on the bit under detection, where the value of N typically depends 
on the product BT. For BT^ 1 , a value of N = 5 is sufficient. 

Making the above simpl if ication in (15) produces the result 
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I ^ erfc 
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i = l 


erfc 
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where E = An is the energy per bit, 


r A 2 
erfc X = 




e"^ dy , 
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From (13), we find that 

1 i r* V r^*T* / l. \T c r* J r r\ -r- / , ^ ( "^ ) 








^0=1 


|s^. [ir BT (1 + n)] + Si [t: BT (1 - d^^^ n)]| 


(i) 


b_y \t^;n) 




^|si [7TBT(2n + l -d^?^yn)]-Si [7rBT(2n-l -d^^''^n)]} ; 

n^,-l 

^ |si [uBT(l +n)] - Si [ttBT (3 + dj^''^)]} (19) 


d(i) 


and 

^0 


A 


b^^\ts;n) 


A 


d^=-l 


do=-i 


(20a) 


= ~ |si [tt BT (1 - n)] + Si [it BT (1 - d^^^■ n)]| 

d(i) 

= - -Tlj- |si [tt BT (2n + 1 - d^^^^^)] - Si [tt BT (2n - 1 - dj 

n= 0,-1 (20b) 


The superscript i on bn(t^;n) and b„(t^;n) refers to the evalua- 
“ u s n s 7 ,- \ ‘a / /,-\ /.-A 


tion of these quantities for the Hh data sequence 
. . . , d_y ^ d/ ^ I d^iV. . . * dj^^) ; i = 1 .2, . . . ,2N+1 . 


v-n ’ 


-N+1 


b.y \ts;n) 
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A 




-^|si [uBT(l-n)] - Si CttBT (3 + d_^^^^)]| .. (20c) 


Substituting (19) and (20) in (18) and simplifying yields the final 
result 

M 


D^^(BT,n) = 


1 


/2BT 


ni-N 0“"'’ ' ^nVl’) s(2n+l-nd„(;,)[ 


d(^U, 


+ d,y>s(2N + l-<'ndj|’>)+ dj^> s(2N + l-ndJ^*')J 


D_^^.(BT,n) 



where 

S(x) ^ ^ Si(TrBTx) . (22) 

Combining (16) and (21), one obtains the average error probability 

performance of the sample detector in the presence of ideal low-pass 

band! imi ting and data asymmetry. This performance is illustrated in 

Figure 2, with BT a parameter and fixed data asymmetry^ namely n = 0.1 

(10%). In this and the remaining numerical illustrations, the number N 

of preceding and succeeding bits whose intersymbol interference is taken 

into account has been set equal to 4. Comparing Figure 2 with the same 

figure in [2], where data asymmetry was absent, we observe a similar 

behavior as a function of BT, namely, the average bit error probability 

in the presence of 10% data asymmetry is again minimum for BT = 0.9. To 

what extent this is true at other values of data asymmetry is demonstrated 

-4 -5 

in Figures 3 and 4 when, for fixed error probabilities of 10 and 10 , 

respectively, the energy-to-noise ratio required to achieve this error 
rate is plotted as a function of percent data asymmetry with BT again 
as a parameter. Relative to the E/Nq values at n = 0, the curves of 
Figures 3 and 4 illustrate the combined degradation due to intersymbol 
interference (bandwidth limiting) and data asymmetry. 
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3.0 PERFORMANCE OF A BAND LIMIT AND SAMPLE DETECTOR IN THE 
PRESENCE OF NRZ DATA ASYMMETRY - D.C. RESTORATION BY 
CAPACITIVE COUPLING 

In [7], it was shown that, when the input signal is coupled 
through a capacitor, and hence the actual d.c. component in the signal 
vanishes, the degradation due to data asymmetry is reduced. The reason 
for this is that the d.c. restoration has the effect of producing an 
artificial shift in the decision threshold away from the shortened 
symbols (the culprits in the degradation due to data asymmetry). We 
will show here that, for the sample detector with bandl imiting, a 
similar reduction in data asymmetry degradation occurs when the input 
signal is capacitively coupled. 

Following the approach taken in [7], let D denote the transition 
density of the input NRZ data stream, e.g.|( for random data, D=0.5. 

Then the average value of this data stream is nDA. Alternately, the 
positive (longer) signal level is now A(1 - nD) and the negative (shorter) 
signal level is now -A(l + nD). These levels replace A and -A, respec- 
tively in (3). The impact of this on, for example, the time response 
bj^(t;n) of (7) is that A is replaced by A(l-d^nD). Thus, 

2d (l-d nD)A 
b„(t;n) = —5 L 

II -TT 



and, since the argument of the cosine function is unchanged by the capa- 
citive coupl ing, the proper sampl ing time is still t^ - 1-(1 +fi) . Using 
(23) rather than (7), one should have little difficulty in directly , 
seeing how to modify D^.(BT,n) and D^j.(BT,n) of (21) for use in the 
error probability expression of (16). Specifically, we now have that 


D^^.(BT.n) = 


/2BT 



d J '■ ) (l - . 0 d„< ’ )) - (l - n D d„(,V)] S (2N+1 - n d^'J/ ; 
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d(’%, 


+ d_<„* > (l - n D d_y >) S (2N+1 + n d_^ >) + dj ' ’ (l - „ D d^' ' >) S (2N+1 - n dj;,>) 


D^^^(BT,n) = - 


/2FT 



dn<’)(l - no «„<'■)) - d„y,)(l - nDd„y>)] s(2n+l - nd^y,’) 


(i). 


+ d_y>(l -nDd_y>) s(2N+l + nd_y>) + d|’'>(l - n Dd^ *>) s(2N+l - n d,yy)| 


(24) 


with S(x) still defined in (22) 

Analogous to Figure 2, the error probability performance of the 
filter-sample detector with data asymmetry and d.c. restoration by 
capacitive coupling is illustrated in Figure 5 with BT once again a 
parameter, n = 0.1, and transition density D=0.5. Comparing Figure 5 
with Figure 2, v/e observe that, for a fixed error probability, and data 
asymmetry, the required E/Nq at the optimum BT=0.9 is virtually inde- 
pendent of whether d.c. restoration is present or not. On the other 
hand , for other values of BT, e.g., 0.5, the required E/Ng can be much 
smaller when d.c. restoration is employed. Finally, Figures 6 and 7 
demonstrate the corresponding results to Figures 3 and 4 for the d.c. 
restoration case. Here, the value BT = 0.9 appears to be optimum for 
all data asymmetries in the sense of minimum required E/Ng at a fixed 

Pe- 

4.0 CONCLUSIONS 

The combined degrading effects of intersymboT interference due 
to bandl imiting and data asymmetry have been assessed for a filter-sample 
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type detector. For small values of asymmetry (less than about 105^), the 
optimum filter bandwidth- bit time product remains equal to 0.9 and the 
corresponding amount of energy-to-noise ratio degradation at a fixed 
error rate is virtually independent of whether or not d.c. restoration 
is present. In particular. Tables 1 and 2 summarize typical results for 
the two cases when BT=0.9. The quantity AE/N^ represents the additional 
E/Nq required at a given value of data asymmetry relative to its value 
at n = 0. 

On the other hand, for values of BT other than the optimum, d.c. 
restoration has a beneficial effect in reducing energy-to-noise ratio 
degradation due to data asymmetry. 

The other conclusion which can be reached by comparing the numer- 
ical results of this appendix with those in [8] is that the sample detector 
is much less sensitive to data asymmetry than the integrate-and-dump 
detector. This is not surprising when one realizes that the degradation 
due to data asymmetry for the integrate-and-dump detector comes about 
because of a reduction in the signal energy as a result of integrating 
only over a fraction of the total bit interval. By comparison, the 
reduction in the peak of the filter response at, the sampling instant 
due to data asymmetry should have a considerably lesser effect. In 
fact, in the absence of band! imi ting, the integrate-and-dump detector 
would still yield considerable E/Nq degradation (see [3-7]) due to data 
asymmetry, whereas the sample detector would show none. Of course, the 
sample detector without bandlimiting would produce infinite noise power 
and thus, on an absolute E/Nq basis, would be far worse than the integrate- 
and-dump. 
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DETECTION OF ASYMMETRIC NRZ DATA USING A 
GATED INTE6RATE-AND-DUMP FILTER 

by 

Marvin K. Simon 


1.0 INTRODUCTION 

It is well known that the optimum detector (in the sense of 
minimum error probability) of symmetric NRZ data is an integrate-and- 
dump filter whose integration time extends over the entire data symbol 
interval. When the input data stream possesses asymmetry, then the 
ideal inteprate-and-dump filter is no longer the optimum detector; thus, 
it is possible that an alternate (possibly simpler to implement) detector 
could yiel^. superior performance. In Appendix C (see also [1]), the per- 
formance of a filter-sampler type data detector when the input was asym- 
metric NRZ data was investigated. It was shown there that, by optimiz- 
ing the ratio of filter bandwidth to data rate, one could, for sufficient 
data asymmetry, outperform the ideal integrate-and-dump filter. An 
alternate approach is to use a gated integrate-and-dump filter (Figure 1) 
where the gating interval can be selected to minimize tne error proba- 
bility for a given amount of data asymmetry. Such a selection results 
in the best achievable error probability performance using a gated 
integrate-and-dump at each value of asymmetry. Alternately, the gating 
interval may be fixed at a value which represents a compromise between 
the additional SNR degradation (over that of the ideal integrate-and- 
dump) which can be tolerated when no asymmetry exists and the SNR 
improvement obtained at the maximum value of data asymmetry expected. 

In either case, the selection of a gated integrate-and-dump filter has 
the additional implementation advantage of operation at high speed with 
smaller circuit losses since the constraint on its switching times at 
the symbol transition instants can now be considerably relaxed. 

The purpose of this appendix is to present the performance of 
the gated integrate-and-dump filter as a function of data asymmetry and 
gating interval and demonstrate the relation between these two param- 
eters which optimizes this performance. As in previous studies, the 
results will be given for both the case of no d.c. restoration and the 
case of d.c. restoration by capacitive coupling. The numerical results 
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presented permit the above-mentioned tradeoff between additional SNR 
degradation at no asymmetry and SNR improvement at maximum asymmetry 
to be made. 


2.0 PERFORMANCE OF GATED INTEGRATE-AND-DUMP DETECTOR IN THE 

PRESENCE OF NRZ DATA ASYMMETRY - NO D.C. RESTORATION 

When the input data is random with equiprobable symbols, then 
the gated integrate-and-dump (GI&D) depends, in general, on the polarity 
of the symbol over which it is integrating and that of the preceding 
and succeeding symbols. Thus, one must compute these outputs for each 
of the eight possible three-symbol sequences (analogous to the approach 
taken in [2] for the ideal integrate-and-dump) and their corresponding 
conditional error probabilities. Then averaging these conditional 
error probabilities over the equal probabilities of the eight three- 
symbol sequences gives the average error probability performance of 
the GI&D detector. 

When asymmetry is present, the GI&D outputs can be grouped into 
four distinct categories. Arbitrarily assuming that the +1 symbols 
are elongated by nT (relative to their nominal value of T sec) and 
the -1 symbols are shortened by the same amount, then for all four 
three-symbol sequences whose middle symbol is +1, the mean GI&D output 
X is 

X = Jl^Td -2e) , (1) 

where E^ is the symbol energy, e is the fractional (relative to T) gate 
interval at each end of the symbol, and n is the fractional data asym- 
metry. Similarly, for the three-symbol sequence in which all symbols 
are -1, the GI&D output is 

X = . (2) 


The remaining three possible sequences are split up as follows. For 
the sequence whose middle symbol is -1 and both adjacent symbols are 
+1 , we get 



-jE7T(l-2n + 2e) : 

(3) 


1 


while for the two sequences whose middle symbol is -1 and only one 
adjacent symbol is +1, we get 


X = 


-jEg T ( 1 - n) ; 0 <- e ^ 


-JEJ (l-2e) ; 


2 < E < 2 


( 4 ) 


For all eight possible three-symbol sequences, the standard 
deviation a of the GI&D output is easily determined to be 

A 




^0-2e) 


(5) 


Finally, then, the average error probability is given by 


1 erfc 


'Pi 0 . 


( 6 ) 


where 


erfc X = 




exp (-t*^) dt 


(7) 




and the overbar denotes statistical averaging over the eight equi- 
probable three-symbol data sequences. Substituting equations (1) 
through (5) into (6) and the fact that each three-symbol sequence 
occurs with probability 1/8 gives the final result: 
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16 


erfc 


(.1 - 2 e ) 


+ ^ erfc 


+ 16 erfc 


/ I - 2n + 2e \ 

^ PV^ / 



\PT^i 


0 ^ E i 






erfc 


(l-2e) 


2 - ^ ^2 


( 8 ) 


Figure 2 is a plot of Pp versus e with n as a parameter and E /N„= 9.6 dB 
(corresponding to P^=10”^ when E=n=0). We observe from this figure 
that, for a given value of data asymmetry n, P^ is minimized by choosing 
E = n/2. Figure 3 is an illustration of the symbol energy-to-noise 
ratio (in dB) required to achieve an average error rate of 10’^ in the 
presence of data asymmetry. The curve labeled e=0 corresponds to the 
performance of the ideal I&D (also see Figure 3 of [2]). The remaining 
curves indicate a constant E^/Nq for values of n<.2e in accordance with 
the second equation of (8) followed by an increase in E^/Nq with n as 
required by the first equation of (8). Note that each of these curves 
cross the e=0 curve at some value of n, say tiq, which means that, for 
ti>tiq, the GI&D outperforms the ideal I&D in the sense of requiring less 
Eg/No for a given average probability of error. The dashed curve in 
Figure 3 represents the performance corresponding to selecting e=n/2 
at each value of n and is thus the best achievable with the GI&D. 

3.0 PERFORMANCE OF GATED INTE6RATE-AND-DUMP DETECTOR IN THE 

PRESENCE OF NRZ DATA ASYMMETRY - D.C. RESTORATION BY 

CAPACITIVE COUPLING 

When the GI&D is capacitively coupled, then d.c. restoration of 
the input data stream takes place. For NRZ data of transition density D 
(e.g., D-0.5 for random data), the average value of the data stream is 
nDjEg/T . Thus, after d.c. restoration, the positive (elongated) symbol 
now has a level ^E^/T (1 - nD) and the negative (shortened) symbol has 
the level -Je^/T (1 + nD) . The impact of this unbalance in positive and 
negative symbol levels on the output of GI&D for the eight possible 
three-symbol sequences is summarized below. 

X = /E~T~ (1 - 2e)(l - riD) ; all four three-symbol sequences 

^ whose middle symbol is +1 (9) 

X = -JE T (1 - 2e)(l + nD) ; sequence in which all three 

symbols are -1 (10) 


-jEgT [l-2n+nD + 20(l -nD)] ; 0 < s <. | 

-JVr(l - 2e)(H-nD); 2 ® - 1 

sequence in which middle symbol 
is -1 and both adjacent symbols 
are +1 (11) 
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-^E^f [1 -n +nD- 2enD] ; 

X = 

-JEJ- (l-2c)(l + nD) ; 

two sequences in which middle 
symbol is -1 and one adjacent 
symbol is +1 (12) 

Since the standard deviation of the GI&D output is unaffected by the 
d.c. restoration, (5) still applies and thus the average probability 
of error is given by 



ierfc y^(l -2e) (1 -hD) + erfc ^^(1 - 2e) (1 + nD) 


+ ig erfc 


+ g erfc 


/^s /I - 2n+ nD+2e(l -nD)\ 

/rrrr /_ 

rs /I - n + nD - 2e nD\ . 

>oV /rrzT ij’- 


(1 - 2e) (1 - nD) + ^ erfc 


0 < e < 


(1 - 2e) (1 + nD) ; 


2 < e < 2 


Figure 4 again illustrates versus e with n as a parameter and 
Es/Nq= 9.6 dB, where is now computed from (13). We observe from 
this figure that, for a given value of data asymmetry, there exists 
a value of e which minimizes P^; however, unlike Figure 2, this value 
of e, namely is not equal to n/2. Figure 5 is the analogous 

figure to Figure 3 when d.c. restoration is present. Again, the dashed 
curve corresponds to e = which represents the best achievable per- 
formance using the GI&D as a data detector. Comparing Figure 5 with 
Figure 3, we observe the considerable reduction in SNR degradation due 
to data asymmetry when d.c. restoration is employed. This improvement 
is ancilogous to that achieved when other types of data detectors are 
used [2]. 



Figure 4. Average Error Probability versus Gate Interval at Symbol Edge with Data 
Asymmetry as a Parameter - D.C. Restoration by Capacitive Coupling 
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4.0 CONCLUSIONS 

The use of a gated Integra te-and-dump filter for detection of 
asymmetric NRZ data can, depending upon the amount of data asymmetry 
present, produce significant improvement in SNR degradation due to 
asymmetry relative to that of an ideal integrate-and-dump filter. 

Such a filter has the further implementation advantage of allowing 
finite switching times at the symbol transition time instants. 

Since, for no data asymmetry, the performance of the gated integrate- 
and-dump is inferior to that of the ideal integrate-and-dump {matched 
filter), a tradeoff exists between the two filter types depending 
upon the amount of asymmetry present in the data. Whether or not d.c. 
restoration is present, there exists a value of gate interval for a 
given value of asymmetry which minimizes the average- error probability 
performance of the receiver and thus represents the best achievable 
performance using a gated integrate-and-dump. The particular value 
of gate interval which achieves this minimum does, however, depend 
upon whether d.c. restoration is used or not. 
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TRACKING PERFORMANCE OF UNBALANCED QPSK DEMODULATORS 
PART I - BIPHASE COSTAS LOOP WITH PASSIVE ARM FILTERS 

by 

Marvin K. Simon 
Uaddah K. Alem 

INTRODUCTION 

Unbalanced quadriphase-shift-keying (QPSK) is an attractive means 
for transmitting two digital data streams which have different average 
powers. The two data streams are not constrained to have identical data 
rates nor must they have the same data format; e.g.* one might be an NRZ 
sequence and the other a Manchester code. In fact, it is the difference 
in data rates which causes the unbalance of power when it is desired to 
have symbol energies and therefore error rates on the tvio channels within 
the same order of magnitude. 

Previous results [1,2] have indicated that when the unbalanced 
power ratio is large, e.g., approximately 4:1 or greater, a biphase 
Costas loop is a more efficient demodulator than a fourth-power tracking 
loop. These results, however, accounted only for the filtering effect 
produced by the loop's two arm filters on the equivalent additive noise 
perturbing the loop. When the bandwidth of these filters is selected on 
the order of the data rate, as is typical of optimum Costas loop design 
[3,4], the filtering degra< ations of the data modulations themselves and 
the cross-modulation noise produced by their multiplication in the loop 
often cannot be neglected. 

The purpose of this report is to incorporate these additional 
filtering effects into the analysis of a biphase Costas loop demodulator 
of unbalanced QPSK. Many of the results obtained herein are in the form 
of closed-form expressions which can easily be evaluated numerically for 
design and performance prediction purposes. The generality of the results 
enables them to be applied to a wide variety of applications such as the 
performance of a subcarrier tracking loop for the three-channel Space 
Shuttle Orbiter (SSO) Ku-band return link [5]. 

SYSTEM MODEL AND LOOP EQUATION OF OPERATION 

Consider the Costas loop illustrated in Figure 1 when input x{t) 
is an unbalanced QPSK signal plus noise, i.e.. 




















x(t) = cos «>{t) + /2P, m,(t) sin ^(t) + n,.(t) , (1) 


2 "' 2 ' 


where 4>(t) - ixi^t + e(t) with the radian carrier frequency and 
0(t) = 0g + «gt the input phase to be estimated, m,(t) and m,(t) are 


1 


( 2 ) 


independent data modulations, and n^.(t) is the additive channel noise 
which can be expressed in the form of a narrowband process about the 
actual frequency of the observed data, i.e., 

n^.(t) = v^|N^(t) cos 4>{t) - Ng(t) sin <t(t)| . . 

In (2), N (t) and N (t) are approximately statistically independent, 

w 2 ) 

stationary, white Gaussian noise processes with single-sided noise spectral 
density Ng w/Hz, and single-sided bandwidth B^<Wg/2TT. 

The input signal x(t) is demodulated by the quadrature reference 

signals 


r^(t) = sin ${t) 

r^(t) = /2 cos $(t) 

to produce the corresponding phase detector outputs (ignoring second 
harmonic terms): 


(3) 


ejt) = x(t) r (t) = K, K 
s m s I m 


ni 2 (t) - Ng(t)J cos <^(t) 

■ *^1 ^1 

e^(t) = x(t) r^(t) = ^ - N^{t)J sin (p(t) 




v^m^(t) + N^(t) 


cos ^(t) 


(4) 


where (p(t) - 4>(t) - $(t) is the loop phase error. After lowpass filter- 
ing ^^(t) anc 
respecti vely. 


ing £g(t) and £^(t) with arm filters G(s), these signals become. 


ZjCt) = G(p) Ej(t) = K, l<^ 


m 


/P^ m2(t) - f)g(t), 


cos 9 ?{t) 


- 

I m 


/P^ tfi^{t) + f)^(t) 


sin <;o(t) (5a) 
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2 ^(t) = G(p) e^(t) = K, jj*^ra 2 (t) - Nj(t)J sin Ip(t) 

+ K, [#j" m,(t) + cos <p{t) (5b) 

where the "hats" denote filtering of the corresponding signals; e.g. , 
m^(t) ^ G{p)nii(t)- Thus, the output of the'third multiplier is the 
dynamic error signal 


Zo(t) S ,^(t) z^(t) 

1^2 2 

= — L — Pg (t) - P-j m^^ (t)J sin 2<p{t) 

+ 2 /P^ nii{t) m 2 (t) cos 2(p{t) 

* Ifs^ (*) - \^(t) - 2-^ "l2{t) 

- 2*^ m|(t) sin 2i;i?(t) 

2^ m^{t) N^(t) - 2 /p 7 m^(t) n^(t) 

- 2 Ngltjj cos 2<^(t)j . 

The instantaneous frequency of the VCO output is related to ZQ(t) by 

= Ky [F(p) zjt)] + 


dt » — • “0 

and hence the stochastic equation of loop operation becomes 

= 2n„ _ kf(p) 


dt 


P2i"2^(t) - Pim/(t) 


sin 2(p{t) 


+ 2/P^P^ iTii(t) rriglt) cos 2cpit) 

+ V2[t,2(^(t)] 


(6) 


(7) 


( 8 ) 


0 


where 


5 


V2[t,2(^{t)] = 


- N^^(t) - 2.^ni2(t) N^(t) 


2^ m^(t) Nj.(t)J sin 2(p{t) 


N^(t) - 2^/^ ri1^(t) NgCt) 


- 2 N^(t) NgCtjj cos 2(p{t) (9) 


A ? 7 

and K - V *^V' previously done [3,4] in problems of this 


type, we shall now decompose each of the signal terms in (8) into its 
mean value plus the variation about this mean, e.g., 


m^^ (t) sin 2{/?(t) = <nii^ (t)> sin 2<p{t) + 


m. 


(t) - (t)>J sin Z<p(t) 

no) 


where the overbar denotes statistical expectation and < > denotes time 
averaging.* It is easily sho^vn that 


(t)> ^ D|^ 


Sn,.(f)lG(j2Trf)r df ; k=l,2 


( 11 ) 


and since m-j(t) and m 2 (t) are independent, 


m^(t) iTi 2 (t) = 0 


02 ) 


In (11), Smj^(f) denotes the power spectral density of m|^(t); k=l,2; and 
lG(j2Trf) |^ is the magnitude .squared of the arm filter transfer function. 


Rewriting (8) using the above decompositions gives 

— - KF(p)|[P 2 D 2 -P^D^] sin 2cp{t) + r)^{t,2cp)^ (13) 


where the total equivalent additive noise n^{t,2(p) is given by 

A 


Og(t,2(^) = y^{t,2(p) - P^ n^(t) sin 2(^(t) + P 2 n 2 (t) sin 2(^(t) 

+ 2/P| P 2 n^2(^) 2</?(t) (14) 


The additional time averaging is necessary since ni] (t) and m 2 (t) 
are cyclostationary processes. Alternately, one may include a random phase 
in these two processes and statistically average them over the uniform 
distribution of this phase [6]. 




6 




m 


o 

o 

.1 

( 

© 


' • 

i 

i 

I 

i 

i 




with the self-noise nj^(t); k = l,2; and the cross-modulation noise n^g^^) 
defined by ‘ ' 


nk(t) - tn^(t) - (t)> ; k = 1,2 


n^2(t) = m^(t) ni2(t) 


STATISTICAL CHARACTERIZATION OF THE EQUIVALENT ADDITIVE NOISE 

The self-noise and cross-modulation noise processes all have zero 
mean and a continuous power spectral density component. Since the band- 
width of these processes is very wide with respect to the loop bandwidth, 
it is sufficient to find for each one only the power spectral density at 
the origin which, when multiplied by the loop bandwidth, gives the con- 
tribution to the total noise power of that modulation component. From 
(14), the autocorrelation of T\^{t,2(p) is easily shown to be 

Rg(T,2<;^J) ^ <ng(t,2</j) ng(t+t,2^)> 

= P^^sin^ 2(;pRn^(T) + p| sin^ 2<^ Rn2(T) 

+ 4 P-j P 2 cos^ 2^Rni2('f) 

where 

(t) i <iS (t) m.(t+T)> = [ S~ (f)|8(j2irf)|^ df ; k=l,2 

k •'-00 K 

Rnj^(T) = <n(^(t) n|^(t+T)> ; k=l,2 

~ *Cn-j2(t) ~ Pm2^^^ 

RjCjit) ^ N^(t) N^(t+T) = Ngit) N^Ct+x) 

= ^ [ lG{j2Trf)|^ e'^^'"‘*'Vdf . (17) 

The equivalent noise power spectral density at the origin' is then 


7 


# 


m 




n(2(p) = 2 R^{t. 2(^) dT . 

■Loo 


(18) 


Substituting (16) together with (17) into (18) and using Parseval's 
Theorem yields 

Ng(2(^) = sin^ 2<^ Sn^(O) + P^^ sin^ 2<;^? Sn2(0) 

^ ^00 

+ 8 P^ ?2 COS^ 2<^ Sn,^(f) Sm2(f) |G(j2Trf)Pdf 

00 

«00 *00 

+ 4{p^Nq Sm^(f) lG(j2Trf)r df + P 2 NQ S„,2(f) |G(j2Trf)t'^ df 

' ^*.00 'fvOO 

lOO 

|G(j2irf)|'^ df} . (19) 

^—00 ' 


N(f f 


2 i 


where [3] 


' L ^ i K = 1.2 

with = 1/Tj^ the data rate of the digital modulation m|^(t); k=l,2. 
In [3], it was shown that, for cases of practical interest, the effect 
of the self-noise power on loop tracking performance was negligible. 
Thus, for all practical purposes, we may concern ourselves only with 
N^(2(p) evaluated at^=0, i.e.. 




(p-0 


= 8 P, 2^0,2 - 4 P,NoD,Kd,^P2N„D2Kd2 


(20) 


where Dj^; k=l,2; is defined in (11), 


12 


^m/^^ ^012^^) |G(j27rf)r df 


4 Sm|<(f) |G(j 2 rf)r df _ 

/.“ %fc(f) |G(j2rf)l^ df ’ 
/” |S(j2rf)l'’ df 


k = l,2 


^ •'-TO 


r 


II |G(J2’tf)r df 


(21) 


# 
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m 




1^1 




# 


and 


B. = [ lG(j2Tif)|^ df 

^-00 

is the two-sided arm filter lowpass bandwidth. 


( 22 ) 


EFFECTIVE LOOP SIGNAL-TO-NOISE RATIO AND PHASE JITTER 

From the loop equation of operation given in (13), we can see that 
the effective loop signal-to-noise ratio is given by 

,2 


(P2O2-P1O1)' 




(23) 


where is the single-sided loop bandwidth. Substituting (20) into (23) 
and simplifying results in 


^e " 4 '^L 


(24) 


where p = P^/N^B^ with the total power P^ = ^1 + ^ 2 * 
"squaring loss" which is given by 


(P2D2- P^Di)' 


¥i 2 ^ 1 ^ 2 '^! 2 

P^DiKn + P2D2KD. + -T^ Kl + ^ ^ 


(25) 


N, 


Defining the modulation indices, 




(26) 


then (25) simplifies to 


^^2^2 ” '^T^l ^ 


n,D,Ko^ + + 2n,n2 ( 1 ,^) 0 


12 


Ordinarily, to optimize tracking performance, i.e., maximize p^, 
the arm filter bandwidth is selected relative to the higher of the two 
data rates. Thus, arbitrarily assuming that i ?2 > i?i > we rewrite (27) in 
the final form 
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m 


(h2D2 " *^101 )' 
TJb. 


n,D,KD^ + '^L ^"l "2 ’*12(12^) 


(28) 


where Rj 2 ■ ’ 

What remains is to characterize the tracking phase jitter perform- 
ance. In the linear region of operation, the variance of the loop phase 
error Zcp is given by 


1 


Z(p 


4 

P s] 


(29) 


Since the demodulation reference signals of (3) are at rather than 2 (Oq, 
then the tracking jitter on the output data streams is 


<P 


1 2 
4 ""Zip 


p5, 


(30) 




We conclude then that characterization of the loop squaring loss in terms 
of the system parameters is sufficient for predicting the loop's tracking 
performance. 

SQUARING LOSS EVALUATION FOR PRACTICAL FILTERS AND DATA MODULATION SPECTRA 

To graphically illustrate the theory previously developed, it would 
be convenient to obtain closed-form expressions for the parameters defined 
in equations (11), (21), and (22) for practical filters and data modulation 
formats. Indeed, such is possible for certain special cases of interest 
as we shall now illustrate. Let us consider a simple RC filter with 3 dB 
cutoff frequency f^ for the Costas loop arm filters. Then, 


m 


|G{j2TTf)!' 


1 + 


2 • 


(31) 


Substituting (31) into (22) and carrying out the integration gives the 
relationship between two-sided noise bandwidth and cutoff frequency, 
namely. 


B. - 


(32) 


10 


Also, from (21), substituting the square of (31) into the numerator of 
the expression for and recognizing that the denominator is merely 
equivalent to of (32) gives, upon integration, 

\ 2 • '' 

Evaluation of D|^, k = l,2, and 0^2 requires that we further specify 
the forms of the two data modulations. For NRZ data with power spectral 
density 


= T. 


sin^ irf Tj 


it has been previously shown [3] that 

Dj^ = 1 - 20 tl “ 6xp (“2B^/i?j^)] j k=l,2 

3 - (3 + 2B./fl(.) exp 

Kp = — , !_!S ; k = 1 ,2 

^ ’ - 2B^ C’ - 

For Manchester codes, the power spectral density is 

sin^ (irf T. /2) 

Sm.(f) = T;, — ; k=l,2. 

Recognizing the trigonometric identity 

.4 212 

sin X = sin X - sin 2x , 

then (36) can be rewritten as 


sin^ (TTf T|^/2) sin^ uf T 


/ ‘ i,\ ^ V " I I u 

S„ (f) = 2 U) — \ 

k \2/ (xfy2)^ 


T. ; k = 1,2 

<"”k> 


which spectrally looks like the difference of two NRZ spectra of different 
data rates. In view of this, it is obvious from (21) that, for Manchester 
codes, 

ORIGINAL PAGE IS 
OF POOR QUALHia 


n MANCH . . NRZ 
^ ^k 


_ n • k = 1 2 

^k=2/?k k > k 1 ,2 


„ MANCH _ 1 I o „ NRZ .. NRZ 

'^D,. ■ . MANCH 1^ ^k '^Dr 


MANCH 


^k 


n NRZ „ NRZ 
V^^k " ^k 


k^ 1,2 . 


( 39 ) 


Substituting (35) into (39) and simplifying yields the following results 
for Manchester codes: 


= 1 




^ b ;^" [3-4 exp + exp (-2B ./flj^)] ; k = 1 ,2 

1 k 

9 - 4(3 + 8./^’,^) exp (-8.^) + (3 + 2B./i?j^) exp (-28.//?,^) 
3 - 4 exp + exp (-2B^./^j^) ~ 


1 - 




k = 1,2 . 


(40) 


All that remains is to derive expressions for the cross-modulation 
coefficient 0^2 corresponding to the various combinations of modulation 
formats and relative data rates. Starting out with the simplest case 
(at least algebraically) , we shall assume that rai(t) and m 2 (t) are both 
NRZ data streams with respective bit rates and /?2 > i?i . Then, 

from (21), (31), and (34), we have that 

/2 

1 


'12 


T^T2 


.2 .2 

■” sin TTfT,| sin irfl 2 

1 

f 

J + (f/fc> - 


df . 


(41) 


Using a partial fraction expansion followed by considerable integral 
evaluations and manipulations, we arrive at the final answer, namely. 
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+ 4- To Cr) ^ T 5 - (5 + 2B./i?J exp (-2B./i?J 

® 2 V'W (B,/^o) 


- (5 + 2B./i?i) exp {-2B./i?^) 


i[5.2(^-.l) B./B^] exp [- 2 (^- 1 ) B./b^ 




m^(t) is NRZj jn 2 (t) is NRZ . 


If m-j(t) is now a Manchester code while ni 2 (t) is still NRZ, then using 
(36) in addition to (21), (31), and (34) gives 

j, i-oo sin^ trf Tt/ 2 sin^ To F i 

D = — ^ L , £ ! ^ df . 


Once again using the trigonometric identity of (37), D -|2 of (43) simplifies 


2.2 r- 

2 r“ sin Trf T -|/2 sin itf T2 1 

(T.,/ 2 ) T2 f"^ b + (T/fc)^- 


4 T T 

TT T^T2 


2 2 

“ sin irf T.| sin trf T 2 




Comparing (44) with (41), we immediately arrive at the following results; 


MANCH 

NRZ 


MANCH 

NRZ 


■ ^^12 \r^=ZR^ ^ H2 ’ h - 


“12'"' • - 2 »ir B,= 2 B - “12 - • ■ 

«r-2 
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Letting m^(t) now be NRZ and ni 2 (t) a Manchester code, we get the remaining 
result: 


NRZ 

MANCH 


'12 


NRZ 

n NRZ . 


* •■^2 ~'^1 •* 


Finally, when both m^(t) and m' 2 (t) are Manchester codes, we have 


4 4 

16 r“ sin irf T^/2 sin ir.f T 2/2 


16 r 


1 + (f/fj,) 


To simplify (47) in terms of previously given results, we require the 
trigonometric identity: 

4 4 2 2 1.2 2 

sin ax sin bx - sin ax sin bx - ^ sin 2 ax sin bx 

12 2 12 2 
- ^ sin ax sin 2 bx + jg sin 2 ax sin 2 bx . 


Making use of (48) in (47), we see, for example, that 

MANCH NRZ NRZ 

n MANCH , n NRZ on NRZ 

*^12 “ ^^2 Ry2R^~ i?^=2i?^ 

i?2=2i?2 




MANCH 

MANCH 


NRZ 


NRZ 

- 2 D- 

i?^= 2 i?| ^nz 

^ 2 = 2^2 ■ 

2 D 

^Ui 2 

i? 2 = 2 i ?2 






2i?i >. /?2 — ^1 


Expressions for D.j 2 , corresponding to (45), (46), and (49), in terms of 
system parameters can be obtained by substituting (42) into these equa- 
tions and simplifying. For example, corresponding to the first of the 
two cases in (45), we get 
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D,2 = Tj (1 - Tj/T,) - STj (Tj/T,) 








15 - 3(5 + 2B./^2^ (-2B./i?2) 

- 4(5 + B.//?^) exp (-B./i?^) 

+ (5+2B.//?^) exp 


+ 2 

1 

■ 2 
+ 2 


5+2 


r- .T 

5+2 


fe- ’) ^ 


e- ') 


/Br 


iBr 


exp 


/T, 

2 


exp [" 2(27 ■ ^) 


fe- 0 




- T T 

5 + 2 ( 2 ^+’) exp 

5 . 2 (^. 1 ) B./il^] «P 


whereas for the second case in (45), 

rr, -,/T,x2 


•^12 " ^2 


_i i/_i\ 
Iz ' ^V 2 / 


1 + 


1 /^ 

At,I 


- 212(1 -■of') 2 

2 \ {B./B^r 


i ?2 ^ 2 i?i ; 


1 


(50) 


1 ' (B./^ 2 ^' 


15 - 3(5 + 2B./i?2) exp (-2B./i?2) 
- 4(5 + B^/?^) exp (-B./^^) 

+ (5+2B.//?^) exp (-2B./i?|) 


+ 2 



) ^^''2 

exp 

’2T2 

-)Bi/R 

2 

1 1 

■ 2| 

r 1 



■ 

- /^i \ 



I 5 + 2 

L ' 



exp 

:Kii- ’) 

B./B^ 





(’ ^ 2 T 2 , 



r / ■^1 



+ 2 

5+ 2| 

) S/''2 

exp 


-) B./B 

1 

2 

1 

" 2 


/T, V 

. 


-ft-') 

_ 


5+ 2| 

fe"') 


exp 

Bi/flj 



Finally, when m^(t) is NRZ, ni 2 (t) is Manchester, and i?-| , then corre- 
sponding to (46) we have 

Di2 = ^ 




+ g T2(T2/T-|) 


(B./i?2)' 


15 - 3(5 + ZE./R^) exp (-28.//?^) 


- 4(5 + B^//? 2 ) exp (-B^/i? 2 ) 

+ (5 + 2B./i?2) exp (-2B./i?2) 

. 2[b.2(i -^) exp [-2(1 -^) B/.,] 


1 

■ 2 

[5 + 2 I 


I 

exp 


B./fl,_ 


r ' 

/ h 

\ 

— j 


\ 

+ 2 

5 + 21 


j 

1 exp 



1 

p- 

/ 0 N 

— I 


r / 

“1 

1 

" 2 

5 + 2| 


I 

exp 


Bi/ff,_ 




i?2 ^ 

• 

(52) 




Obtaining a detailed closed-form expression corresponding to (49), i.e., 
the case of m^(t) and m 2 (t) both Manchester coded streams, is left as an 
exercise for the reader. 


NUMERICAL RESULTS 

Illustrated in Figures 2 and 3 is the behavior of the normalized 
cross-modulation distortion D^2/^2 ^ function of the ratio of the two 

data rates with the ratio of two-sided arm filter noise bandwidth to high 
data rate as a parameter. In Figure 2, the results are for the case where 
both m^(t) and m 2 (t) are NRZ data streams and are thus obtained from (42). 
The results in Figure 3 are obtained from (50) and correspond to the case 
where m|(t) is NRZ and m 2 (t) is a Manchester code, and i ?2 Also 

illustrated as dotted curves are the limiting curves which ignore the 
bandl imi ting effects of the arm filters. These curves which are obtained 






nil(t) is Manchester code at rate 1/T^ 
ni 2 (t) is NRZ at rate 1/T2 


Cross-Modulation Distortion vs. Ratio 
of Data Rates; Ratio of Arm Filter 
Bandwidth to High Data Rate is Parameter 
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from (42) and (50) by letting correspond to the results given 

previously [1 ,2]. 

We now turn to numerical evaluation of the loop squaring loss (28) 
which, for a fixed loop bandwidth and ratio of total power-to-noise spec- 
tral density, is a direct measure of mean square tracking jitter. Figures 
4 through 6 illustrate (for fixed and 112 ) the behavior of as a func- 
tion of the ratio of two-sided arm filter noise bandwidth to the higher 
of the two data rates ^^2 with the ratio of data rates i? 2 /-ff-| and PjT2/^0 
as parameters. Assuming Py/N^ to be fixed, then ther variation of squaring 
loss with PyT 2 /Np directly reflects the effect of changing the high data 
rate i? 2 * Furthermore, at low values of B^//? 2 , we observe from Figures 5 
and 6 that additional interesting peaks and valleys of the squaring loss 
characteristic occur. These extrema represent tradeoffs between SxS 
distortion and cross-modulation noise or SxN power, depending on which 
of the latter dominates the total noise. 

The numerical evaluation of the tracking jitter, for a fixed ratio 
of arm filter noise bandwidth to loop noise bandwidth (B^-/B|^) [(30)], is 
shown in Figures 7 through 14. Figures 7 through 10 illustrate the 
behavior of a^, for fixed n-| , H 2 » 3s a function of the ratio of the arm 
filter bandwidth to the higher data rate i ?2 for all possible signal 
format combinations with 7?2/i?^ and PyT 2 /Ng as parameters. The minimum 
values of for some ratios of represent best design points when 

the combined effect of Nx,N distortion and cross-modulation noise or SxN 
power is minimal. Assuming Py/Ng to be fixed, the variation of with 
PtT 2 /Nq is shown in Figures 11 and 12 for the case when m^(t) and m 2 (t) 
are NRZ and the case when m^(t) is Manchester code and m 2 (t) is NRZ, 
respectively. As is intuitively true, the tracking jitter performance 
improves with the increase of PyT2/^0" Furthermore, for a fixed value 
of B^/7?2» the variations of 0 (p with P 2 /Py “ '^2 shown in figures 13 
and 14 for the two previous cases. The figures show that, as the powers 
of the two data rates become equal (n 2 =l/ 2 ), the loop is not able to 
operate. 


Squarii 

PyT^/N, 

i?2 £.ff' 





15 dB 12.5 dB 





Figure 5b. Squaring Loss vs. Ratio of Arm Filter Bandwidth to High Data Rate 
PtT2/Nq and R2/R'\ are parameters; m,(t) is Manchester code, mr,(t) 
is NRZ; ' 














Squaring Loss vs. Ratio of Arm Filter Bandwidth to High Data Rate 
PTt 2 /Ng and parameters; my(t) is Manchester code, m^Ct) 
is NRZ; 27?-, 






Figure 7. Tracking Jitter Standard Deviation versus the Ratio of Ann Filter 
Bandwidth to High Data Rate; PyTo/Nn and Ro/R-i are parameters; 
m^(t) is NRZ, m 2 (t) is NRZ; 
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Figure 8 , Tracking Jitter Standard Deviation versus Ratio of Arm Filter Band- 
width to High Data Rate; PjT2/Nq and are parameters; m-|(t) 

is Manchester code, m2(t) is Manchester code; — ^1- 
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TRACKING PERFORMANCE OF UNBALANCED QPSK DEMODULATORS 
PART II - BIPHASE COSTAS LOOP WITH ACTIVE ARM. FILTERS 

by 

Marvin K. Simon 


INTRODUCTION 

In Part I, we investigated the performance of a biphase Costas 
loop with passive arm filters insofar as its ability to track an unbal- 
anced QPSK type modulation. In a previous Costas loop study for biphase 
modulation [1], it was demonstrated that considerable improvement in 
tracking performance could be obtained by employing active arm filters 
of the Integra te-and-dump type as opposed to passive arm filters. The 
purpose of this part, then, is to investigate whether a similar perform- 
ance improvement can be obtained for an unbalanced QPSK modulation. 

Before going into the details of such an investigation, we point 
out that the use .of active arm filters in a Costas loop carries with it 
the assumption that symbol synchronization of the two data streams is 
known. This assumption, namely, that symbol synchronization be avail- 
able before carrier synchronization, appears to involve one of those 
unstable situations in which carrier acquisition cannot be achieved 
without first having symbol synchronization, and vice versa. If, how- 
ever, the range of frequencies over which the carrier loop must be 
acquired is kept small (on the order of the data rate or less) by doppler 
correction, then the presence of the beat note, which causes the symbol 
synchronization input to be amplitude modulated by the cosine of the 
carrier loop phase error, does not destroy the information carried in 
the symbol transitions from which the synchronization clock is extracted. 
Hence, for a maximum frequency uncertainty less than the data symbol 
rate, the symbol synchronizer can be designed to acquire rapidly as the 
carrier VCO is swept through the region around the zero beat frequency. 
Thus, a coarse estimate of symbol synchronization is obtained before we 
have actually acquired the carrier loop or the lock detector has decided 
to kill the sweep. . 


2 




# 














9 


9 


SYSTEM MODEL AND LOOP EQUATION OF OPERATION 

To avoid the duplication of effort, we shall, in this section, 
draw heavily upon the notation, definitions, and results given in Part I. 
Consider the Costas loop with integrate-and-dump arm filters illustrated 
in Figure 1, where again the input x(t) is an unbalanced QPSK signal plus 
bandlimited Gaussian noise of constant spectral density. Then, by analogy 
with (5) of Part I, the arm filter outputs are now given by 

Zg(t) = Kj^T 2 [/^m 2 (t) - N^(t)] cos (p(t) 

■ *^1 *^m^2 ^*^^1^^^ 

^c^^^ Kj^T 2 [/P^m 2 (t) - Ng(t)] sin (p{t) 

+ Kl Kjj]T 2 [/^m^(t) + Nc(t)] cos c^(t) ; 

(k-l)T2 < t < kT2 , (1) 

where* 

kT2 

ni,-(t) - ni.(t) pp [t- (k-l)Tp] dt ; i = 1 ,2 

’ . '2 J(k-1)T2 ^ ^ ^ 

kT2 ■ 

^,(t) = f N (t)pp[t- (k-l)Tp] dt ' 

^ ^2 J(k-1)T2 ^ ^ 

kTp 

N (t) pp [t - (k-l)Tp] dt ; 

^ '2J(k-l)T2'' ^ ^ 

(k-l)T2 t .<. kTp , (2) 

and p^.(t), i = l,2, is the basic unit power symbol pulse in the J_th data 
stream defined to be nonzero only on the interval (0,T|). 

Forming the product of z (t) and z (t), gives the dynamic error 

signal , 

•k 

We assume here, as in Part I, that Ro>.R-\{T 2 and that the 
integrate-and-dump arm filters are identically designed to accommodate 
the higher data rate signal. 
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|(' 2 2 T 2 / - - “ 

^0^^^ "" " ~" I “ P^<ni]^(t)> sin 2 (p{t) 

2 V ^ - 


+ 2 ^/P^ P2 <iiii (t) m2(t)> cos 2 q){t) 


+ Njt,2(^(t)] j (3) 

where the total equivalent noise Ng[t,2 </?(t)] is defined by 
Ng[t,2<^(t)] = [N5^(t)-.N^?(t)-2/^m2(t)N^(t)-2/^m^ sin2</7(t) 

+ [2v^m2(t)N^,(t)-2/^ni^(t) N^(t)-2N^(t)N3(t)]cos 2 c/ 7 (t) 

+ |P2 ni2^{t) - <m2^(t)> - P-, lii^^(t) - <m^^(t)> | sin 2 (pit) 

+ 2 /P^pp^ (t) iii 2 (t) - <m.j (t) m 2 (t)> cos2(/?(t). (4) 

In (3) and (4), the overbar denotes statistical expectation and <> 
denotes time average. Since the VCO phase estimate e(t) is related to 


ZQ(t) via 


Kw'^(P) 

e(t) = Zgit) , 


then the stochastic equation of loop operation becomes* 

2 p(p = -KF(p)jT2^ P2<^2^(t)> - P-j <m^(^)> sin 2<^ 

+ 2 t|.^P^ <m,(t) iSjit)) cos 2 <p + T| Ng(t,2(^ )j (6) 


where 


Since the digital modulations can be expressed in the form 


mpt) = I a^. ppt- nip ; i = 1,2 


For simplicity of notation, we drop the dependence of 99 on the 
time variable t. 
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where is the polarity (±1) of the iith data symbol in the Uh data 
stream, then assuming that the and 3 ^ 2 '^ independent, we have 

from (2) and (7) that 


<m^ (t) m2(t)> = 0 . (8) 

Furthermore, the third term in (4) represents the self noise of the 
filtered modulations, which has been shown previously [1] to have a 
negligible effect on loop performance. Thus, ignoring these self-noise 
terms and using (8) in (4) and (6) gives the simplified results 

2p^ = -KF(p)|t 2 ^ P 2 <>^ 2 ^(t)> - P-|<m^^(t)> sin 2cp + l 2 ^ Ng(t,2(/5 ) j (9) 

where 

Ng(t,2</7) = [N 3 ^(t)-N^^(t)- 2 ^m; 2 (t)N^(t)- 2 /F^m^(t)N^^ sin 2cp 

+ [2 v/P^ if> 2 (t) - 2/^ m^(t) N^(t) - 2N^(t) N^(t)] cos Zcp 

+ 2 P 2 iTii(t) m 2 (t) cos 2cp . (10) 

Note that the first two terms of the total equivalent noise in (10) repre- 
sent the traditional signal xnoise and noisex noise terms characteristic 
of the error signal in Costas loops. The third term of (10) corresponds 
to the cross-modulation noise of the two data symbol streams. 

Linearizing the loop, we have that 


2 <P = Hj^(p) 


f2<jn/{t)>-P,<fn/(t)> 


where H 2 (s) is the closed loop transfer function defined by 


= 


_P2<>^2^(^)> - P-,<ffiv(t)>_'T2^KF(s) 

•i- rP2<m2^(t)> - P^<m^^(t)>]j2^KF(s) 


Assuming that the loop noise bandwidth is much narrower than the noise 
bandwidth of Ng(t), the variance of 2cp can be approximated by 
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K\ 


^ 277^■.12 


(13) 


where 


2B. 


P2<ni2 (t)> “ Pi<V^t)> 




(14) 


and 


k 2 


<RN (t.t+T)> dx 


(15) 


with Rf^ (t,t+x) denoting the autocorrelation function of N (t), i.e., 
6 ” 


R„ (t,t+t) ^ N„(t)N„(t+T) . 


Nq' " e' '' ’e' 

TRACKING PERFORMANCE - THE CASE OF SYNCHRONIZED SYMBOL CLOCKS 

We assume first that the two data clocks are synchronous in that 
they both have a positive-going pulse starting at t = 0 and their pulse 
repetition frequencies are integrally related, i.e., 


(16) 


4 


= n , 


(17) 


where n>il is an integer. Since from (2), m 2 (t) is the output of a 
matched filter matched to the basic pulse shape of m 2 (t), then regard- 
less of the data formats of m.|(t) and m 2 (t), we have from (7) that 


<m/(t» = 1 


(18) 


On the other hand, evaluation of <m-] (t))> depends quite heavily on the 
data formats of m^(t) and m 2 (t). For example, if m.|(t) is NRZ data and 
m 2 (t) is Manchester coded data, then 

1 ; 0 :S. t T-i 


P^(t) = 


Pjit) = 


1 


0 ; Otherwise 


1 ; 0 < t < 

-1 ; T 2/2 < t < T 2 


(19) 


# 
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f ® 




D 


and 


<ra^(t)> .= 0 


( 20 ) 


If both m^(t) and m 2 (t) are NRZ data, then 


<ni-, (t)> = 1 . 

These results, together with those for the remaining combinations of 
formats for m-|(t) and m 2 (t) , are given' in Table 1 below. 


( 21 ) 


Table 1. Evaluation of <m-|^(t)> ; ~ 


ni2(t) 


m^(t) 


1 

Manchester 

NRZ 

Manchester 

0 ; n^2 

1 ; n even 

; n odd 
n 


1 ; n= 1 

NRZ 

0 

1 

t 


Since we have assumed a linear system in writing (13), then inso- 
far as computing the equivalent noise spectral density N’ , we may set 
2 2 G 

sin 2^ = 0 and cos 2(jo = 1 in (10). Doing so, we get the simplified 

result 

Ng(t,0) = 2/P“ m2(t) fl^(t) - 2/F^nT^(t) 

- 2N^(t) Ng(t) + 2 /P^ m.|(t) m2(t) . (22) 

Taking the autocorrelation function of the nonstationary process Ng(t,0) 
and averaging the resulting expression over t yields 


<R^ (t,t+T)> = 4 




+ <Rj;j fi (t,t+T)> + ffl.lt. t+T» 

c s 1 z ~ 


(23) 


Mti 


where 


N (t,t+T)> 


<ifi2(t) Nj,(t) m2(t+T) f^^(t+r)> 


<R~ 1^ {t,t+T)> = <m^(t) l^g(t) tfi^ (t+i) f^5(t+T)> 

1 s 


<R^ ^ (t,t+x)> ^ <N^{t) Ng(t) N^(t+r) N^lt+r)) 

s' 


<R^^jj^^(t,t+T)> = <m^(,t) m2(t) ffi^(t+T) m2(t+T)> 


(24) 


Since, from the definitions in (2), the noise processes N (t),f!l (t) and 

0 s 

the signal processes m-j (t) ,ffi 2 (t) are piecewise constant over a T 2 -second 
i nterval , then 


where 


<R,x M (t,t+T)> 


<R|fl^f;j^(t,t+T)> 


^ R^~i^*TtTt+Ty^ 
X s 




'^rii2 ’ 

Ul 

< T 

0; 

U1 

> T 

a.2 2 _M) . 

m^ N ' T 2 ' 

Ul 

< T 

0 ; 

Ui 

i — 
A 


Ul 

^ T 

0 ; 

Ul 

> T 

2 2/1 UK . 

o~ a~ { 1 - -tr^) i 
mi m 2 T 2 

Ul 

< T 

0 ; 

Ul 

> T 


(25) 


% = <N(?(t)> = <Ns^(t)> 


2T. 


and 


" ’ 






(26a) 


(26b) 
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is evaluated from Table 1. Substituting (25) and (26) into (23) yields 


<Rj^ (t,t+T)> 


4 

0 ; 


^2 ^1 '^m 


2To 


^n\2 


m 


+ Pl Poa/ (I --^4 ; 


< T, 


r2"m.('' ■ T. ' ? - '2 


• 1 ) '2 
III > T, 


(27) 


Finally, integrating (27) over t between and +«■ [see (15)] gives the 
equivalent noise spectral density 


«e' ' 


2 ^0 ^ ^1 ^2 ^2 2 
Po + P, 0 - + ^ + — a~ 

2 1 m-j 2 T 2 Ng m-j 


(28) 


whereupon, from (13), the variance of the 2<p process becomes 


4NoB, 


°Z<p - (P^ + Pj) 


(Pl^Pj) 


(P2-P,a/)2 


I _1 L_ 


P + n 2 ^0 , ^ ^1 ^2^2 2 

^2 “^l '^m^ 2 T 2 Ng °m^ 


Alternately, 


.2 _4_ . 

""lip ■ P.S|_ ’ 


A 


'0 \ ’ 


(29) 

(30) 


where »?|^ is the "squaring loss" defined by 

(P 2 -P 1 V)' 


'1 


p + p 2 _^g_ ^*"1 ^2^2 2 

^2 “^l ^m^ 2 T 2 Ng '^m^ 


(Pl+Pg) 


In terms of the power ratio Yp-P-]/P 2 * (31) can be put into the final 
desired form 

n 242 

n — I 


(31) 


P 1 + Yp ^ ^P ^T2 '^ni 
1 +Yn + -^=^7^ — + 


U + Yp): 


P ^m| 2 R^2 


tty: 


(32) 


where, as in Part I, Rj 2 - (P-j + P2)T2/f^0 total power-to-noise 

ratio in the high data rate bandwidth. Note that, except for the depen- 
dence of a- on the data rate ratio R^/R^ (see Table 1), the squaring 


loss of (32) is, unlike the results in Part I for passive arm filters, 
independent of this ratio. 

Figure 2 illustrates versus yp with fixed at values typical 
of coded and uncoded system's, arid as determined from Table 1 corre- 
sponding to the specific data formats of m-j(t) and m 2 (t). Note that, 
while the bi-phase Costas loop is unable to track a conventional balanced 
QPSK signal since, regardless of the value of 5^ = 0 when n = Vp = 1 and 
m^(t) and m 2 (t) are of the same format , it is indeed able to track a 
balanced QPSK-type signal in which m^(t) and m 2 (t) have identical data 
rates and powers but are not the same format. For example, if m, (t) is 
NRZ and m 2 (t) is Manchester coded, or vice versa, then when n = 1, = 0, 

and furthermore, when Yp = 1 , we have from (32) that 





T2 


(33) 


TRACKING PERFORMANCE -- THE CASE OF UNSYNCHRONIZED SYMBOL CLOCKS 

More often than not, the two data clocks will be unsynchronized 
since, typically, the two data streams are generated from totally inde- 
pendent sources. In this situation, there is, in general, no relation 
between the epochs of the synchronization pulses corresponding to the two 
data clocks and, furthermore, the ratio of the two data rates need not be 
an integer. Nevertheless, still assuming that the active arm filters are 
matched to the pulse shape- P 2 (t) , then <iri 2 (t)> is again given by (18) . 
Computation of <ifii^(t)> for the unsynchronized clock case is a bit more 
tedious than was the case for synchronized clocks. The nature of the 
mathematical model required is very much like that used in computing the 
cross-modulation distortion 0^2 fact, the results for 

<ifii^(t)> corresponding to the various data format combinations for m-j(t) 
and m 2 (t) can be obtained directly by taking the limit of ®i^^2 

approaches infinity [see equations (42), (50), (51), (52), for example]. 
These results are summarized in Table 2. 

Thus, to compute the loop squaring loss for the unsynchronized 
clocks case, we use (32) with now determined from Table 2. Note that, 
in the limit as Yj-+«, Table 1 and Table 2 are equivalent; i.e., for a 
sufficiently large ratio of data rates in the two channels,- the loop 


11 


Also: nii(t) is Manchester 


m 2 (t) is Manchester 


I n 


or: (t) is NRZ 


iTi 2 (t) is Manchester 


all n 


Also: 

m^(t) 

is Manchester 


ni2(t) 

is Manchester 

or; 

m^(t) 

m2(t) 

is NRZ) 

> all n 
is NRZ/ 

or: 

m^(t) 

is Manchester 


ni2(t) 

is NRZ 




Yp- P|/P2 


Figure 2a. Squaring Loss versus Ratio of Powers in the Two Data Channels 
With Yx = T]/T 2 = n as a Parameter; Rj2=-3 dB; m](t) is Manchester, m 2 (t) 
is NRZ; n odd; Synchronized Data Clocks 





•Also: nii(t) is Manchester 

m2(t) is Manchester 

op: is NRZ 

in^it) is Manchester 


Also: m^(t) is Manchester ^ _ A 

m2(t) i's Manchester I \ 

or: m, (t) is NRZ) _ 
m2(t) is NRZi 

or: iTii(t) is Manchester') 

' , , n even 

mo t is NRZ / 


Figure 2 b. Squaring Loss versus Ratio of Powers in the Two Data Channels 
With Yj = T-]/T2 = n as a Parameter; Rt 2=""'0 clB; m-](t) is Manchester, m2(t) 
is NRZ; n odd; Synchronized Data Clocks 




tracking performance is approximately the same for synchronized and unsyn- 
chronized clocks, as one would intuitively expectv^^ Figures 3 through 6 
illustrate versus Yp“P-j/P 2 ’ Rj 2 = PjT 2 /Nq fixed at the same values 
as in figure 2. The curves are distinguished from one another by the 
formats of the modulations m-|(t) and Note that, when the higher 

data rate modulation m 2 (t) is fianchester coded, the squaring loss is sig- 
nificantly improved relative to when this same modulation is NRZ. Again, 
the reason for this stems from the fact that the active arm filters are 

matched to the pulse shape of m,(t) which, for a Manchester code, results 

2 

in much smaller values of ajy,^ (see Table 2) than for an NRZ pulse shape. 

It is important to note that, here again, as for the synchronized clock 

case, the bi-phase Costas loop is quite capable of tracking a balanced 
QPSK signal (yj = Yp = i ) when m 2 (t) is Manchester coded., . 


Table 2. Evaluation of <mT (t)> ; ^ 





Manchester 

NRZ 

1 

- 1^+12 -6y^] i 

^T'6’'T 

Manchester 

< 2-R^ 

1 

< Fl^ < 

• 

2V "2 

1 - i ; fl2 >. 2fl, 

NRZ 

1 

■6v^ 

1.1 

3yt 


COMPARISON OF SQUARING LOSS PERFORMANCE - ACTIVE VERSUS PASSIVE ARM FILTERS 

We are now in a position to answer the question raised in the first 
paragraph of the introduction, namely, is there a tracking performance 
improvement obtained by employing active integrate-and-dump arm filters 
as opposed to passive arm filters and, if so, how much. To begin this 
investigation, we borrow a result from Part I, namely, (28), which charac- 
terizes the loop squaring loss for passive arm filters. In terms of the 













Figure 4a. Squaring Loss versus Ratio of Powers in the Two Data Channels 
With Yj = T-]/T 2 as a Parameter; R-j-p = -3 dB; m-](t) is NRZ, m 2 (t) is Man- 
chester; Unsynchronized Data CIocks 
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notation developed here in Part II, (28) can be rewritten as 

{^2 - YpD,) 


( 1 +v 


“2 Kd^ . Vp . K, (-^) {^) . (i;^) 


( 34 ) 


where , D^. Kd-j , and D-j^ are all defined in Part I. Note the simi- 
larity in form of (34) and (32), corresponding to the squaring loss for 
integrate-and-dump arm filters. 

Since , D 2 , Kd-| , Kd 2 » and D 12 /T 2 of (34) are only functions of 
B ^.//?2 and Yj> then analogous to the biphase case [1], for given values of 
Yj, Yp» and R.j- 2 , an optimum value of B ^//?2 exists in the sense of maxi- 
mizing of (34) (i.e., minimizing the squaring loss). Using that value 
of B^/i? 2 , namely, ('^-j/^2^opt’ defining the corresponding value of 5^ 
by then the minimum improvement in tracking performance (or, equiva- 

lently, in squaring loss performance) is given by the ratio of (32) to (34): 


T k _fk 


“opt 


' - > “m 


2\Z 


Oj.Kp^+Vp 0^ Kp^ +K|_ 


l^Ypw^^ 


^2 ^opt 


2 Yp Rj 2\/Pi2 

1+Y, 


(D2“ YpD^)' 


2-1 


2 1 + Yp ^ ’^T2 ^m. 

I + Y ■ J- 4 . L 4- d 

°m, ^ 2R^2 1 + Yp 


(35) 


Assuming single-pole (RC) arm filters as the basis of comparison, wherein 
K|_=l/2 [1], and unsynchronized symbol clocks (i.e., , determined from 

Table 2), Figures 7 through 9 illustrate I in dB versus the channel power 
ratio Yp with the data rate ratio Yj as a parameter and values of total 
power-to-'noise ratio R.p 2 typical of coded and uncoded systems. We observe 
from these figures that the improvement in squaring loss performance of 
using integrate-and-dump arm filters as opposed to single-pole arm filters 
is an*increa$ing function of Yp and depends heavily on the choice of 
data formats for m-|(t) and m 2 (t). Also note that, in the limit, as Yj 
approaches infinity, D-| = Kq.| = D^2/^2 " ^ ^ifi] “ ^ or 0, depending on 
the data formats of m.|(t) and m 2 U)- Thus, I of (35) simplifies consider- 
ably for this special case. 
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Figure 7. Squaring Loss Improvement (in dB) Using Integrate-and-Dump Arm 
Filters as Opposed to Single-Pole Passive Arm Filters; m.|(t) 

^ and ni 2 (t) are NRZ 
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Figure s. Squaring Loss Improvement (In dB) Using integrate-and-Dump Arm 
Filters as Opposed to Single-Pole Passive Arm Filters; m-|(t) is^ 
NRZ, MgCt) is Manchester 
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CONCLUSIONS 

^ We have evaluated the performance of a biphase Costas loop with 

active arm filters when tracking an unbalanced QPSK signal. The cases 
of synchronized and unsynchronized clocks have been treated and results 
are given for all possible combinations of NRZ and Manchester formats, 
for the two information channels. It has been shown that, while for a 
fixed ratio of data rates and total power- to-noise ratio in the higher 
data bandwidth the squaring loss does indeed increase as the ratio of 
powers in the two channels increases, the rate at which this loss increases 
(tracking performance deteriorates) depends heavily upon the data formats 
in each channel. Thus, it is concluded that, when the ratio of data rates 
is on the same order of magnitude as the inverse of the power ratio, 
i.e., approximately equal signal energies in the two channels, the 
biphase Costas loop can be used as an efficient demodulator of QPSK. 

On the other hand, if the energy in the two channels is very unbalanced, 
e.g., one channel is coded and one is uncoded, then it is still possible 
to efficiently use a biphase Costas loop for demodulation of unbalanced 
QPSK provided that the higher data rate channel is Manchester coded. It 
is understood that the foregoing conclusions are quite general and are 
not intended to rule out specific design situations in which sufficient 
total power-to-noi se ratio is available to tolerate large squaring losses. 
In an individual situation, one must resort to the specific numerical 
results given in the illustrations to determine the suitability of employ- 
ing a biphase Costas loop for demodulation of unbalanced QPSK. 

As a final note, the results derived in this part will play an 
important role in assessing the error probability performance of unbal- 
anced QPSK receivers where the data detector is, in fact, an integrate- 
and-dump fi 1 ter . The results of such an investigation are- the subject 
of a future paper by the author. 
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TRACKING PERFORMANCE OF UNBALANCED QPSK DEMODULATORS 

PART III - BIPHASE POLARITY-TYPE COSTAS LOOP 
WITH PASSIVE ARM FILTERS 

by 

Marvin K, Simon 
Waddah K. Alem 

1.0 INTRODUCTION 

In Parts I and II of this study, we addressed the problem of 
tracking an unbalanced QPSK signal with a conventional biphase Costas 
loop with analog input phase detectors, and an analog third multiplier 
(the one that forms the loop error signal). Because of de offsets 
associated with analog multipliers, it is common practice to hard-limit 
the in-phase* channel arm filter output and replace the analog third 
multiplier with a chopper-type device (switched multiplier) which 
typically exhibits much less offset. While it is also possible to 
replace the input in-phase and quadrature analog phase detectors 
with switched multipliers, the impact of doing so on the resulting 
tracking performance is minimal since the arm filters will pass only 

the first harmonic of these phase detector outputs. Thus, aside from 

2 

the 8/ir power loss associated with the first harmonic of a square 
wave, the performance of the loop would be identical to that given 
in Part I for an analog third multiplier or that to be presented here 
in Part III for a switched third multiplier. For ease of terminology, 
we shall refer to a conventional biphase Costas loop with a switched 
third multiplier as a "biphase polarity-type Costas loop" or, even 
simpler, a "polarity-type Costas loop." 

Generally speaking, introduction of a limiter (hard or soft) 
into a system results in signal suppression, the amount of which is 
a function of the signal-to-noise ratio,at the limiter input. This 
signal suppression, in turn, reduces the total loop gain and, as a 
consequence, the loop bandwidth. Another potential problem with the 
limiter under strong signal conditions is that it may increase the 
tendency of the loop to false lock. This tendency was demonstrated 


* 

For unbalanced quadriphase, we shall refer to the in-phase 
channel as that corresponding to the point of data extraction for 
the higher power signal. 
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in [],2] for a polarity-type Costas loop with a biphase data input. 

The purpose of this appendix is to address the tracking behavior 
of the polari ty-type Costas loop with unbalanced QPSK input and com- 
pare its performance with that of the conventional Costas loop dis- 
cussed in Part I of this same study. In particular, for the case of 
single-pole Butterworth (RC) arm filters and a particular combination 
of NRZ and Manchester coded data on the two channels, the squaring 
loss (tracking jitter penalty relative to a linear loop) is evaluated 
and illustrated as a function of the ratio of arm filter bandwidth to 
higher data rate and total signal power- to-noise ratio in this higher 
data rate bandwidth. Also numerically illustrated is the correspond- 
ing mean-squared tracking jitter performance as a function of these 
same receiver parameters. 

2.0 SYSTEM MODEL 

Consider the Costas loop with hard-limited in-phase channel 
illustrated in Figure 1. The’ input x(t) is, as in Parts I and II, an 
unbalanced QPSK signal plus noise, i.e., 

x(t) = /2'pYni-j(t) cos <t(t) + m 2 (t) sin <i'(t) + n^. (t) , (1) 

where 4>(t), mi(t), m 2 (t) and n^. (t) have been previously defined. To 
avoid unnecessary duplication, we shall not rederive the expressions 
for the in-phase and quadrature arm filter outputs, z (t) and z_{t), 
but rather refer the reader directly to the results in equation (5) of 
Part I. Thus, hard-limiting z^(t) and multiplying by z^{t) produces 
an error signal ZQ(t) at the switching multiplier output (assuming 
this multiplier has unit gain) which is given by 

Zg(t) = z^(t) sgn [Z 5 (t)] 

= K-| m 2 (t)m(t) sin (p{t) + m^(t) m(t) cos ^(t) 

+ K^m(t) [Nc(t) cos f/?(t) - Ng(t) sin <;o(t)] , (2) 

m(t) = sgn m 2 (t) cos ^(t) - t^m^(t) sin cp(t) 

- [Ng(t) cos (p(t) + Aj,(t) sin (/?(t)]| (3) 


where 


ro m 2 (t) Data Detector 


In-Phase Arm ("I" Channel) 


e (t) Low-Pass 

►-i ^ Filter 


r,(t) 




-1 

+1 


Acquisition Voltage 


r (t) 
c 



Loop 

J-U 

Filter 

• 

F(s) 


e^(t) Low- Pass 
? ^ Filter 


Quadrature-Phase Arm 
("Q" Channel) 


To m-j(t) Data Detector 


Figure 1, Costas Loop with Hard-Limited In-Phase Channel 
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and sgn (x) = x/jx] denotes the "signuni" function of the argument x. 

The instantaneous frequency of the VCO output is related to ZQ(t) by 

= K^,[F(p)zQ(t)] + (4) 


and hence the stochastic differential equation of operation for the 
loop in Figure 1 becomes 

2 = 2f?Q - K F(p) ^2 m 2 (t) m(t) sin (/?(t) + 2 /Fj~ m-i (t) m(t) cos (p{t) 

+ 2m(t) N [t,(j9{t)]| (5) 

v/here K= Ki K^, and 

N[t,(^(t)] = N^(t) cos (p{t) - N^(t) sin (p{t) . (6) 

Consider first decomposing m 2 {t)m(t) sin (p{t) into 

m 2 (t) m(t) sin 9 ?{t) = <m 2 (t) m(t) sin <;p(t)> 

+ ni 2 (t) m(t) sin (jc)(t) - <m 2 (t) m(t) sin 9 ?(t)^ (7) 

where the overbar denotes statistical expectation with respect to 
both the noise and data modulation processes, and < > denotes time 
average. Similarly, 


m^ (t) m(t) cos <;o(t) = <m^ (t) m(t) cos (j?(t)> 


+ jjm-j (t) m(t) cos (pit) - <iTii (t) m(t) cos <^(t)> 


(8) 


The first terms of (7) and (8) represent the normalized dc components 
of the error signal at the switching multiplier output. Since these 
dc voltages are, in general, a function of the loop phase error 
(specifically, in this case, a function of twice the loop phase 
error), it is traditional to refer to their weighted sum as the loop 
S-curve or tracking characteristic . With this in mind, we can write* 


* 

For convenience, we drop in our notation herein the dependence 
of (p on time. Also, the loop bandwidth is assumed to be much less than 
the data rate and thus (V5(t) is essentially constant over a symbol time. 
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!*• 






2 m(t) sin </ 9 (t)> + 2 <irii (t) m{t) cos (/ 7 (t)> = 

2 y^<ni 2 (t) m{t)> sin (/? •<• 2 <rn-| (t) ni(t)> cos 9 ? - 

= (9) 

where f-j(x) is a nonlinearity which is periodic in x with period 2 tt 
and has unit slope at the origin, i.e., x = 0. Thus, a defined by 




0^2 ^^1 

is the signal amplitude suppression factor which results due to the 
combined distortion effects on the input modulations m|(t) and ni 2 (t) 
of the hard limiter and finite arm filter bandwidth. In the next sec- 
tion, we evaluate this suppression factor for modulation data formats 
and an arm filter type of current interest in practice. 

The second terms in (7) and ( 8 ) are zero mean processes which 
represent the -self noises of the filtered, hard-limited data modulations 
at the switching multiplier output. Thus, defining these self noises 
by 

n^^ [t, 2 c/?] = 2 /P^[m^(t)m(t) - <m.|(t)m(t)>] cos 


= 2 /P^[ni 2 (t) m(t) - <ni 2 {t) m(t)>] sin 9 ? ( 11 ) 

and substituting (11) together with (9) into (5) gives the resulting 
equation of loop operation 

2 ^ = 2no-K’F(p) |/Pj a (2?>) + 2 i 5 (t) +n^2^t,Z<p) + . 

( 12 ) 


Finally, since for good tracking performance the loop operates in the 
region of small variations of (p (i.e., its linear region), then 
since is proportional to sin (p> its self-noise power in 


1 


6 


the narrow loop bandwidth will be vanishingly small when compared with 
that of the dominant noise component m{t)N{t,(p). Thus, as has been 
done in similar analyses [3,4], we shall ignore the effect of this 
self-noise term. Under this assumption (12) reduces to 

where 

Ng(t,2(/p) = 2m(t) N(t,<v5) + n^^(t,2(^) (14) 

is the equivalent additive noise perturbing the loop. 

3.0 CALCULATION OF THE SIGNAL AMPLITUDE SUPPRESSION FACTOR 

The input digital modulations m^ (t) ; i = l,2; may be represented 
as 

CO 

" I ^n2 ^2^'^" '''’^ 2 ) ’ 

n=-« 

00 

iii-](t) ^ -I a^^^ p (t- nT^ - e) (15) 

n=-“ 


where a^^ = +1 is the polarity (±1) of the ]ith data symbol in the i_th 
data stream, T^= 1/i?^. is the symbol time and p^-(t) is the symbol pulse 
shape of the j_th data stream (e,g,, for NRZ coding) 


P^(t) 


1 ; 0 ^ t < T. 

0 ; otherwise 


(16) 


and e is the random epoch between the two data streams which accounts 
for the fact that, in general, the two modulations are asynchronous. 
After passing through the Costas loop arm filter, tnese same 


modulations become 



m.it) = 

00 

1 a^. p.(t- nT.) ; i = l,2 

(17) 

where 



P|(t) - 

rT 

P,-(t:) g(t - t) dx • 

" SCWSh 

(18) 


with g(t) the impulse response of the arm filter G(s). Thus, from (10) 


= <m2(t)m(t)> 
which, from (3) and (17), becomes 


1 

T 

r^2 

m^(t) m(t) dt 

9?=0 '2 J 

0 


<p=0 


( 19 ) 


2 “ T, 


T, 


I a^gPgit-nT,) 


n=-co 


sgn I a^2 P2^^“ ^^2^ 


n=-c 


( 20 ) 


The statistical average over the filtered Gaussian noise process N^(t) 
in (20) can be performed quite easily with the result 

/p£ 


S3" j'^2 - "s<t> 


= erf 


2 ^ ®n2 ^2^^ " '^"^2^ 

n=-" 


( 21 ) 


where p 2 = 2 P 2 /NqB^. is the high power channel signal-to-noise ratio in the 
arm filter bandwidth, 


. A 2 
erf X = 

/tT 


p 

exp ("t") dt , 


( 22 ) 


and denotes the two-sided noise bandwidth of the low-pass arm filter 
G(s), i.e., 




|G(o 2nf)r df , 


(23) 


Substituting (21) Into (20), we get the somewhat simplified result 




I ^n2 ^2^^' '^^2^ 


n=-oo 


erf 


^ %2P2<*‘"'^2> 


n = -ra 


dt, (24) 


In general, evaluation of as given by (24) is a difficult, if 
not impossible, task due to the presence of the error function non- 
1 inearity which impedes performing the statistical average over the 
doubly infinite data sequence. Fortunately, however, the problem can 
be somewhat simplified by noting that, in practical system design, the 
3 dB cutoff frequency of the arm filter is chosen on the order of the 


data rate or greater to assure good tracking performance. This implies 
that the tails of P 2 (t) beyond a single synibol interval decay rapidly, 
and thus it is sufficient to consider only the contribution of the 
data b.its 3 q and a_-j insofar as computing the argument of the error 
function in the interval 0^tiT2** Assuming then that the data bits 
are random and take on values ±1 with equal a priori probability, we 
obtain the result 


1 


2T. 


where 




P 2 ^.(t) erf 



dt + 


- p2(t-i-T2) 


2T, 


fT, 


p'o ( t) erf 



0 ^ t i T 


2 ' 


2 P2-(^^ 


(25) 


dt 


(26) 


1 

Computation of a-j defined in (10) is slightly more complicated 
since the time origin of both data streams is involved and thus an 
additional average over the random epoch e between them is required. 
Averaging first over the filtered Gaussian noise process N[t,<^(t)]- 
N^(t) cos i;(9(t) + N^(t) sin (p{t) with variance = N^B ./2 gives 


_d_ 

d(p 


m, (t) erf 

V I 


/pj ”12 ^ ^ ^ ^ ~ *^1 ^ 

) sin (p 

\ 

“n 


/ 


(p-Q 


'^Pl /. 2,.. ) P2 . 2,, V 

m-] (t) exp j- 


(27) 


where 


1 


2P 


N.B. 


(28) 


Again, because of the independence of the two data streams and 
the asynchronism of their timing, (27) can be simplified to 


Extension of this approach to take into account intersymbol 
interference from more than one previous pulse interval follows in a 
straightforward manner. 



9 


~~ <(m/(t))>/exp j- -y in 2 (t)j 


The first factor in (29) is easily evaluated as [1] 

^ = f S^^(f) |G(j2,if)j^ df , 

Loo 1 

where Sm^(f) is -the power spectral density of the modulation m-|(t). 
The second factor in (29) can be simplified by once again using the 
approximation of considering only intersymbol interference from the 
previous pulse as was done in reducing (24) to (26). Thus, 


1 f^2 


^^2 Jo 


^2 - 2 

yP2+< 


1 f^2 


TT '^m, 2T 


1 "'2 JQ 


2 a2 
2 ^2- 


■_(t) 


Finally, letting 


A ^ 


denote the ratio of powers in the two data channels, then from (10), 
(26), (31) and (32), the signal amplitude suppression factor a is given 
by 


1 

ypl-\w;\, -p 


1 (^Z 


’'2 A 2 ^ 4 - 

2 P2+^^^ 


^2 ji 2 / _ V J j. 

2 P2- ^ ^ dt 


Note from (33) that the effect of the quadrature data channel (Channel 1) 
on the signal amplitude suppression factor is to reduce its value rela- 
tive to that [the first two terms of (33)] when the input is a biphase 


10 


m 


m 


9 


9 


9 


9 


m 




modulation [4]. 

As an illustration of the application of the result in (33) to 
a particular case of practical interest, we consider a Costas loop with 
single-pole (RC) Butterworth arm filters and input modulations specified 
as follows. The high data rate modulation m 2 (t) is characterized by 
NRZ data with pulse shape p 2 (t) specified in (27), while the low data 
rate channel has a Manchester coded modulation m-j(t). 

As is well known, an RC filter has a transfer function 

G(ji») = (34) 

where to = 2irf is the radian 3 dB cutoff frequency. Taking the inverse 
Fourier transform of (34) and substituting this result, together with 
(16) i into (18) results in 

p 2 (t) = [1 - exp (-iij^t)] u(t) - [1 - exp (-o)^(t - T 2 )) ] u(t - T 2 ) (35) 

where u(t) is the unit step function. Thus, from (26), (33), and (35), 
we obtain - 



where we have made the substitution x-t/T 2 and taken note of the fact 
that, for an RC filter, the two-sided noise bandwidth and 3 dB 
radian cutoff frequency co are related by w -2B.. Also, the mean- 
squared value for a Manchester coded modulation after arm filtering 
is given by [5] 

\ ^ ^ exp (-B./i?^) + exp (-28^:/;?^)] . (37) 


Before illustrating the variation of a with we must char- 
acterize the signal-to~noise ratio P 2 in terms of the system parameters. 
From the definition of p^, we have that 


i!L 


'P T 




- R. ^ ' 
ir 2 




or in terms of the ratio of arm filter noise bandwidth to data rate. 


2R2 

^2 " bTT^ * 

Since, in a given system design, the data detection signal-to-noise 
ratio R2=P2’r2/^0 typically fixed by the requirement on bit error 
probability performance in the high data rate channel, we shall param- 

9 

eterize R 2 and plot the signal power suppression factor a as a function 

of B^./i? 2 - Figure 2 illustrates such a plot. The curves are drawn for 

By/fp il.5 or, equivalently, t /T ^ 1/3, where t = 1/w is the arm 
1 ^ c c c 

filter time constant. 

9 

For large B^./i? 2 s an appropriate asymptotic expression for a 
can be obtained from (36) and is given by 

T / «2 \ 

(- Byflj) • 

Also, in the limit as R 2 ”, we obtain the asymptotic expression 

= 1 - 2 fj^nn 2 - exp (- 2 B^./i? 2 )] , 

which is identical to the result obtained for a biphase modulation 
input [4]. 

4.0 CALCULATION OF THE EQUIVALENT NOISE SPECTRAL DENSITY 

In practical applications, the bandwidth of the Costas loop is 
ordinarily designed to be narrow with respect to the equivalent noise 
bandwidth of Ng(t,2^) defined in (14). Thus, we can further approxi- 
mate N^(t, 2 f/?) as a delta correTated process with equivalent single- 
sided noise spectral density 
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N, 2 


,co 


R|^ (t) dx 
e 


(42) 


where 


Rfi (T) = <N^(t,2(^.)N^(t+T.2¥))> . 


(M) 


Since, as previously mentioned, for good tracking performance, the loop 
will operate in the neighborhood of (f-O, vie shall evaluate (43) and 
hence (42) only at the point cp = 0. For a conventional Costas loop 
without the hard limiter in the in-phase arm, the equivalent noise 
spectral density is indeed independent of cp [3]. Proceeding then under 
the foregoing assumption, we have from (3), (6), (11), (14) and (43) 
that^ 


"" 4<^N^(t) N^(t+x) sgn |/pj m2(t) - Ng(t)| sgn |/P^ jri2(t+x) - N^(t+x)J^ 

+ 4 <(m^ (t) m-j (t+x) sgn Fi2(t) - N5(t)| sgn |/pj m2(t+x) - N3(t+x)|^ 

= 4 [Rj^jCx) + R-^(x)3 <^sgn m2(t) - Ng(t)| sgn |/^ m2(t+x) - N^(t+x)|^ 

(44) 


where 


and 


Rr(t) = Njt)N (t+x) = |G(j 27 rf)l^ df 


R. (x) ^ <m^(t)m^(t+x)> = S^(f) !G(j‘2irf)|^eJ^’''^'' df 

1 ■'-CO 1 

For the RC filter of (34), it is simple to show that 

RfiU) = 


(45) 


(46) 


(47) 


Furthermore, for Manchester coded data on Channel 1 with power spectrum 


Note that, because of the independence of m^(t) and 1112(6), the 


term <mi (t) m{t)> = <m^ (t) sgn m2(t) - in n^^[t ,2( 

evaluates to zero. 
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. 4 

sin ' 


V(f) = T, 

1 / ir f Tt 


( 48 ) 


{^) 


we have that 


1 - 3 JpJ' + 

h 


4e - e 


- 2 B 4 ItI 

cosh 2B^ T “ 3e 


2B./i?-j 


; 0 ^ T I ^ 


-281x1 ^ -2Bi/i?-| 

/ i^i\ e '[4cosh B./^i - 3] - e ^ ,'cosh2B,x T, 

S'" ■ V<”^) 


-2Bi I T I 

e [4 cosh - cosh 2B ./^’-j - 3] 


; T-j .< T 


(49) 


Thus, what remains is to compute the correlation function of the hard 
limiter output as required in (44). 

The general problem of computing the correlation function at the 
output of a zero memory nonlinearity (e.g., a hard limiter) when the 
input is signal plus Gaussian noise can be attacked by a variety of 
techniques, among which are the transform approach [6], Price's Theorem [7] 
and the series expansion method [8], Letting 


y{t) = sgn {/P^ m 2 (t) - N^(t)} , 

then to compute the equivalent noise spectral density defined in (42), 
we must calculate' 


(50) 




= 8 


[R^(t) + P^ Rfn(t)]Ry(x)dT . 


(51) 


m 


As was true for the signal suppression factor and even more so 
here, evaluation of R^(x) when mp(t) is specified by (15) is difficult 
if not impossible. Even when one is willing to make the simplifying 
assumption of requiring that the statistical average on the data be 
taken only over two data symbols, additional simpTifying assumptions 
are necessary if one is to obtain meaningful numerical results for the 
equivalent noise spectral density. r'C' IS 
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It is convenient to start with the power series approach of [8] 
to evaluate R (t), namely, 

V-o = 




where H|^(x) is the j<th degree Herniite polynomial defined by 



HJx) = -k2/2 

dx ■ 

(53) 

and 

4 «n(") _ -2Bi|t| 

(54) 


The case of interest here is where the data on Channel 2 is uncoded, in 
which case, the signal -to-noise ratio would be large to achieve a 
small error probability on this channel. Thus, since is directly 
proportional to for a fixed ratio of arm filter bandwidth to data . 
rate then for large p^, the first term (k=l) in the series expan- 
sion of (52) will dominate. Thus, keeping only this term and again 
assuming that the statistical average on the data sequence is computed 
only over two symbol intervals, we obtain 
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"‘^'2 


erf 


ni2(t) 


erf 


Z ^ / V 

■y ni2{t+x) 


dt 


dr 


+ 4N„B.(f) 


rZT ^ - 4 B. |t| 
e 

-2To 


T. 


exp -^[ni 2 ^(t) + ni 2 ^(t+t)] I dt 


2T, 




1 f^2 


2 JQ 


+ 4N„B. 0 , . 

0 T 1 j_2t^ 

,? -2B-1 1 t| 



"/pp 


/pp 

erf 

Jt '^'2^^^ 

erf 

j—^ m2( t+x) 


dr 


dt 


dt 


1 r'^2 


2 ■'0 


exp -^[m 2 ^(t) + ni 2 ^(t+T)]| dt 

(55) 


where m 2 (t) and rfi 2 (t+T); 0^t^l2j 0ji|t|:<2T2; are computed from (17) 
and (35) and R|n-|(T) is given by (49). Carrying out the statistical 
average, v/e obtain after much manipulation 


Ne - 


where, analogous to Equation (57) of [4], 3^, is defined by 




-2Bz/R2 

f](y) + e 


dy 




+ (2B^.//?2) 


fT -(4Bi/i?2)y 

e 

0 


“4B.j/i?2 


dy 


r1 






+ (2B./i?2) 


2p-, |-1 ..(2B.j/i?2)y 
0 


(56) 


R. 


dy 

2Bi/-!? 2 /y+lN 




fa(y) 


dy 


with 


( 57 ) 


I 


fi(y) 


p ,1 p /pp ry 

= erf J-y- - g^(x) dx + g^(x) dx 


1 


1-y 


+ ^ g^(x) g^ix+y) dx - g^(x) g-,(x+l-y) dx 
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D 




fgly) 


f3(y) 


(i^) exp^ (“P 2 / 2 ) + 92(x) dx + ^ g 2 ^x) 


1 f'i-y 1 fV 

? Iq ^ 2 ^^^ + 4 92^^^ 92 (x+^-y) ^x 


ry 

0 


1 V ? r^? r 

(1^) ~ J] 


1 

4 


(y 


g-| (x) g|(x+y) dx 


© 


^4(y) 


1 2 If 1 f-^ 

6 xp (-P2/2) + ~2 g2(x) dx + 92 (x ) g2(x“^1 ~y ) dx 




1-y 


g 2 (x) g 2 .(x+y) dx 


(58) 


© 


© 


and 


g^ (x) = erf 


g.(x) = exp 


^ (1 - 


1 - 2e 


-(2B^./ff2)^~ 


(59) 




■© 


# 


5.0 CALCULATION OF SQUARING LOSS AND TRACKING JITTER 

We have already observed in previous analyses of this type [3-5] 
that the performance of Costas loops in either their linear or nonlinear 
tracking regions is characterized by the behavior of a parameter referred 
to as squaring loss . This loss represents the additional degradation in 
equivalent loop signal -to~noise ratio relative to a linear loop, i.e,, 
one in which there exists no noisex noise or signal x noise error signal 
components, More specifically, if in (9) the normalized S-curve f‘|(2c^) 

page 




is replaced by Zip, its linear approximation, then one obtains the 
linear equation of loop operation from which the effective loop signal - 
to noise ratio p' insofar as computing the behavior of the Zip process 
is 

a)^ 




( 60 ) 


where is the single-sided loop noise bandwidth and represents 
the noise power of the equivalent noise process N^(t,2^) in this 
bandwidth. 

Using (56), we can express (60) in the form 

'■ • iCw)(^)©- 

and, in addition, for fiQ=0, the tracking jitter is given by 


where 


p = 


^Zip 

Pi + P, 

“¥1: 




and 


(61) 


(62) 


■ (63) 


^ + Yp / 

is the so-called squaring loss with a given by (33), given by (57), 
and Yp defined in (32). 

Since the modulation reference signals are at odq rather than 
2 uq, then the tracking jitter on the output data is 

2 1 2 1 • 


(64) 




1 2 
4 ""Zip 


pS^ 


Using (38) along with the definitions of a, 3 ^^ and Figure 3 
illustrates the variations of versus B./i ?2 with PjT 2 /Nq as a parameter. 

Superimposed on these curves (in dashed lines) are the corresponding 
results obtained from [5] for the biphase Costas loop with passive arm 
filters. Figure 4 shows the phase jitter aqj as a function of B^./i ?2 
for the hard-limited loop with P.J.T 2 /NQ and i? 2 /ff-j as parameters. The 
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Figure 3, Squaring Loss Variations Versus vvith 

as Parameters r in.j (t) is Manchester code, m 2 (t) is NRZ. 







phase jitter of the conventional loop [C] for R^/H-^=Z0 is also included 
for comparison. The computations for the phase jitter are based on a 

4 

typical ratio of B^/Bj^=10 . We observe from these results that, for 
high signal-to-noise ratios, the hard-limited loop actually outperforms 
the conventional loop. 

6.0 CONCLUSION 

The squaring loss of biphase polarity-type Costas loops (with 
a switched third multiplier) has been analyzed in this report. An 
improvement in the tracking performance of these loops at high signal- 
to-noise ratios over the biphase Costas loop with passive arm filters 
has been demonstrated. The numerical .evaluation shows that the improve- 
ment (in terms of the squaring loss) might be as high as 2.8 dB. It 
has also been shown for a given signal-to-noise ratio and a given arm 
filter bandwidth to high data' rate ratio that the squaring loss does 
not change significantly with data rates when the ratio of the data 
rates is high. This is particularly true for small values of arm filter 
bandwidth to high data rate ratio. 
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SUBCARRIER TRACKING ANALYSIS FOR THREE-CHANNEL ORBITER 
KU-BAND RETURN LINK 

by 

Marvin K. Simon 
Waddah K. Alem 

INTRODUCTION 


In a previous report .(Ref . 1, Appendix C) , we studied the perform- 
ance of a Costas loop for recovering the carrier from the three-channel 
quadrature multiplex signal transmitted on the Orbiter Ku-band return 
link. We observed there that the output of the quadrature phase detector 
of this loop serves as the input to the subcarrier tracking loop (see 
Figure C-2 of Ref. 1, Appendix C). Upon establishing a subcarrier refer- 
ence signal, the two lower rate modulations, ni 2 (t) and m 2 (t), can then 
be demodulated. This appendix presents the tracking performance of a 
Costas loop used for subcarrier recovery (Figure 1). 


SYSTEM MODEL AND LOOP EQUATION OF OPERATION 


Based upon the above comments, the input signal is of the form 
’see equation (C-15) of Ref. 1, Appendix C): 


x(t) 


1^2^ S2(t) + S3(t) + N^(t) 


cos (^^(t) 


+ |/2P^ s^(t) - s^{t) s^lt) S 3 (t) - N^(t)j sin <^^(t) 


Sd(t) 


0 


where (p^{t) = carrier tracking loop phase error 
s-|(t) - m.|(t) - high rate modulation 

Sp(t) = mp(t) Sin w t = low rate modulation ^i^ on square-wave 

subcarrier 

S 3 (t) - ni 3 (t) Cos u) t = low rate modulation i¥2 on quadrature 

square-wave subcarrier. 

The square-v/ave subcarrier frequency f^^ = u)^^/2Tr is nominally 8.5 MHz, 
while the modulation m 3 (t) is operation data at a rate - 192 kbps and 
the modulation m 2 (t) has a rate up to 2 Mbps, Thus, assuming that the 
input bandpass filter H(s) has amplitude and phase characteri sties which 


( 2 ) 
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are symmetric around its center frequency f , and a bandwidth which is 

S 0 

relatively wide with respect to the two data rates, e.g,, on the order of 
the subcarrier frequency, then letting 


r^(t) = >^ cos [o)g^t - 

r^(t) = PZ sin [u^^t - 


(3) 


and ignoring second harmonics, we get for the quadrature phase detector 
output:* 

Cc(t) = K, cos <p^ [-f + ±/p^ *3(1) cos <p, 

+ /2 N^(t) cos (ai^^t - 


sc 


sc 


+ K, K„ sin <p^ 


^ "'lc(f>“sc'-‘Psc) "dc(‘'“sc’- 

W 


- >^ N^(t) cos - (P^^) 


In (4), K is the phase detector (multiplier) gain, cp is the subcarrier 

m SC 

phase error, and the remaining quantities are defined (in operator form) 
by 

ni2(t) ^ Hjp) iTi2(t) 
m3(t) = Hj^(p) m3(t) 

mic(t;a>^c’-'^sc) “ “^'sc)^ 


(5) 


and 

Njlt) ^H(p) Nj(t), 

N^(t) ^ H(p) N^(t) 

where H^(jo)) is the lowpass equivalent of H(jui) and is defined _by 




1 

2 


H[j(to + w^^)] + H[j(w *- Wg^)] 


( 6 ) 


(7) 


m 


on t . 


For convenience of notation, we omit the dependence of (p and (p^ 


sc 


4 


Expanding the filtered noise processes N (t) and N (t) as 

N^(t) = Jl [N^^(t) cos Wg^t - N^ 5 (t) sin u^^t] 

N^(t) = /2 CN 5 ^(t) cos ui^^t - Ngglt) sin m^^t] , (8) 

v/here N (t), N (t), N (t), and N (t) are approximately statistically 

CC .CS Sw ss 

independent, stationary, lowpass Gaussian noise processes with power spec- 
tral densities 

Sj^(iu) = y |H^(joj) I , (9) 

then substituting (8) into (4) yields 

E^(t) = K, cos <p^ iSjCt) sin iSjd) cos 

+ cos < 0 ^^ - N^^(t) sin 

+ K, sin cp^ /F^ 'idc<*’“sc'-‘^sc) 


,(t) COS + N^^(t)-sin . (10) 


Similarly, the output of the in-phase multiplier can be shown to be 

^s^*' ° 'Sn ‘f’c 4'‘^^2 1’sc ■ 'V''% ®sc 

- fi„(t) sin - N^^(t) cos 

V’c "is<*-"sc'-‘^sc> - ™d5‘‘-V»'sc> 


+ Ng^(t) sin ^sc ’ 


where (in operational form) 


™ls(*’“sc’-‘^sc> <"sc‘ -»’sc>] 

'"ds<ti“sc’-‘'’sc> - t'd(*) <“sc‘ "^sc>^ 


After lowpass filtering with in-phase and quadrature arm filters G(s), the 
phase detector output signals z^(t) and z^(t) are given by (10) and (IT), 
respectively , where the "hats" now denote filtering by the cascade 
Gg(s) - H^(s) 6(s). Multiplying these two low pass filter outputs 


5 


(assuming, for simplicity, that tin's multiplier has unit gain) gives the 
cfynaim'c error signal: 


= Im 2 


n ^ 


— j^cos - p3m3^(t)] sin 2 cp^^ 

32 2 / — — 

t^cos /p^Pjfnjtt) m3(t) cos 2g.^ 

"'^i"'»’c^^ls(tlV-?’sc)™lc(t->“sc--?’sc) 

" 2 ='"''^c^™ds(‘’V-^sc)”dc<‘^“sc'-?’sc) 

- 2 sin^^a^ ■^^"'is(t:V’-'?’sc) '"clc<‘=“sc'-'^sc’ 

+ — sin 2 <p^ /pj I'n^ttl^sin ^^(t) sin 

■ Ss(‘>“sc>-^sc> - ^ "ds<t>sc--^scr ■ 

„ 2 c), P*2^h) cos - /p^ cos (p^^ 

• "’lc(ti"sc’-'^sc) - "dc(‘-“sc'“'^sc)_ 


V 2 (t;(;o^, 2 ^^^) = j\^ sin 2<p^^ + A^ cos 2cp^ 




col t= coi OJ CM . CM col {= cof t= 


COS (p^ - NjT(,(t)] Sin 2cp^^ 


+ sin ^ cp^ 


+ 2 sin (p^ cos cp^ ^sc^^^ ‘ ^cs^^' ”ss^^^^ 

Ncg ( t ) tri 2 ( t ) + / p ^ 

Ns 3 ( t ) m ^ lt ) + / P ^ m ^( t ) 


2 

cos ip^ 


s^n (p cos <p 


.2 


sin w 


1 


mT ^( t ; o)^J ~ / P ^ ’^ dc ^^’"' sc ^ 


1 '" lc'''»‘"sc 




d "'ds'"’'“sc' 


ss 


(t) 


sin <p cos (p 


nilc ( t ; o ^ sc ^ - /h 


d "'dc'^’^^sc' 


Ss^^’“^sc^ " '^ds^^’‘^sc^ 


^cs(^) 




sin cp^ cos (p^ [Ncs(t) N 3 ^(t) + N^^(t) N^ 5 (t)] 
/^ Ncc(t) 


2 

cos 


sin (;c? cos (p 


~^2 ^ sc ^^) ^-( t ) m ,( t ) 


3 ’’ss' 


+ 2 sin ^ 


/F7 m, (t;cu^J - m. (t;t._) 


1 '"lc'''’^sc^ " ’ ’d '"dc'"’‘^sc 
/^ Ss «=^“ sc > - ■^ * dr .< t ’“ sc ) 


ss 


(t) 


sc 


(t) 


- 2 sin (p cos ip 
c c 


/P7 Slc ( t 5 “ cr .) - ™ H ,( t ; i»,J 

1 Ss<t;“sc 


1 "'1g''-’“sc' ■ ' "d "'dc''-’“sc 

- /Pj ’nt|3(tiWj|.) 


Ncc<‘> 
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m 


& 


and 

[m^(t) cos w^^t] 

iS, 5 (t;«.sc) = m, 5 (t;o.j^.O) - Gg(p) [in,(t) s1n 

™dc<‘'V> - ™dc<*’“sc’“> ° ®0<P> 

i«ds(t;“sc> - '*ds<*-“sc-°> ' G(,(p) [Sj(t) sin (17) 

Each of the signal terms in (13) can be broken up into its mean value and 
a self-modulation component around tht mean, e.g., 

PjiSj^it) - P3in3^(t) = <[P2iS2^(t) - p3iS3^{t)]> 

+ P2[m2^(t) - <iS2^(t)>] - P3[ra3^(t) - <rfi3^Ct)>] 

(18) 

As before in analyses of this type (Ref. 1, Appendix C; Ref. 2), we shall 
neglect the effects of the self- modulation noise components. Furthermore, 
for the cross -modul ation noise components, we may set (p^ ~ ~ which 

eliminates all but that due to the m^lt) (^^(t) product. Finally, by 
analogy with (21) and (33) of Ref, 1, Appendix J, 


<"'ls(tiV-»'sc) ®lc(‘’“sc’-‘^sc» ' 
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(» 


C 


« 


where = 1/T-j is the data symbol rate of nii(t). Since (2tiI<)/T-j >> 
for any k 0, and for k = 0, 


„„ = 0 , 

""sc 'o 


then 


Similarly, 


( 20 ) 


( 21 ) 


<"'ls(*-“sc’-f’sc) '"dc(*’“sc”‘Psc» ' ° 


( 22 ) 


Also, since m.(t) and m.(t), i j , are independent, we have, for example, 

1 iJ 


ni 2 (t) m 3 (t) = 0; ni 2 (t) = 0; etc, 

Thus, (13) simplifies to 


(23) 




2 2 


2 

(^) cos^ [Pg Dj - P 3 D 3 ] s1n 2<p^^ 


2 

+ 2(4") /(\ P 3 ^^23^^^ ^ 


where 


Di - <iSf(t)> = S„.(ffl) Gg(jB) ^ g- . i = 2,3 (25) 

•'-CO 

with Sin^-(w) denoting the power spectral density of m|(t) and n 23 (t) = 
m2(t) ni2(t) . 

The instantaneous frequency of the VCO output is related to ZQ(t) 
by ' 


2 \ 


2 du) 


(24) 


iiM = K^(F(p) Zg(t)] 


(26) 


where F(p) is the loop filter transfer function and Ky is the VCO gain in 
rad/sec/v. Thus the stochastic integro-differential equation of loop 
operation becomes 


- K F(p) j(-t/ cos 2 [P^ Dj - P3 D3] sin 

'^^2^3 " 23 ^*' ''2(’'’‘Pc’^''sc'j 

(27) 

where ^V’ 

CHARACTERIZATION OF THE EQUIVALENT ADDITIVE NOISE 


Vie now proceed to evaluate the noise spectral density, N^, of the 
equivalent delta-correlated noise process, n (t;^ ,2^ ): 


f^e = 2 ^(v) ^2^3 ^23^^) ^ 

' •'—CO CO t- 


(28) 


where 


Rv 2(^} = ^2( ^ ’‘^c ’ ^*^sc ^ ^^2^ ’^c ’ ^*^50^ 


R 23 ( "^ ) = <np^(t) n^^(t+T:)> 


23 V w '>23' 


(29) 


Substituting (14), (15), and (16) into (29), we obtain, after considerable 
algebraic manipulation 

,2 

rp- f. , 
m 2 


Rv^(-c) = Rm^(0 + 4(-J-) CP2Rrii^(T) P3Rfn3(^)3 %(r) 


+ 4 


Ri (Rm, (-Tiw,,) + R^;, (t;i 


sc 


iV'mic' ’ sc^ ' ''m-,s 

+ PdK'dc‘"’V> Ss*'’“sc>) 


f-aU) 


R23(t) = Rm.('r) Rm_('r) 




'm- 


(30) 


where 


lU 



A 

<ni|(t) m^(t+T)> ; i = 2,3 

(■'■•“sc) 



("•“sc’ 

A 

<“dc<*’“sc’ '®dc(‘'-"’“sc» 


and similar definitions for the autocorrelation functions Rrii] 

Rifids^^’^sc^ ■ 

r, f \ A f” r n \ ^ ^ (' 

^—00 

Integrating (30) between -« and ~ and noting that 

r“ 9 /^n\^ /■” 4 ri/i 

R^^(t) dT = {-f) Gg(jco) ^ 


Rrt.(T) Rg(r) dT = f S M G(,(j^) 2 E; ,- = 2,3 

' — cn “• po 

■” ^0 f” 4 


■to /“ ^ H 1 

RiT|3(t^) S„^2(w) GQ(joj) ^ 


Sni "" power spectral density of m-|(t) cos u^^t 

1 c 

(tjj;w ) = power spectral density of s.(t) cos uj t 

^dc 


etc., 

we get an expression for given by 
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z:zr:zzj:zz 


# 


n 




= 4 N 


0 


Nq 


Gq(jw) 






Gq( jw) 


,03 


+ P, 


+ P 


Sni3(>^) 


v(-) 

d 


Gq( joj) 


duj 

271 




+ p, 


du) 

2ti 




GoEj(“-“5c>^, 


doj 

2ir 

4 


. . ^4 

+ 8 (vJ ^2 ^3 ^»I3(“) 

'—00 


Gq( jw) 


do) 

2it • 


duj 

2tr 


(35) 


In arriving at (35), v/e have used the additional relations 


(36) 


where S|^^(w) and Smj(u)) are the power spectral densities of m^(t) and s^(t), 

respectively. Also, we have taken note of the fact that Sm (w), Sm.(tj), 

4 1 ^ 

and jGQ(jw)[ are even functions of w. 

Finally, from (27), the equivalent loop signal -to-noise ratio , p', 

is 


cos^^^c 0'[P2D3-P3D3]f , „ _ 

“ . .4 ’ 


N B, 
e L 


(37) 


where p = P/NqBj^ is the loop signal -to-noi se ratio of a phase-locked (CW) 
loop operating on the total power P M-j + P 2 + P^ + P^j » and 


4N(,P 


2 /4 

cos (p - 

^C \7I 


^2 ”^2 ' ^3 ”^3 


(38) 


is the loop squaring loss . Defining the modulation indices n^> i =1,2,3; 
and H(j by 


Pi 


i = 1 ,2,3 
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( 39 ) 
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we obtain, after much simplification, 


2 2 
_■ (^) [’'2°2-"3°3^} 


where 




Smo(‘^) Sm^(co) Gn(jo)) 


" (tt) ^2 j 2 tt 

^—co 


f . h3 J O^IM f 


4 H 

4] S|,i^(u)) GQ[j(a)-a)^^}] ^ 


f” ! 

+ GQ[j(oj-u)^^)]! 


and p. = 2 P/NqB. is the signal -to-noise ratio in the "arm filter" bandwidth* 


defined by 


tA 1% 


G„(j«i)r d«3 • 


Assuming arbitrarily that ^ '*^^en (40) can be rewritten as 

(cos [Ho Dp “ Ho Dn] \ 

^ ^ 'i/^2,p/4f fh\^3 ■ 

2P^/Nq Hit / *^2 43\'Nq j Tp 


The tracking phase jitter performance can now be specified, In 
the linear region, the variance of the loop phase error 2 o? is given by 


Note that the bandwidth defined in (42) actually characterizes the 
cascade of the subcarrier loop arm filter and the lowpass equivalent of the 
input bandpass filter. 




13 


~ p' “ ’ 

v/here p' is given in (37). However, since the demodulation reference 
signals if (3) are at u rather than 2(u , then the tracking jitter 

on the output data streams is 

2 1 2 J_ , 

4 ^2(p p 5^ ' ' 

We note that S_, {(o) is similar to S (co) except that it is spread 
md 

slightly at the high frequencies because of the multiplication of m-j(t) 
by S 2 (t) and S 2 (t). However, since « 27r/T-| , then and Spp|^(u)) 

are approximately the same in the vicinity of the subcarrier frequency. 
Thus, to a good approximation. 






k = 2,4 


Additional practical assumptions can be invoked to simplify (44) still 
further. Since the bandwidth of Gg(jo)) is typically selected on the 
order of the larger of the two data rates l/T^ and l/T^ and since 
1 /T^ » I/T 2 or I/T 3 , then 

(w) GQ[j(a)-o)^^)] (to^^) I GQ(jw) ^ 


^m-| *^i ’ 

As an example, suppose nii(t) is an NRZ modulation of rate /?-| = l/T-j 
Then, 


S,„ (») = T 


/sin wT -|/2 

1 ( wT .,/2 ■ 


sin 


1 \ 


For a subcarrier frequency of 8,5 MHz, and data rate 1/Tn 100 Msps, 


In mode 1 of the Space Shuttle Orbiter-TDRSS Ku-band return link, 
iii-j(t) is a rate 1/2, constraint length 7, convolutional code with data 
rate 100 Msps. 


.1 


1 




lOOx 10 


/sin ti(8.5/100)\ 
\ 7t(8. 5/100) ) 


0.988 
lOOx 10^ 


The numerical evaluation of the subcarrier tracking jitter (44) 
for a fixed ratio of arm filter noise bandwidth to loop noise bandwidth 
is shown in Figures 2 through 5. Three assumptions are made: 

(a) The carrier loop is tracking perfectly; thus, cos cp^ = 1. 

(b) The bandwidth of H(w) is much wider than that of arm 
filters G(w). 

(c) The arm filters G(to) are assumed to be one-pole Butterworth 
(RC) filters or 

|e(j2’rf)l^ = ^ 2’ (“'S' 

where f^ is the cutoff frequency. The two-sided noise bandwidth (42) of 
G(j 2 tt f) is given by 

B. = irf^ . ' (49) 

Figures 2 and 3 illustrate the case when m.|(t) and ni 2 (t) are NRZ, while 
m 2 (t) is a Manchester code. Figures 4 and 5 illustrate the case when 
m^(t) is NRZj v/hile both m 2 (t) and ni 2 (t) are Manchester codes. The even 
figures have the signal-to-noise ratio (Pl' 2 /Ng) as a parameter, while 
the odd figures have the high subcarrier data rate as a parameter. 

It is observed that the changes in the subcarrier tracking jitter 
as a function of B^./i ?2 are more obvious when m 2 (t) and mg(t) are both 
Manchester codes than when iTi 2 (t) is NRZ and m^^t) is Manchester. 


'L radians) 










i ^ Figure 3. Subcarrier Tracking Jitter versus Ratio of Arm Filter Bandwidth 

to High Subcarrier Data Rate a parameter; mJt) and 
in^Ct) are NRZ» m^Ct) is Manchester'^code; 

■ 

” -7 

: i 




radians} 



radians) 


B,/B^ 


0 Mbp 
192 kbps 

10'^ 


(Perfect carrier tracking) 


Figure 5. Subcarrier Tracking Jitter versus Ratio of Arm Filter Bandwidth 
to High Subcarrier Data Rate/? 2 l a parameter; (Tii(t) is NRZ, 

m 2 (t) and in 3 (t) are Manchester cocfes; i?-j 
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APPENDIX I 

ERROR PROBABILITY PERFORMANCE OF UNBALANCED QPSK RECEIVERS 

by 

Marvin K. Simon 

1.0 INTRODUCTION 

Previous analyses by the author [1-5] and others [6-9] discuss 
the tracking performance of unbalanced QPSK receivers with particular 
emphasis on determining the mean-squared tracking phase jitter on the 
reconstructed carrier (or subcarrier) reference signal. As is true in 
any phase-coherent receiver, this lack of perfect coherence between the 
received signal and the reconstructed reference causes an additional 
degradation in system error probability relative to that due to the 
additive channel noise alone. Frequently, this increase in error proba- 
bility due to noisy synchronization references is translated into an 
equivalent increase in signal power required to produce the same error 
probability as obtainable in a perfectly synchronized system. This 
increase in required signal power is referred to as the noisy reference 
loss . 

Evaluation of noisy reference Toss in BPSK, QPSK, and offset 
(staggered) QPSK systems has been extensively covered in the literature 
[10-15]. While the techniques used there are certainly applicable to 
demodulation of unbalanced QPSK, the complexity of the evaluation when 
the ratio of data rates in the two channels is large (this will be 
explained in more detail in the next section) prompts one to look for 
a simpler calculation procedure. Indeed, such an approach is possible 
when the noisy reference loss is small or equivalently the effective 
signal -to-noise ratio in the tracking loop bandwidth is large, i.e,, 
the loop operates in its so-called linear region. Making such an 
assumption for purposes of error probability performance evaluation 
is quite reasonable when one realizes that this very same assumption 
has already been implied in the previously mentioned tracking perform- 
ance studies of unbalanced QPSK. The approach taken is to expand the 
error probability conditioned on the loop phase error <p in a power 
(Maclaurin) series in ^ and then, keeping only the first few terms of 
this series, average this conditional error probability over the proba- 
bility density function (p.d.f.) of By doing this, we obtain the 


additional error probability due to noisy synchronization references as 

an additive term directly proportional to the mean-squared phase jitter 
2 

directly associated with the receiver's tracking loop. Similar 
arguments can be advanced to give us closed-form results for the noisy 
reference loss itself. 

Using the above approach, this appendix presents results for the 
error probability and noisy reference loss performance of the unbalanced 
QPSK receivers corresponding to various combinations of data formats in 
the two channels. Specifically, for the subcarrier tracking loop asso- 
ciated with the two lower rate channels of the three-channel Space 
Shuttle Orbiter (SSO) Ku-band return link, wherein the two data channels 
are typically 192 kbps Manchester coded data and 2 Mbps NRZ data, 
numerical evaluation of this performance in these two channels is given. 

2.0 ERROR PROBABILITY PERFORMANCE FOR SMALL NOISY REFERENCE LOSS 

Common to all of the Costas-type demodulators of unbalanced QPSK 
previously considered [1, 2, 4,5, 6,8, 9] is the fact that the input 
signal plus noise is first demodulated with in-phase and quadrature 
phase detectors, the outputs of which represent the inputs to the data 
detectors for the two information channels. An example of this state- 
ment is illustrated in Figure 1 for the case where a conventional Costas 
loop is used for tracking the unbalanced QPSK signal; We shall pursue 
this case in detail and later indicate how the results can be extended 
to other Costas-type demodulators. 

The input x(t) to the Costas loop of Figure 1 is the unbalanced 
QPSK signal plus noise, namely, 

x(t) = )'^2P^ It-] (t) cos ii>(t) + v^P^ m 2 {t) sin <J>(t) + n^(t) , (1) 

v/here <J>(t) = w^t + OQ, with the radian carrier (or subcarrier) frequency 
and 0Q the input phase (assumed time independent) to be estimated, m-j(t) 
and m 2 (t) are independent data modulations assumed, in general, to be 
asynchronous, and n^.(t) is the additive channel noise which can be 
expres;sed in the form of a narrowband process about the actual frequency 
of the observed data, i.e., 

n^.(.t) = /2 (N^(t) cos ^>(t) - Ng(t) sin <l)(t)} . 


( 2 ) 
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In (2), N (t) and N (t) are approximately statistically independent, 
stationary, white Gaussian noise processes with single-sided noise 
spectral density Nq w/Hz, and single-sided bandwidth Bj^<WQ/2n. 

Assuming that the input in-phase and quadrature phase detectors 
are of the analog type, then when x(t) is demodulated by the quadrature 
reference signals 

r^(t) = ^ sin $(t) 

r^(t) = ^ cos $(t) , (3) 

the corresponding phase detector outputs (ignoring second harmonic 
terms) are* 

^^(t) - x(t) r^(t) = cos ip{t) - /^m^(t) sin (p{t) 

- N^(t) cos (pit) - N^(t) sin (p[t) 

E^(t) = x(t) c^(t) = /P^ m 2 (t) sin cp[t) + m^ (t) cos (p[t) 

- Ng(t) sin (p{t) + N^(t) cos cpit) , (4) 

where <;!?(t) = qi(t) - $(t) is the loop phase error. In addition to their 
forming a loop error signal when multiplied, these phase detector out- 
puts also, as previously mentioned, serve as inputs to the channel data 
detectors which we shall assume are matched filters. The outputs of 
these matched filters at the end of the (k-1) signaling interval in 
each channel are respectively given by 

1 fkTp+ep 

Qp = 7" E_(t)Pp[t-(k-l)Tp-ep] dt 

^ ^2 Mk-DTg+ejf 2 2 2 


/•kTi+e-j 


(k-l)T^+G 


e^{t) p^ Ct - T] - dt , 


( 5 ) 


v/here p^.(t); i = 1 ,2, is the basic unit power symbol pulse in the ith 
data stream ni^.(t) defined to be nonzero only on the interval (0,T^.), 


gain. 


r 

For simplicity, we shall assume these phase detectors have unit 


5 




and e^.i i = l,2, is the arbitrary pulse epoch in these same modulations 
Obviously then, the channel modulations m^.(t); i = l,2, are related to 
p.(t) and e. by 




m^(t) = I p.(t - nT. - e^.) ; i = 1 ,2 

n=_oo 


( 6 ) 


where i = 1,2,. are independent ±1 sequences with the properties 


m 


a . a . = 5 

ni mi mn 


f’ 


m = n 


0 ; ni 5 ^ n 


a . a . 

ni mj 


0 all m,n 

i 


(7) 


where the overbar denotes statistical expectation. 

Substituting (4) into (5) separates the matched filter outputs 
into their signal and noise components, i.e.,* 




aj^_.| 2 cos (p - /p 7 irii9 sin <j) - cos <t> - sin (p 


1 '"12 


s2 


’c2 


a |^_1 1 cos sin + cos ip , (8) 


where 


m 


A 2 

12 “ T. 


|>kT2+e2 

( k-1 )T2 +e:2 


m^ (t) P2[t - (k-1) T 2 - G 2 ] dt 


m 


A ± 

21 " T, 


^kl^+El 

(k-l)T.,+E^ 


ni2(t) p-|[t - (k-1 ) - e.|] dt 
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A 1 f'^i^^i 


Si T. 


N„(t) p.[t - (k-1) T. - E_.] dt ; i = 1 ,2 
(k-l)T,+c, ^ T ^ ^ 


1 1 


^ci T. 


. 1 /-kite. 

A 1 ( 1 1 ^ 

(k-l)T,+e ^ 


(t) p.[t - (k-1) T^. + e.] dt ; i = 1,2. (9) 


1 1 


We shall hrop the dependence of (p{t) on t and further assume 
that (p{t) = Is constant over a symbol interval of either data rate. 


6 


€> 


The data estimates a|^_i -j and aj^ ^ 2 the Channel 1 and Channel 2 
data bits aj^_-| i and aj^_i 2 transmitted, respectively, in the (k-l) sig- 
naling interval are then given by 


^k-1 ,i " ^^i^ V i = '• >2 


( 10 ) 


where sgn x = x/|x[ is the signum function. Note from ( 8 ) that the pre- 
sence of a loop phase error due to noisy carrier synchronization refer- 
ences causes crosstalk, i.e., the data in one channel affects the 
detection of the data in the other channel. This crosstalks which is 
proportional to sin <fi, in combination with the degradation of the desired 
signal amplitude by cos tfi together account for the noisy reference loss 
associated with the detection of an unbalanced QPSK signal. 

Letting 

^2 ^ " *^s2 “ ^c2 


N, = - N , sin A + N 1 cos A 
1 si ^ cl 

then, from (9), one can immediately show that 

N, 


2 4 - "0 

"i - - ?f7 


i = 1,2 . 


( 11 ) 


( 12 ) 


Thus, the probability of error in each channel conditioned on a fixed 
phase error is given by 

P|,_((;)) = Prob 5 ^ aj^_^ ^ ; i = l,2 


(13) 


which, from ( 8 ), (TO) and (12) becomes 


P^^(ci.) = i erfc 


' Vl - . 

cos <j> + — rn^T sin <i> 


•0 


0 


‘21 


+ j erfc 




cos <p 


IP,T, 


m 2 i sin (j) 


(14a) 
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<• 


Pr (^) 


erfc 


^^2^2 ri '^2 . . 

cos ^ + /-|r-mi 2 sin ^ 


'0 


0 


(14b) 


05) 

and the overbar denotes averaging over the data sequence and uniformly 
distributed random epoch dependences in m 2 i and rn-j^* We shall assume 
throughout the remainder of this report, as was done in [1,2], that 
the data rate i?^ - l/T^ of ni^ft) is arbitrarily chosen to be higher 
than or equal to the rate = 1/T-j of m.|{t). Then, when the ratio 
of these rates is large, the statistical average over the data 

sequence required in Pei(<{>) of (14) involves a comparably large number 
of complementary error function evaluations since m 2 i must be computed 
for each possible sequence of m 2 (t) in the interval (k-1 )T-j+E.j ^t <.kT^+e^ . 
Furthermore, for each of these contributing terms to P^^((|)), we must 
perform an average over the p.d.f. p(^) of the loop phase error (j) to 
determine the avev-age error probability performance of Channel 1. More 
will be said later about this additional average over ij>. 

In the meantime, to get around the computational bottleneck, v/e 
propose a scheme based upon expanding Pe^‘{<|)) into a power series in 
which, for small provides credible results. In particular, we note 
from (14) that P^.{(J)) is of the form 

Pe;(<I>) - erfc [A^ cos (f. + C- sin <j)] +1 erfc [A- cos sin ; 

1 = 1,2 (16) 

where A. is a constant and is a random variable in the sense that 
1 r 

it is both data sequence and random epoch dependent. Differentiating 
(ifi) once and twice with respect to ((> gives the results 


+ erfc 




cos <{) 


- . 
m^2 Sin ({) 


N 


0 


where the complementary error function (erfc x) is defined by 


erfc X = ^ 


/IT 


exp dy 




(C. cos - A. sin <P) exp [-(A. cos <P + sin <})) ] 


+ (C. cos <P + A. sin i|)) exp [-{A. cos <f) - sin (ji) ] 

2>^ ^ ^ ^ V 


(A., cos (t> + sin <|i) exp [-(A^ cos <j> + ?• sin <}>) ] 


+ 2 (A. cos iji+ sin cos <{) - A. sin (j)) exp [-{A. cos i|)+ sin cf)} ] 


4- (A. cos (j) - c,- sin exp [-(A. cos ^ . sin (Ji) ] 


+ 2(A. cos (l>- sin cos (ji + A. sin cj>) exp [-(A. cos <l>- sin tj.) ] 


which, when evaluated at (|) = 0, become 


d Pe^. ((|)) 
d ij) 


exp (-A?) + exp (-A.^) 

2/ir ^ ^ 2vV ' ' 


d^ Pe.(<i>) 


-V~ = “ tA. exp (-Af ) + 2A_. exp (-a/)] 

^ exp (-Af) p ^ ^ 

/tT 


Furthermore, from (16) 


Pp (f ) = 7 erfc A. 

i <fi=0 ^ 

Thus, combining (18) and (19), we get 
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it 


Pp M 


Pc (0) + 


PEi(<l>) 


d <f) 


4>=0 


(j) + 


PS 


d tf>' 


j erfc A^. + 


<p=0 


A. exp (-Ar) „ 

^ [1 + 2 ?/] 
2/^ ^ 


( 20 ) 


Comparing (16) with (14), we can immediately identify A- and in terms 
of the receiver parameters, whereupon (20) becomes 




Pp (<}>) = ^ erfc ^ exp (-R^) [1 + 2 Yj 


1 


P^^((|)) - 2 /^ + 2 J 


exp (-R 2 ) 


1 + 2R 


'm 


12 


n Yi 




( 21 ) 


with signal-to-noise ratios in each channel defined by 

A P • T. 

^ _J — L . i = 1 2 


(22) 


and the ratio of data rates defined as in [1,2] by 


'T 


c 

i?. 


L 

T. 


(23) 


# 


Note from (21) that, as promised, the error probabil ity consists of a 
term representing the performance for perfect carrier synchronization 
plus a term proportional to the crosstal k. 

What remains insofar as evaluating P[;^((|)) is to compute the 
normalized mean-squared crosstalk terms m^^^ and m-j ^2 f°*' various com- 
binations of data formats in Channels 1 and 2. Evaluation of m^^-, and 
m -]2 depends not only on the specific data formats in the two channels 
but also on the relative synchronism of the two data streams. More 
often than not, the two data clocks will be unsynchronized since, typi- 
cally, the two data streams are generated from totally independent sources. 
In this situation, there is, in general, no relation between the epochs 
E| and ^2 of the synchronization pulses corresponding to these clocks 
and, furthermore, the ratio of the two data rates need not be integer. 
Recognizing, however, certain similarities between the definitions of 
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-22 

and m-j^, and the computation of cross-modulation distortion 0^2 
associated with the Costas loop given in [1], it is relatively straight- 
forv/ard to show that 


~2 

"^12 

7F 

‘^21 




(24) 


where 0^2 given by (42), (50) » (51) and (52) of [1], and is also 
defined in that reference. These results are summarized in Table 1. 



Finally, averaging Pf^.((|)) of (21) over the p.d.f. of as deter- 
mined by analys.is of the Costas loop (see [1]) and assuming that the 180° 
phase ambiguity is perfectly resolved, we obtain the average error proba- 
bility performance of each channel, namely, 

p/2 

Prr = Pp (<l>)p(<l’)di> . (25) 

^i tt/ 2- ‘'i ' 

Substituting (21), combined with (24) , into (25) gives the desired simple 
result 
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*1 1 j •'I 

4r erfc ^ ^^2 '^f2^ 


'] 


^Eo " 2 zJ ^ ^“^2^ 


1 + 2 R 


^12 


(26) 


Alternately, in terms of the total signal-to-noise ratio R-j.^ the high 
data rate bandwidth and the transmitted modulation indices n-j and 
defined by 

_ (P, + P2)T2 ^ _ Pj 


T2 


P + P ’ 
^1 ^2 


i = 1,2 


(27) 


the results of (26) become 


1 ^T2 ^1 


2 J^T2 ^T ^1 2 J ' tt' ' (- Ry2^T ’^l^ ^T2^2'^12^'^ 


^Eo ” 2 n/^T2 ^2 2 


1 /'^T2'^2 


exp ( “ Rj 2 ^2 ^ ^ ^ ^ ^T2 '^l ^^12^ ^<f> 


(28) 


The mean squared phase error af may be related to the receiver 


parameters already defined by noting that, for the Costas loop of 
Figure 1,* 


1 




pS, 


(29) 


where p = (^i ■*■ ^2^'^'^O^L equivalent total signal-to-noise ratio 

in the single-sided loop bandwidth of a linear loop, and is the 
loop "squaring loss" which reflects the fact that the Costas loop error 
signal is formed by a nonlinear operation, namely, multiplication of the 
in-phase and quadrature filtered phase detector outputs. In terms of 
R .^2 defined in (27) and the ratio of loop bandwidth to higher data 
rate the loop signal-to-noise ratio p can be expressed as 

Rt2 

P = 0-41- . (30) 




Actually, the relationship in (29) is not restricted only to 
the conventional Costas loop of Figure 1 but also holds for a wide 
variety of other types of unbalanced QPSK demodulators (see [4,5,8, 
9]). 
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An expression for can be obtained from Eq. (28) of [1], namely, 

,2 


{r\2^2 


"I 


12 


"^1 °1 \ ’^2 ^2 ^2 2 * Rj 2 ^ *^1 *^2 T2 


(31) 


where the parameters D-j , D 2 J , Kq 2 » ^nd D-j 2 depend on the magnitude 
squared of the arm filter transfer function lG(j2-rrf)]^ and the power 
spectral densities Snn (f) , Sni 2 (f) of the modulations mT(t), m,(t) 


according to the relations 
D 


S^(f) lG(j27Tf)r df ; 
"’k 


k= 1 ,2 


D 


A 


12 


S„(f) S^(f) |G(j2TTf)r df 

'_cn .'"1 ^2 


C,S„,|^(f) |G(J 2itf)rdf 


k=l,2 


C |G(j2»nf tif 

r |G(j2,f)|2 df 


(32) 


Specific closed-form expressions for the above parameters, when the 
Costas loop has single-pole (RC) arm filters and the modulations m.j(t) 
and m' 2 (t) are variously NRZ and Manchester coded data streams, can be 
found in [1 ] . 

3.0 EVALUATION OF NOISY REFERENCE LOSS IN EACH CHANNEL 

Letting denote the value of when = 0 (i.e., ideal PSK 
performance) and approximating the complementary error function of (15] 
by the first term of its asymptotic expansion, i,e., 

effc X = 

/7x 


then, from (26), we have 
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1 + [1 + 2 R 2 o ^ 


<f> 




= 1 + R, 


1 + 2 R 



(34) 


or 


log^O Pe, - ^E, = '°9l0 + R, H + 2 l?2 '’/I 


^£2 " ^°^io ^£2 ~ ^°^io r ^^2 


1 + 2 R 



< (35) 


For small noisy reference losses and small error probabilities, a curve 
^°9 iq P£^- versus R^. in dB would be parallel to the corresponding 
ideal performance curve, namely, log^^ versus R^. in dB. Thus, 
letting r^. be the dB equivalent of R^., i.e., 


= 10 log.|Q R^. ; i = 1 ,2 


(36) 


then 


dlog^Ei ~ ■'°9lo''Ei - ■>O910'’Ei 


dri 


r : - r . 


(37) 


where 


r; = 10 log,g r: 


(38) 


and R^ is the signal -to-noise ratio required to achieve an error proba- 
bility Pe^. in the presence of the noisy synchronization reference. Thus, 
the quantity r.1 - r. represents the noisy reference loss (in dB) at the 
given value of error probabil ity Pe* . Jo compute an expression for it, 
we proceed as follows. Taking log^^Q of Pe^ and differentiating with 
respect to r^. gi ves 


d log,(,PE, 


1 /R~ exp (-R. 


d r^. 


1 O/tT 


erfc 


(39) 






(40) 


Using the approximation of (33) in (39) gives 




0.1 R. 


Finally, combining (35), (37), and (40) and solving for the noisy refer- 
ence loss - r^. in Channel i; i = T,2 gives 




(41) 


4.0 NUMERICAL RESULTS 

As an example of the application of the previous results, con- 
sider the performance of the two lower rate channels of the three-channel 
SSO Ku-band return link wherein m-j(t) is a Manchester coded data stream 
at i?i = 192 kbps and m2(t) is NRZ data with a maximum rate - 2 Mbps. 

We further assume that the power allocation is chosen so that, for the 
given data rate ratio Yj = ^ 10.42, the signal -to-noise ratios 

R.J , in the two channels are made equal, i.e., both channels operate 
at the same error rate. Thus, R-| = R2 implies P.jT^ = P^T^ and the modula- 
tion indices n-|, 02 become 


^1 " P-i + Pz ^ 

A '"z _ ^T 

*^2 - P^ + P 2 ~ 1 + Yj 


0.0876 

0.9124 


(42) 


Further, it is typical to design the Costas loop bandwidth on the order 
of i?2/100 (or less) since most of the power is in the high rate channel 
which controls the performance of the tracking loop. Thus, assuming 
B|^/i?2 = 0.01 ,* Figures 2 and 3 illustrate L-| and L2 of (41) versus B^/i?2 

Smaller values of BL/i?2 as would be typical in practical receiver 
design v/ould yield insignificant losses in L-| and L2. 
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“ 0.18 


0.16 


P^ = io- 


= 1 O' 


Figure 2, Channel 1 Noisy Reference Loss (in dB) vs. the Ratio of Two-Sided 
Costas Loop Ann Filter Bandwidth to Channel 2 Data Rate; m,(t) is 
192 kbps Manchester coded data, ni 2 (t) is 2 Mbps NRZ data. 


sbs 








Figure 3. Channel 2 Noisy Reference Loss (in dB) vs. the Ratio of Two-Sided Costas 
Loop Arm Filter Bandwidth to Channel 2 Data Rate; m.|(t) is 192 kbps 
Manchester coded data, nu(t) is 2 Mbps NRZ data. 
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-4 -5 -6 

for error probabilities of 10 ,10 , and 10 , corresponding respec- 
tively to R] “ ^2 ” Several conclusions may be 

drawn from these figures, First, the noisy reference loss on Channel 2 
is considerably smaller than that on Channel 1, The principal reason for 
this can be easily explained in terms of the result in (41), where it 
is observed that the effective cross-modulation loss on Channel 2, 
is divided by Yj which in this case has a value equal to 10.42. Secondly, 
for either channel, the noisy reference loss decreases with increasing 
error probability. This is intuitively satisfying when one realizes 
that the slope of the error probability versus signal -to-noise ratio 
curves becomes steeper as becomes smaller and thus, for a given a^, 
the parallel ideal and noisy sync error probability curves become closer 
together. Finally, we observe that there exists an optimum arm filter 
bandwidth (for fixed i?„) in the sense of minimizing L . , i = l,2. Since 

2 ^ i 

only depends on this bandwidth, it is clear that this bandwidth 

9 2 

choice is identical to that which minimizes or, equivalently, the 
loop squaring loss. Note that, if is decreased, then the noisy 

reference loss will also decrease, since the equivalent loop signal - 
to-noise ratio p of (30) increases. 

We conclude by noting that, while the results of this appendix 
have been directed principally toward the demodulation of unbalanced 
QPSK by a conventional (single-channel) Costas loop such as in Figure 1, 
the expressions for average error probability [see (26) and (28] and 
noisy reference loss [see (41)] apply in a much broader sense. In par- 
ticular, the two-channel type Costas loops discussed in [8,9] have a 
mean-squared phase jitter given by (29) where, however, is a much 
more complex function of the various system parameters such as data rates 

and channel power ratios. Nevertheless, once S, (and thus ) is deter- 

L (J) 

mined, (26) and (41) apply directly toward evaluation of the noisy syn- 
chronization reference effects of these loops on error probability 
' performance. Other possible applications of the results pertain to 
demodulation of unbalanced QPSK using a biphase Costas loop with 
switching type multipliers. Once again, (26), (41) and (29) apply, 
provided that an expression for the squaring loss can be found (see, 
for example, [5,16]). 


5.0 CONCLUSIONS 


This appendix has demonstrated a simple technique for calculating 
the error probability performance and associated noisy reference loss 
of practical unbalanced QPSK receivers. The result for error probability 
is in the form of a leading term representing the ideal (perfect syn- 
chronization references) performance plus a term proportional to the 
mean-squared crosstalk. For the Ku-band return link subcarrier demodu- 
lation of 192 kbps Manchester coded data and 2 Mbps NRZ by a conventional 
Costas loop, it is concluded that the crosstalk degradation due to noisy 
subcarrier demodulation references is quite small (on the order of 
tenths of a dB or less, depending on the particular channel ar.d the 
ratio of loop bandwidth to data rate in that channel). When the higher 
data rate channel is 1 Mbps Manchester coded data, then since both 
channels are now Manchester coded, the crosstalk loss would be even 
smaller yet (see Table 1). The general results obtained can also be 
applied to other unbalanced QPSK receivers with more complex CostaS" 
type loop structures. 
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PRACTICAL DESIGN CONSIDERATIONS ASSOCIATED WITH OPTIMUM 
CARRIER RECONSTRUCTION TECHNIQUES FOR UNBALANCED QPSK 

by 

Marvin K. Simon 


INTRODUCTION 

In a recent paper [1], the authors address the problem of carrier 
synchronization of an unbalanced QPSK (UQPSK) signal format, starting with 
the well-known maximum a posteriori (MAP) estimation technique as motiva- 
tion for deriving closed loop tracking configurations. In particular, 
such carrier reconstruction implementations are suggested upon examining 
the gradient of the likelihood function whose solution is the MAP esti- 
mator of carrier phase. While the MAP approach leads to closed loop 
implementations with active arm filters which are matched to the signal 
pulse shape (integrate-and-dump circuits for digital signals) and thus 
require knowledge of the data timing clock, one can, in practice, replace 
these filters with passive lowpass filters with, however, an attendant 
loss in performance. Nevertheless, it is common practice to accept this 
performance penalty so as to avoid the necessity of having to lock up the 
bit synchronizer prior to achieving carrier lock. 

Before going into the specific details of the performance of these 
so-called "optimum" carrier reconstruction loops for unbalanced QPSK, it 
is perhaps worthwhile reviewing these same considerations first for simple 
PSK and then for balanced or conventional QPSK. In the latter case, we 
will come upon a new structure which, performance-wise, is theoretically 
identical to the well-known quadriphase Costas loop [2]. Also, as we 
shall see, the nature of the approximation to the true MAP solution, 
which allows us to arrive at this new configuration, has application in 
the unbalanced QPSK situation, thereby permitting us to carry the work 
initiated in [1] one step closer to an "optimum" solution. Indeed, the 
generalization proposed here of the linear in-phase channel configuration 
given in [1] allows carrier reconstruction from a UQPSK signal at all 
ratios of data rates and powers in the two channels, even in the 1 imit 
as these ratios simultaneously approach unity, i.e., balanced quadriphase . 

Finally, we shall discuss the practicality of the unbalanced QPSK 
carrier reconstruction loops found in [1], paying particular attention 


to their sensitivity to variations in channel gains. It will be shown 
that, whereas the channel gains were "optimally" chosen in [1] from MAP 
estimation considerations, a different selection of these gains based 
upon directly optimizing the loop's tracking performance can yield as 
much as a 10 dB improvement in this performance. The significance of 
this statement is not so much the fact that the tracking performance can 
be enhanced by a better choice of gains, particularly since these gains 
are now theoretical ly signal-to-noise ratio dependent, but rather the 
high degree of sensitivity of the performance of the MAP estimation-type 
loop to variations in the channel gains themselves. This point will 
be e;<plored numerically to provide quantitative verification of this 
sensitivity behavior. 

CARRIER RECONSTRUCTION LOOPS MOTIVATED BY MAP ESTIMATION THEORY 
PSK Modulation 

Consider the problem of finding the MAP estimate of carrier phase 
from an observation made on a biphase-modulated, suppressed-carrier 
signal'in additive Gaussian noise. In particular, the signal s[t,e(t)] 
is of the form 

s[t,0(t)] = /^ m(t) sin [wgt + e(t)] , (1) 

where S is the average signal power, m(t) is a binary modulation (a ±1 
digital waveform) of rate i?=l/T and pulse shape p(t), Wg is the radian 
carrier frequency, and e(t) = 0 g+Qgt is the received carrier phase. 

The total received signal x(t) is then 

x(t) = s[t,0(t)] n.(t) , (2) 

where n^.(t) is the additive channel noise which can be expressed in the 
form of a narrowband process about the actual frequency of the observed 
data, viz. , 

n^(t) - /2 (N^(t) cos [a3gt+6(t)] - N^(t) sin [togt+ 0 (t)]} , (3) 

where N,(t) and N (t) are approximately statistically independent, sta- 
tionary, white Gaussian noise processes with single-sided noise spectral 
density Ng w/Hz (see [2]) and single-sided bandwidth B^ < a)g/2ir. 
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The MAP estimation problem may now be stated as follows. Based 
upon observation of x(t) over the interval 0^t:< KT and knowledge of S, 
m(t), Wq, and the precise time instants at which the modulation can change 
states, we wish to choose that value of 6(t) = e (assumed constant over 
Oit<.KT) v/hich maximizes the a posteriori probability p{e(t)lx(t)) . 

The solution to this problem, i .e. , the best estimate of 0 in the MAP 
sense, is well known [2] to be that value that maximizes the function 


K 

o 

rkT 

f(0) = yj cosh 


x(t) sin (w.t + 0 ) p(t) dt 
(k-l)T ^ 

k=l 


or, equivalently, 

K 

In f(9) = I In cosh 

k=l 


0 


kT' 


(k-l)T 


x(t) sin (u)Qt + e) p(t) dt 


. (5) 


The physical interpretation of (5) implies an open loop estimation in the 
sense that, for each value of 0 in the continuum (-Tr,7r), we form'the func- 
tion lnf{0) and then choose that particular value which corresponds to 
the maximum of this function. A discrete (in the sense of a finite number 
of trial values of e) version of this interpretation is illustrated in 
Figure 1 . 

To go from open loop to closed loop configurations, we consider 
an alternate interpretation of the MAP estimate solution. Defining 

rkT 


g(e) ^ 


K , 

I f 

k=l ‘^0 


X tanh 


(k-l)T 
"o rkT 


/2S x(t) cos (o)gt +. 0 ) p( t) dt 


'0 


(k-l)T 


x(t) sin (co„t + 6) p(t) dt 


( 6 ) 


then the MAP estimate § is also the solution to 

g(e) = 0. 


QY POOR (7) 


For any value of 0 other than e, g(§) will be either positive or negative, 
depending on the sign of e-0. Thus, g(e) is intuitively an appropriate 
choice for an error signal in a carrier phase tracking loop. Using this 
motivation. Figure 2 illustrates such a closed loop implementation which 
herein will be referred to as the MAP estimation loop . 
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The MAP estimation loop of Figure 2 is impractical, primarily 
because of the difficulty of implementing the hyperbolic tangent non- 
linearity. To arrive at practical realizations, one must first approxi- 
mate this nonlinearity with simpler, more easily implementable functions.- 
In particular, we have that 


tanh X = sgn x *, 

X large 

(8) 

tanh X = X ; 

X small . 

(9) 


Since the input to the nonlinearity is a monotonic function of signal- 
to-noise ratio (SNR), then the approximations of (8) and (9) correspond, 
respecti vely , to conditions of high and low SNR. Using these approxi- 
mations, then, in (6) leads to the practical implementations of Figures 
3 and 4. In Figure 3, the arm filters are, as in the MAP estimation loop, 
of the Integra te-and-dump type and thus require a knowledge of symbol 
sync. In Figure 4, these same filters have been replaced by lowpass 
filters and, in accordance with this change, the digital filter and NCO 
are replaced by an analog loop filter ancj VCO. The tracking performances 
of these familiar loops, i.e. , the conventional Costas loop (Figure 4a) 
and the polarity-type Costas loop (Figure 4b) or, equivalently, the Costas 
loop with hard-limited in-phase channel, have recently been documented 
in the literature [3,4] wherein the bandlimiting effects of the arm filters 
on both the input data modulation and the noise have been accounted for. 

In both cases, it is shown that, by properly selecting the arm filter 
bandwidth, one can optimize the loop's tracking performance in the sense 
of minimizing its squaring loss or, equivalently, in the linear region 
of performance, minimize its mean-squared tracking jitter. 

QPSK Modulation 

Quadri phase-shift-keyed (QPSK) modulation, as is well known, offers 
the opportunity of achieving a given bit error rate performance with half 
the channeT bandwidth required' by PSK. Carrier reconstruction loops for 


•jUf 

We shall see that, in the case of quadriphase modulation, this 
approximation of the hyperbolic tangent function for small values of 
its argument is not sufficient. 
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such a modulation are variously referred to as the f ourth-power loop and 
the quadri phase Costas loop [?}. Assuming analog (four quadrant type) 
multipliers, these two loops have been shown [2] to be stochastically 
equivalent. 

In this section, we derive the MAP estimation loop for a QPSK 
signal, and by suitable approximations to the nonlinearity which arises 
as a consequence of the MAP theory, reconfigure this loop with practical 
realizations which are valid for high and low signal-to-noise ratios. In 
particular, for the low signal-to-noise ratio case, we are led directly 
to a new structure which, like the fourth-power loop, will be shown to 
be stochastically equivalent to the quadriphase Costas loop mentioned 
above. This new structure, however, has decided implementation advan- 
tages over the conventional quadriphase Costas loop, particularly in 
applications where the receiver may be operated in either a biphase or 
quadriphase mode. The significance of this statement will become clearer 
after we have had an opportunity to examine this new structure in finer 
detail . 

As a generalization of (1), a quadriphase-modulated , suppressed- 
carrier signal takes the form 


s[t,9(t)] = m-j(t) sin [cogt + 0(t)] + ^ m 2 (t) cos [uugt + e(t)] , (10) 


where S is the average signal power in either of the two data channels, 
m^(t) and m^lt) are the corresponding independent, binary modulations, 
both (for balanced quadriphase) of rate /?=1/T and pulse shape p(t), and 
0(t), as previously defined, is again the received carrier phase. Based 
upon an observation of this signal in additive noise as described by (2) 
and (3), then by an analogous approach to that used in deriving (4), we 
find that the MAP estimate of carrier phase is that value that maximizes 
the function 


K I 2 fkT ] 

f( 0 ) = 77 cosh 1^ /2S x(t) sin (u„t + 0 ) p(t) dt 

k=l To Mk-1)T ^ 


X cosh 


■ rkT 


N. 


0 Mk-1)T 


x(t) cos (w^t + 0 ) p(t) dt 


( 11 ) 


or, equivalently, 


T ‘■'~H ■ ■ ' ' 
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K 

In f(0) = I 

k=l 


In cosh 


x(t) sin (o)„t+0)p(t) dt 


+ In cosh 


’0 


A. 

1^0 J(k-1)T 
kT 

x(t) cos (wr,t+ 0 ) p(t) dt 
(k-l)T ° 


(12) 


Once again, to go from an open loop MAP estimate to a closed loop track- 
ing configuration, we introduce the function g(0) = (dlnf( 0 ))/d 0 and use 
g(§) as an error signal in such a closed loop. Differentiating (12) with 
respect to 0 gives 
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kT 
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x(t) cos (ojQt+ 0) p(t) dt 
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x(t) sin {tUpit + 0) p( t) dt 


kT 


(k-l)T 


x(t) sin ((j-t + 0) p( t) dt 


rkT 


(k-l)T 


x(t) cos (o) t + o) p(t) dt 


(13) 


The corresponding MAP estimation loop for QPSK is illustrated in Figure 5. 

To achieve a practical realization of Figure 5 for high SNR, we 
again approximate the hyperbolic tangent function as in (8). The result- 
ing loop, using passive arm filters, is illustrated in Figure 6. If one 
tries to use the approximation of tanh x given in (9) to construct a low 
SNR implementation of the MAP estimate loop, then substituting (9) in 
(13), we immediately see that g(o) = 0 for all e; hence, no error signal 
is generated. To circumvent this apparent dilemma, we go one step 
further by approximating the hyperbolic tangent function by the first 
two terms of its power series, namely, 

tanh X = X - -y ; x small . (14) 

Substituting this approximation into (13) results in 
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g(®) = ^ ^ x{t) sin (tjj t+ 0.) p(t) dt 
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which is of the form 


rkT 

x(t) cos ((o„t+ e) p(t) dt 


x(t) sin (w^t + 0) p(t) dt^ 


g(0) = T i (AB^ - BA^) , ( 16 ) 

k=l 

where A and B are immediately identified from (15). Rewriting (16) as 

' • 

g(e) = I I [AB(B^-A^)] (17) 

k=l 

allows us to draw the practical realization of the MAP estimation loop 

for low SNR (see Figure 7).* Several interesting observations can be 

made from both (17) and Figure 7. The signal product AB represents the 

error signal generated in a conventional biphase Costas loop. The signal 
2 2 

difference B - A represents the output of a lock detector associated with 
the same biphase Costas loop. Thus, the quadriphase error signal of (17), 
namely, AB(B - A ), is the product of the error signal and lock detector 
output of a conventional biphase Costas loop . The accumulation required 
in (17) is, as before, accomplished in the analog version (Figure 7) by 
the loop filter F(s). 

Implementing a quadriphase carrier tracking loop in Figure 7 has 
several implementation advantages over the conventional quadriphase Costas 
loop. First, it allows for either biphase or quadriphase operation, since 
the biphase errorsignal AB is already available. This choice of operation 


This particular quadriphase Costas loop configuration was origi- 
nally suggested to the author by J. C. Springett of Axiomatix, apart 
from its motivation by MAP estimation theory. 









mode can be accomplished using a simple switch (see Figure 7). Thus, 
relative to the hardware required to implement a biphase Costas loop and 
its associated lock detector, all that is needed to generate a quadri- 
phase error signal is one additional analog multiplier. We hasten to add 
that the large SNR loop of Figure 6 can also, as shown in this figure, 
be made to operate either as a biphase or quadriphase polarity-type Costas 
loop. Second, even if the loop is designed strictly for carrier tracking 
of balanced QPSK, the implementation of Figure 7 requires only two quad- 
rature reference signals and two arm filters, as opposed to the four 
reference signals (spaced by ir/4 radians) and the four arm filters 
needed to build the conventional quadriphase Costas loop. 

The only point that remains is: How does the performance of the 

loop in Figure 7 compare with that of the conventional quadriphase loop? 
The answer to this question is that the two loops are stochastically 
equivalent , i .e. , it can be shown [5] that the two have identical sto- 
chastic differential equations of operation. Thus, by starting with 
MAP estimation theory as a mathematical basis, we have proven, via .the 
introduction of a new configuration, that the conventional quadriphase 
Costas loop and its equivalent (the fourth power loop) are low SNR prac- 
tical realizations of the MAP estimate loop for QPSK. 

Unbalanced QPSK 

With the results for PSK and QPSK as background, we are now in a 
position to better understand, interpret, and augment previously obtained 
results [1] for the case of interest here, namely, carrier reconstruction 
techniques for unbalanced QPSK (UQPSK) as motivated by MAP estimation 
theory. While it is true that, depending on the ratio of pov/ers in the 
two channels of the UQPSK signal, a simple biphase or quadriphase Costas 
loop can be employed for carrier reconstruction [6,7,8], the purpose of 
the study in [1] was to search for "optimum" carrier reconstruction tech- 
niques in the sense that the loop would make full use of the total power 
in the two channels to enhance its tracking capability. Also, by applying 
MAP estimation theory directly to the UQPSK signal, the intuitive feeling 
was that some sort of hybrid configuration would arise which would present 
a continuous compromise between the extreme desires of tracking a highly 
unbalanced QPSK signal and the more common balanced QPSK signal in which 


17 


the powers and data rates in the two channels are equal. 

Indeed, the results derived in [1] largely represent this dream 
come true. Unfortunately, however, the low SNR realizations of the MAP 
estimate loop for UQPSK suffer from the fact that they cannot track a 
balanced QPSK signal. The problem herein lies again in the way in v/hich 
the hyperbolic tangent nonlinearity is approximated at low signal-to- 
noise ratios. We shall demonstrate shortly that, by using (14) rather 
than (9) at low SNR, one obtains a practical implementation which allows 
carrier reconstruction from a UQPSK signal at all ratios of data rates 
and powers in the two channels, in particular, in the limit as these 
ratios simultaneously approach unity, i.e., balanced QPSK. In fact, 
in this limiting case, the loop becomes equivalent to the quadriphase 
loop of Figure 7, 

An unbalanced QPSK signal takes the form 
s[t,0(t)] = m-|(t) cos [ojQt+ 0(t)] + m 2 (t) sin [a>Qt+ e(t)] , (18) 


where P^ and P^ are the average signal powers, respectively, in channels 1 
and 2 , and m^(t) and m 2 (t) are the corresponding independent binary modu- 
lations with rates i?-| = 1/T^ , = l/T^ and, in general , different pulse 

shapes p-j(t) and P 2 (t). Applying the MAP estimation theory to an obser- 
vation of the signal (18) in noise, the error signal g(e) for the MAP 
estimation loop is derived in [1] as 
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x(t) cos (oj^t + 0 ) p^(t) dt 
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/2P^ x(t) sin (oj^t + 0 ) p.| (t) dt 
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X tanh 
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( 19 ) 


In (19), T|^ ( k = 1 ,2, . . . ,K-| ) Is the ordered set of time instants at which 
the modulation m^(t) may potentially have a symbol transition in the 
observation interval (0,T), Similarly, 6^ (k = 1 ,2, , . . ,K 2 ) is the ordered 
set of time instants at which the modulation nigCt) may have a symbol' 
transition in the same observation interval. Note that, since the two 
modulations might, in general, be generated from independent data clocks, 
no restriction is placed on the relative synchronization between the taus 
and the deltas. This same assumption could also have been made for the 
balanced quadriphase case with the appropriate modification of (13). 

Illustrated in Figure 8 is the MAP estimation loop for UQPSK with 
g(§) of (19) as its error signal. A practical realization of this loop 
for high SNR, obtained by using the approximation to the hyperbolic tangent 
function of (8) in (19), is illustrated in Figure 9. For small SNR, using 
(9) as an approximation to tanhx, the authors of [1] suggest the loop 
illustrated in Figure 10. Unfortunately, this loop (or its equivalent 
using active arm filters of the type given in Figure 8) has the disadvan- 
tage that, as the ratio of rates arid powers both approach unity, i.e., 
balanced quadriphase, the two lowpass filters LPF-j and LPF^ would become 
identical (or equivaleatly integrate-and-dump filters of equal duration), 
and thus the error signal at the input to the loop filter goes to zero 
for all loop phase errors. If, instead of (9), one were to use (14) as 
an approximation to tanhx, then making this substitution in (19) results 
in the loop illustrated in Figure 11, where again we have drawn the case 
where the arm filters are of the passive type. Note that this two-channel 
Costas loop configuration reduces (except for the one-third gain factor) 
to Figure 7 when the input signal becomes balanced QPSK, i.e., LPF.j = LPF^ 
and, thus, it is capable of tracking such a signal. Rather than go into 
the details of the performance of the loop in Figure 11 at this point, we 
shall turn to an investigation of the sensitivity to gain variations of 
the simpler loop in Figure 10 and use these results qualitatively as being 
indicative of two-channeT Costas-type configurations. Actually, to make 
life even simpler, we shall perform our sensitivity analysis on the equiva- 
lent loop to Figure 10 where active (integrate-and-dump) arm filters of 
the type illustrated in Figure 8 are used in place of the passive arm 
filters. Making this substitution greatly simplifies the analysis, but 
nevertheless allows us to illustrate the conclusions we are after. 



Figure 8. The MAP Estimation Loop for Carrier Phase (UQPSK) 
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SENSITIVITY OF TRACKING PERFORMANCE OF TWO-CHANNEL 
COSTAS-TYPE LOOPS TO CHANNEL GAIN VARIATIONS 


The tracking performance of the UQPSK carrier reconstruction loops 
illustrated in Figures 9 and 10 has been studied in great detail in [1]. 
Part II of this reference treats the passive arm filter case, i.e., Figures 
9 and 10 as they are drawn, while Part I presents the results for the 
equivalent loops with active arm filters. In each case, the phase track- 
ing jitter based on a linear loop model is calculated as a function of 
such system parameters as signal-to-noise ratio in each channel, the ratio 
of powers in the tv/o channels, and the ratio of data rates in the two 
channels. Although the theoretical results were derived as a function 
of the ratio of the gains in the two channels, specific numerical results 
were presented only for the "optimum" choice for this ratio based upon 
MAP estimation theory. 

Unfortunately, as we shall see shortly, selection of this gain 
ratio as motivated by MAP estimation theory does not necessarily optimize 
the tracking performance. In fact, the selection of this gain ratio to 
directly optimize the tracking performance in the sense of minimum mean- 
squared phase tracking jitter, can yield as much as a 10 dB improvement 
in this performance relative to the MAP estimation choice of gain ratio. 

To illustrate this point, we shall begin by extracting pertinent 

results from [1] for the equivalent loop to Figure 10 with active arm 

filters and arbitrary channel gains and In particular, it is 

shown in [1] that the variance of the loop phase error (phase tracking 
2 

jitter 02 ^ for a linear loop model is given by 
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and we have implicitly assumed that the two data clocks are synchronized 
with the ratio of data rates being integer, i.e., 
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Equivalently, the loop "squaring loss" is given by 
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Clearly, from (20) and (23), minimizing the phase tracking jitter 
is directly equivalent to maximizing (minimizing the squaring loss). 
From the form of (23), it is a simple matter to differentiate this equa- 
tion with respect to K 2 /K-J , and thus find that value of channel gain 
ratio which maximizes After some routine algebra, the solution to 
the equation 
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Alternately, defining the ratio of powers in' the two channels by 
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pypi 

we can rewrite (25) in the form 
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where, in addition, we have defined (as in [1]) as the total power-to- 
noise ratio in the high data rate bandwidth, i.e.. 
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Note that the first factor in (27) is the "optimuRi" value of K^/K-j as 
motivated by the MAP estimation theory, i.e., 
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which is independent of signal-to-noise ratio. Only when the channels 
have equal energy, i.e., x=l/n, are the results of (27) and (29) equal, 
i.e., K 2 /K^ = x^. 

Substituting (27) into (23) results in an expression for the minimum 
squaring loss, namely, 
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On the other hand, using (29) in (23) gives the result obtained in [1], 
namely, 
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Figures 12 through 15 illustrate 5^ as computed from (30) and (31) versus 
the power ratio X with n fixed and = -3 dB and 10 dB. The values of n 
selected for these plots range from balanced (n -1 ) to highly unbalanced 
(n=10) data rates. Note that, for the balanced data rate case (n=l, 

X arbitrary) , both (30) and (31) reduce to the same result, namely. 












X=P2/P, 

Figure 14. Squaring Loss vs. Power Ratio for Unbalanced QP5K; Two-Channel 

' Costas Loop; n = To/T-, = 5 








A = P2/P, 

Figure 15, Squaring Loss vs. Power Ratio for Unbalanced QPSK; Two-Channel 

Costas Loop; n = T 2 /T-| = 10 
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despite the fact that the channel gain ratios (27) and (29) are not, in 
general, equal. Furthermore, the result in (32) is also identical to the 
squaring loss of a conventional (one-channel) Costas loop when used to 
track unbalanced QPSK with both unequal data rates and unequal powers [8], 
In this case, the identical in-phase and quadrature integrate-and-dump 
arm filters are assumed to be chosen to accommodate the higher data rate 
modulation. 

Several important conclusions can be drawn from the numerical 
results illustrated in Figures 12 through 15, First, since we have 
already observed. Figure 12 is representative of the one-channel Costas 
loop performance for all n and X, these curves represent lower bounds 
on the squaring loss performance of the two-channel configuration. That 
is, for fixed c-] and X, the two-channel Costas loop will exhibit a smaller 
squaring loss for any data rate ratio n than the value indicated in 
Figure 12. This performance comparison between one- and two-channel 
Costas loops will be made in the next section of this report. 

Second, the difference between the true optimum squaring loss 
[Eq. (30)3 and the value given by (31), as in [1], increases with 
increasing n . Also, for fixed n , these differences are much more 
significant at higher values of total power-to-noise ratio c-| in the 
high data rate bandwidth than at the lower values. As an example, for 
n=10 and ?-| = 10 dB, we see from Figure 15 that, for equal powers (x = l), 
the squaring loss as given by (30) is 10 dB smaller than that predicted 
by (31). The same comparison at = -3 dB only shov/s a 1.4 dB improvement. 

Finally, the true optimum squaring loss [Eq. (30)] decreases with 
increasing c-| , regardless of the values of n and A, except for the case 
of small X and n= T. Using the result of (31) for squaring loss, we 
observe that, -for fixed n , this "optimum" squaring loss can actually 
increase with increasing , depending on the value of X. 

To further illustrate the sensitivity of the two-channel Costas 
loop squaring loss performance to the particular choice of gain ratio 
K^/Ki or, more important, variations in this ratio about some nominal 
design value, e.g., x/n. Figures 16 through 18 illustrate versus 
K 2 /K^ for fixed values of n and , with X as a parameter. In these 
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Figure 16. Squaring Loss vs. Channel Gain Ratio for Unbalanced QPSK; Two 

Channel Costas Loop; ^.j = 10 dB, n=10 
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Figure 18. Squaring Loss vs. Channel Gain Ratio for Unbalanced QPSK 

Two-’Channel Costas Loop; Ci=10 dB, n=2 
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figures, K 2 /K -1 is allowed to vary from 0.1 (x/n) to 2(x/n) with the opti- 
mum values corresponding to (27) clearly indicated by the peaks of the 
curves. The expression used for plotting these curves is. a rewrite of 
(23) in terms of the definitions of n , X, and c-j , namely. 



where 



(33) 
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A COMPARISON OF THE SQUARING LOSS PERFORMANCE OF 
ONE- AND TWO-CHANNEL COSTAS LOOPS 

We have already mentioned the fact that one-channel (conventional) 
Costas loops are inferior to two-channel Costas loops when tracking 
unbalanced QPSK. The squaring loss performance for the one-channel 
loop (assuming both data modulations are NRZ) was given in (32). To 
allow a simple comparison with the performance of two-channel loops, we 
shall assume in the latter case that both channels have equal energy, 
i.e., equivalently, X = l/n. For this case, we have 

already observed that the value of K^/Ki selected by MAP estimation 
theory considerations is identical to that which minimizes the squaring 
loss. Thus, from (30) [or (31)], we find that 



Taking the ratio of (35) to (32) then gives the improvement in tracking 
performance of the two-channel Costas loop over the conventional loop. 



Letting denote this ratio, we have 
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Figure 19 illustrates 6j^ (in dB) versus X with c-j (in dB) as a parameter. 
Clearly, for X approaching unity, the improvement becomes infinite; how- 
ever, we recall that, in this limiting case, neither the one-channel nor 
the two-channel loop is capable of tracking at all, i.e., goes to 
zero for both. 


CONCLUSIONS 

In conclusion, we point out that the MAP estimation theory provides 
good intuition for implementing closed loop tracking configurations for 
BPSK, QPSK, and UQPSK modulations. In the latter case, however, care 
must be exercised in selecting the ratio of gains between the two channels 
of the resulting configuration. In this regard, two options are available, 
namely, that gain ratio which is motivated by MAP estimation theory and 
is independent of signal-to-noise ratio, and that ratio which minimizes 
the loop squaring loss and is dependent on signal-to-noise ratio. The 
difference in squaring loss between these two choices of gain ratios can 
be as much as 10 dB for certain values of data rate ratio and power ratio. 
When the signal energies in the channels are chosen equal, then both gain 
ratios also become equal and yield identical squaring loss performance. 

Even in this case, the two-channel Costas loop can considerably outperform 
the one-channel (conventional) Costas loop, However, neither the con- 
ventional nor the two-channel loop of [1] is capable of tracking balanced 
quadriphase. By an extension of the power series approximation used for 
the hyperbolic tangent nonl ineari ty which arises from the MAP estimation 
approach, we have been able to demonstrate a carrier reconstruction loop 
for UQPSK which should yield better performance than the above-mentioned, 
two-channel loop as the modulation becomes more balanced; in particular, 
it acts like a quadriphase Costas loop in the limiting case of balanced 
QPSK, A detailed analysis of the tracking performance of this loop is 
the subject of a future report. 
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ON THE EQUIVALENCE IN PERFORMANCE OF A PRACTICAL REALIZATION 
OF A MAP ESTIMATION LOOP FOR BALANCED QPSK AND 
A CONVENTIONAL QUADRIPHASE COSTAS LOOP 

by 

Marvin K. Simon 

1.0 INTRODUCTION 

In Appendix J (also [1]), the author considered the similarities 
in structure of carrier reconstruction loops motivated by MAP estimation 
theory for PSK, QPSK, and unbalanced QPSK modulations. In particular, 
for the balanced QPSK case, it was shown that, by approximating the 
hyperbolic tangent nonlinearity in the MAP estimation loop by the first 
two terms in its power series, an interesting practical realization of 
this loop results which applies at low signal-to-noise ratio (see 
Figure 1). Indeed, the error signal in this loop is formed by multi- 
plying the error signal and lock detector output signal of a conventional 
biphase Costas loop. We also note from Figure 1 that such a quadri- 
phase loop can be constructed using only a pair of quad.rature reference 
signals and a pair of arm filters, as opposed to the four reference 
signals (separated by tt/4 radians) and four arm filters required in 
a conventional quadriphase Costas loop (see Figure 2). The loop of 
Figure 1 also has the advantage that it can easily be switched from 
a biphase mode to a quadriphase mode depending on the form of the input 
modulation. With all this in its favor, the only question remaining is: 
How does the performance of the loop in Figure 1 compare with that of 
the conventional quadriphase loop in Figure 2? 

In the next section, we shall derive the stochastic differential 
equation of operation of the loop in Figure 1. Following that, we 
shall make an analogous derivation for the conventional quadriphase 
Costas loop. (The result of this derivation is available in [2] if 
the signal distortion effects of the arm filters are ignored.) Com- 
paring the two results, we shall then reach the conclusion that the two 
loops are stochastically equivalent , i.e., they have identical stochastic 
differential equations of operation. Thus, an alternate conclusion is 
that the conventional quadriphase Costas loop and its previously shown 
equivalent [2], the fourth power loop, are low signal -to-noise ratio 
practical realizations of the MAP estimate loop for QPSK . 
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2.0 DERIVATION OF THE STOCHASTIC DIFFERENTIAL EQUATION OF 

OPERATION FOR THE CARRIER RECONSTRUCTION LOOP (Figure 1) 

As mentioned in the introduction, the error signal for the quadri- 
phase loop of Figure 1 is formed from the product of the error signal 
and lock detector output signal of a biphase Costas loop, namely, z^(t) 
and z^{t) , respectively. Since the performance of a biphase Costas loop 
whose input is an unbalanced quadriphase modulation plus additive Gaus- 
sian noise has been previously treated [3], we shall be brief in our 
presentation here. Assuming then that the input to the loop of Figure 1 
is the balanced QPSK modulation 

x(t) = v^m^(t) cos' $(t) + v^m^Ct) sin $(t) + n^(t) , (1) 


then from [3], letting ^1 = ^2 

Zjlt) = G(p)E5(t) = 

z^(t) = G(p)e^(t) = 


P, we have for the arm filter outputs 
/P ni 2 (t) - Ng(t) cos (p{t) 

sin (pit) 


I m 
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1 m 


/P m^(t) + N^(t) 


I m 


/P ni2(t) - N^(t) 


+ K, K 
1 m 


/P m^ (t) + N^(t) 


sin (pit) 
cos (pit) 



where the "hats" denote filtering of the corresponding signals, e.g., 
iTii(t) = G(p)m^(t). Thus, the output of the third multiplier is the 
dynamic error’ signal 

z^(t) = z^(t)z^(t) 

sin 2(pit) + 2PitIi (t) iTi 2 {t) cos 2(pit) 


1 m 


m2^(t) - m^^(t) 



N 3 ^(t) -N^^(t) - 2/P m 2 (t) N^(t) -2/P m^(t) N^(t) 
2/P m2(t)N^(t) - 2/P ifi^(t) N3(t) - 2N^(t) N^(t) 


sin 2 (pit) 


cos 2(pit) 


Similarly, taking the difference of the squared filter outputs gives 


( 3 ) 


5 


Z2(t) = - z^{t) 


- - P 


m 2 ^{t) -nii^(t) cos 2(p{t) + 2 P m-| (t) ni 2 (t) sin 2(p{t) 


= 2z^(t) 


Ng^(t) - N^^(t) - 2/P tri2(t) N^(t) - 2y¥ (t) N^(t) 

2»^ m 2 (t)N^(t) -2/P rn^(t)N^(t)-2N^(t)N^(t) 


2(p(t) =i> 2(pit) - 


cos 2(p{t) 
sin 2</?(t)| 
(4) 


where the notation "=^" means "replaced by." As before (see (9) of [3]), 
letting 


Si Jt,2(p{t)'\ = 


N 3 ^(t) - N^^(t) - 2/P m 2 (t) N^(t) - 2/P m^(t) N^{t) 


sin 2(p{t) 


2/P ifi 2 (t) N^(t) - 2/P m^(t) N^(t) - 2N^(t) N^(t) cos 2cp{t) . 


(5) 


then the product of z^(t) and ^2^^^ gives the resulting quadriphase error 
signal 


Zglt) = z^(t)z2(t) 




6 m-|^(t) m 2 ^(t) - m^'^(t) - m^(t) 
+ 4 P^ m-| (t) m2(t) m^^(t) - iri2^(t) 
+ v^[t,4(^(t)] j , 


sin 

cos ^(p{t) 


(6) 


where 


and 


v^[t,4<;i)(t)] = 2 U2[t,2(^(t)] V2[t,2(/?(t) - |] 

+ 2 u^lt,2(p{t) ~ |] V2[t,2(;t>(t)] 
+ 2 V2[t,2(^(t)] V2[t,2(^(t) - I-] 
u^[t,2(p{t)'] = P 

+ 2 P nil (t) m 2 (t) cos 2(p{X) 


(7) 


m 2 ^(t) - m-j^(t) sin 2(p{t) 


( 8 ) 


6 


Substituting (5) and (8) into (7) and simplifying yields 
V4[t,4(^(t)] = Ag(t) sin 4</)(t) + A^{t) cos 4c^(t) , 

where 


(9) 


A (t) = 
s 


N^^(t) - N^^(t) + 4/P m2(t) N^^{t) - tii^(t)N^(t) 


+ 6N^^(t)N^^(t)-12/P 
+ 6?^(t) - N^(t) 


m2(t) N5{t)N^^(t)-m^('t)N^(t) N^^(t) 


- 4P’'/^ni2(t)fi5(t) 
,3/2.c; 




ifi2^(t) - m^^(t) 


3 mi^(t) - m^ (t) 
-2. 


- 24Pm, {t)m.(s'-,) N (t) f^^(t) 


A(t) = 4N (t) fi;"{t) ^4N (t)N^^(t) + 4/P 


C' ' s ' s 

+ 12 P (t) ni2(t) 


m 2 (t) N(?(t) +m^{t) N^^(t) 


Nc^(t) - N3^(t) 


+ 12 PN^(t) N^(t) 


m2^(t) - !T)^^(t) 


- 12,^j^i2(t) N^(t)N^(t) + m^(t)N^(t) N^(t) 


^2. 


+ 4P^/^m2(t) N^(t) 3m^^(t) - ni2‘'(t) 

+ 4 P^^^ (t) N^(t) 3m2^(t) - iri^^(t) 

The instantaneous frequency of the VCO output is related to ZQ(t) by 


d $(t) 
dt 


= Ky [F(p) Zp(t)] + u>Q 


( 10 ) 


( 11 ) 


and hence the stochastic equation of loop operation becomes 


Ad<p{t) 

dt 


4fiQ-KF(p) P^ 


6 m-|^(t) ni2^(t) - ifii^(t) - ni2^(t) 


+ 4 P^ m-j (t) m2(t) 

+ V4[t,4(/?(t)] j , 


m^^(t) -ifi2^(t) 


sin 4(;c>(t) 
cos 4(/?(t) 


A 4 4 

where K = S ’ 


( 12 ) 




a 


m 


m 


§€ 




u 


3,0 DERIVATION OF THE STOCHASTIC DIFFERENTIAL EQUATION OF 
OPERATION FOR THE CONVENTIONAL QUADRIPHASE COSTAS LOOP 

For the quadriphase Costas loop of Figure 2, the input x(t) is 
again given by (1) and the demodulation reference signals r^(t); i = l, 

2,3,4, are, respectively 

r.(t) = /2 sin [$(t) + (i - 1)|] ; i = l,2,3,4. (13) 

Multiplying x(t) by r,^. (t) and ignoring second harmonic terms gives the 
four phase detector outputs which, after arm filtering by G{s), become 

Zi(t) = [/P m 2 (t) - Ng(t)] cos [<;^?(t) - (i-l)|-] 

- K^[/P m^(t) + N^(t)] sin [<^(t) - ( i - 1 ) |] ; i = 1,2. 3,4 

" ^1 ^m{^c^^^ - (i - a^{t) sin [(p{t) - (i -1) J]} 

(14) 

Multiplying z-|{t) and Z 2 (t) gives 

z,|(t)z 3 (t) = K^^K,^^ |[a^^(t) - a^^(t)] +a^(t)a^(t) cos 2c^(t)| . 


(15) 


Similarly, 




sin [2(<;^?(t) -})] 


= K 


+ a^(t) a^(t) cos [2((/?(t) - 

[9(f(t) - as^(t)3 + a^(t)a^(t) sin 2<p{t)j . 

(16) 


Thus, the quadriphase error signal ZQ(t) is obtained by multiplying (15) 
and (16), namely. 


4 .4 


Zo(t) ^ n, z,.{t) = ^ {-[a^2(t)- a2(t)]2 


^ 4 -\n 2 , , , 2,. a , 2,^vV sin 4(^(t) 


1^4 w 4 

+ 1 -- a (t) a (t) [a„^(t) - a.^(t)] cos 

2 C S L b 


(17a) 




3 


8 


or 


Z()(t) 




8 


sin 4(^(t) 


5.=0,2,4 


K. ^ 


•*• - I (-1)^^"^ a^^(t) COS 4(^{t) , (17b) 


8 fi.= l,3 


where is the combinatorial coefficient defined by 


n! 


n"^k " k!(n-k)l 


(18) 


From the defining expressions for a (t) and a (t) in (14), we get that 

0 s 


a/{t) = i(-l)*-\cJ^fi 2 (t))"{N 3 (t)) 


Ji-m 


(19) 


Substituting (19) in (17) gives the final results for Zg(t), namely, 

i 0,-1 


Zo(t) = 


■ '^1 ^ 


8 


I I (-1)*''%C I f {-l)*'"‘tC . C 

(a= 0,2,4 ^ ^ m=0 k=0 m 4 £. k 


X m 2 ^(t) m^'^(t) N^^“'"(t) N^^"^''^(t) j sin 4<;<?(t) 


1 m 


S U=l,3 


l (-1)'“'-^%C I V(-l)'‘"'iC 4 
=1 Q ^ m=0 k=0 K. m 4 a K 


X m^{t) N^^"^"'^(t) cos ^ip{t) . 


(20) 


The terms which are independent of the noise components N (t) and 

^ ... . c 

N^(t) are obtained by letting a = m and k = 4-£= 4-m. Thus, 


4.4 


ZgCt) 


- ') V )„2 


£=in 

k=4-m 


mi'^(t) +ni2^(t) - 6nii^(t) iTi2^(t) 


- 4 m^ (t) ^^(t) .iTii^{ t) - m2^(t) 


sin 4i;£>(t) 

cos 4^(t) j . (21 ) 


Comparing (21 ) with the first two terms of (6), we see that, except for 
a factor of two in gain, the two are identical. Furthermore, evaluating 




the remaining terms in the summations in (20), we get the identical 
signal X noise and noise x noise terms as in (9) combined with (10), 
except again for the same factor of two in gain. (Carrying out the 
algebra to prove this identity is left to the reader,) Thus, 




Kv kJ 

8 


6 m-|^(t) m^ (t) - m-|^(t) - m2^(t) 


sin 4(/?(t) 


+ 4 m-j (t) m 2 (t) mi^(t) - m 2 ^(t) cos 4</?(t) + v^[t,49o(t)] 
. 4 4 

and letting K now equal K-j we get the identical stochastic 

equation of loop operation as in (12). 


4.0 CONCLUSIONS 

We conclude by pointing out that, while we have indeed shown that 
the low signal-to-noise ratio realization of the MAP estimation loop for 
QPSK is stochastically equivalent to the conventional quadriphase Costas 
loop, we have not attempted to give the tracking performance of these 
loops nor compare it to that of a biphase Costas loop. Determining the 
tracking performance of the quadriphase loop of Figure 1 (or Figure 2), 
taking into account the bandlimiting effects of the arm filters, is con- 
siderably more complicated than the equivalent analysis for the biphase 
loop. The principle reason for this stems from the fact that the 
signal x signal, signal x noise, and noise x noise terms now require eval- 
uation of fourth order moments of the filtered signal and noise com- 
ponents. While indeed such evaluation is possible, the resulting expres- 
sions and necessary algebra to arrive at them are quite complicated, 
even for the simplest case of a single-pole (RC) arm filter. Neverthe- 
less, these results, upon completion by the author, will be presented 
in a future report. 


J 
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APPENDIX L 


THE EFFECTS OF RESIDUAL CARRIER ON COSTAS LOOP PERFORMANCE 
AS APPLIED TO THE SHUTTLE S-BAND UPLINK 

by 

Marvin K. Simon 


1.0 INTRODUCTION 

Traditionally, a Costas loop is intended for use in receivers 
which must reconstruct a carrier reference .from an input signal whose 
carrier component is totally suppressed, e.g., a biphase modulated 
carrier. In certain applications, however, the Costas loop is called 
upon to accurately track a signal whose carrier component is not com- 
pletely suppressed. One such application occurs on the Shuttle S-band 
uplink during the time when both data and a ranging subcarrier are 
1 inearly modulated on the same carrier and the data modulation index 
is not tt/2. a simple block diagram which, for the sake of analysis, 
characterizes this situation is illustrated in Figure 1. Included in 
this illustration are the means by which the data and the rang.ing sub- 
carrier are extracted using the in-phase demodulation reference generated 
by the loop. 

Several key questions arise relative to the performance of the 
Costas loop under these unorthodox conditions; 

1. Is the loop capable of successfully tracking the input 
independent of the value of the data modulation index? 

2. Is it possible to extract the ranging subcarrier as shown 
if the data modulation is removed? 

3. What is the additional threshold power-to-noise ratio 
required to operate the loop in the PM (ranging and data) mode as 
compared to the PSK (data only) mode? 

4. What tradeoffs exist between power in the ranging channel 
and loop threshold performance as a function of the data and ranging 
modulation indices? 

This appendix attempts to answer these and other questions related 
to the performance of Costas loops in the presence of residual carrier 
.-by establishing a theory for such performance as a function of the key 
system parameters. For example, we shall first demonstrate that there 
exists a critical data modulation index below which the loop will not 
operate, regardless of the value of signal-to-noise ratio. 
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2.0 SYSTEM MODEL 

Consider the Costas demodulator illustrated in Figure 1 whose 
input signal is of the form 

s{t,e) = sin [oi^t + B^d(t) + 3^ sin (w^^t + Q^(t)) + o] , 


sc r' 


where P is the total received power, Wq is the carrier radian frequency 
and 0 the corresponding input phase to be estimated, 3^ is the data 
modulation index with d(t) the data waveform, and 3^ is 'the ranging 
subcarrier modulation index with the radian subcarrier frequency 
and 0^(t) the tone ranging modulation.’ Using simple trigonometry, 
s(t,0) of (1) may be decomposed into its carrier, ranging, and data 
components, namely, 

s(t,0) = sin [(D^t + 3|jd(t) + 6] cos [3^ sin (to^^t + e^(t))] 


sc r 


+ cos [w^t + 3 |jd(t) + e] sin [3^ sin (uj^^t + e^(t))] 


or, making use of the relations. 


cos [3 sin x] = 5! ^ ’ ^n "" 


1 1 ; n = 0 

W \ n > 0 


sin [3 sin x] = 2 ^ (2n+l)x 


where ’is the n_th order Bessel function of the first kind, then 


s(t,o) = /2P cos 3jj sin (wQt + o) ^n‘^2n (3g) cos [2n (wj^t+ a^(t))][ 

I n=n ’ 


+ sin 3nd(t) cos (iD^t + e)] \ cos [2n (a>^^t + 9^{t)] ] 


0 ^ 


+ cos 3j cos (oogt + 0) 2 I J2n+l^^s^ [(2n+l ) (w^^t + e^(t)) ] 

( n=0 


- sin 3jd(t) sin (wgt + 0 )j 2 I d2n+l^^s^ [(2n+l ) (o)^^t + o^(t)) 

( n=0 


( 4 ) 


The additive channel noise n^it) can be expressed in the form of a 
narrowband process about the frequency of the observed data, namely, 

n^.(t) = v^{N^{t) cos (w^t+e) - N^(t) sin (oiQt+e)} , (5) 

where N (t) and N (t) are approximately statistically independent, sta- 
c s * 

tionary, white Gaussian noise processes withsingle-sided noise spectral 
density Nq w/Hz and single-sided bandwidth < Wq/2tt. 

The input signal plus' noise is bandpass filtered where the filter 
bandwidth is narrow enough to exclude the ranging subcarrier and its 
harmonics. Thus, demodulating the bandpass filtered signal plus noise 
by the quadrature reference signals 

r^(t) = K-| cos (wQt+e) 

r^(t) = K.| sin (w^t+e) (6) 

gives the corresponding phase detector outputs (ignoring second har- 
monic terms) 

" *^1 cos 3^ I J 2 n^^s^ 

+ Ki >^ sin 3dT^(t) | | ^ 

+ Kj^[N^{t) cos q) - Ng(t) sin (pi 

S CO ^ 

J2n^*^s^ [2n(a)gc^ ‘^ ^ 

+ /P s1n 3^d(t) j L ‘^2n^^s^ ^ 

+ [N^(t) sin + Ng(t) cos g?] , (7) 

where is the phase detector gain and cp = 0 - 0 is the loop phase error. 
Assuming that '■o-_/2Tr >> B./2, where B. is the two-sided arm filter 
bandwidth, i.e. , 

B. M |6(d 27rf)|^df , (8) 
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then only the n = 0 term in the summations of ( 7 ) pass through these 
filters. Thus, the arm filter outputs are simply 


z (t) 






/F cos 3 ^JQ( 3 g) - Ng(t) 


sin (p 


+ [/T sin 3 jd(t) Jq(B^) + N^(t) 


cos (p 




^ cos 00(^3) - N3(t) 


^ sin B^d(t) Jq( 6^) + N^(t) 


cos (p 

sin <p 


( 9 ) 


where the "hats" denote filtering of the corresponding signals by the 
arm filters, e.g., d(t) = G(p)d(t). Thus, the output of the third 
multiplier is the dynamic error signal 


ZqW = Z^(t)Z3(t) 


Z .,2 


where 


— {jQ^(Bg) PCd^(t) sin^B^j - cos^ B^] sin 2 <p 

- dQ^(Sg) P sin 2 B^d(t) cos 2 cp + V2[t,2</?]|, 


N^^(t) - N5^(t) + 2/r cos Pd'^O^^s^ ^s^^^ 


+ 2/r sin Bjd{t) ^^(Bg) fi^.(t) 


sin 2 (p 


(10) 


2v^ cos Bj 00(85) N(.(0 - 2/r sin 8^ 3 ( 1 ) J(,(d5) lijU) 
COS 2 (p . (11) 


2 N^(t) NgCt) 


The instantaneous frequency (relative to Wq) of the VCO output is 


related to ZQ(t) by 


^ = K^,[F(p)Z|3(t)] . 


( 12 ) 


Thus, the stochastic equation of operation of the loop becomes 


6 




O 


= -KF(p) 00^(33) P 


~ . 2 2 ' 
sin - cos 3 . 
m a a 


sin Z<p 


+ OffCSg) P sin^ 3^[d^(t) - sin Zcp 

V ; — . 

self-noise 
of data 

- P sin 2Bjd(t) cos 2<p + V2(t,2v?)| » (13) 


# 


where 


= <d^(t)> = I S^(f) |G(j 2.f)l^df 


(14) 


« 


represents the power in the data modulation at the arm filter output 
with S_(f) denoting the power spectral density of the unfiltered data 

m A 9 p 

d(t), and K= Ky. 

The first term in the braces of (13) Is the loop S-curve, while 

the remaining terms are all zero m.ean and thus contribute to the total 

2 2 

noise perturbing the loop. Note that, when cos 3^j = D^sin 3 ^. or, 
equivalently. 


cot^ 3, 


= D, 


m 


the loop S-curve vanishes and thus the loop will not lock at any loop 
signal-to-noise ratio. For values of 3^j less than the critical value 
satisfying (15), the loop locks at = ± (2n+lj tt/ 2, n= 0,1,2,.... For 
values of 3 ^ greater than this critical value, the loop locks at 
= +mr, n = 0,1 ,2 

As a specific illustration of (15), consider the case of Man- 
chester coded data (all of the S-band uplink low and high data rate 
tracking modes) and single-pole Butterworth arm filters (typical of 
Costas loop in network transponder). Then the mean- squared filtered 
data power, 


m' 


is given by [1] 

3 - 4 ex p ( -B I /R^ ) + ex p ( - 2 B . /R^ ) 


m 




(15) 


( 16 ) 


where .ff^ = l/T is the data (coded or uncoded) symbol rate. For a 
single-pole Butterworth filter with transfer function 

^ (17 
and 3 dB cutoff frequency f . the noise bandwidth B. is given by 

C 1 

= 1-f^. (18 

Using (15), (16) and (l8), Table 1 tabulates the critical values of 
data modulation index for the various data symbol rates of interest 
on the S-band upl ink. 


Table 1 





D 

m 

^Pd^crit 

32 kbps (Low 

Data Rate Uncoded) 

134 kHz 

13.155 

0.886 

0.8157 rad 

96 ksps (Low 

Data Rate Coded) 

134 kHz 

4.385 

0.6636 

0.8872 rad 

72 kbps (High 

Data Rate Uncoded) 

308 kHz 

13.44 

0.8884 

0.815 rad 

216 ksps (High 
Data Rate Coded) 

308 kHz 

4.48 

0,670 

0.8848 rad 


3.0 STATISTICAL CHARACTERIZATION OF THE EQUIVALENT ADDITIVE NOISE 

As mentioned above, the equivalent additive noise terms in (13) 
all have zero mean. Furthermore, each has a continuous power spectral 
density component. Since the bandwidth of these processes is very wide 
with respect to the loop bandwidth, it is sufficient to find for each 
one only the power spectral density at the origin which, when multiplied 
by the loop bandwidth, gives the contribution to the total noise power 
of that component. 


8 


9 


In previous analyses of this type [1,2], it was shown that, for 
cases of practical interest, the effect of the self-noise of the data 
modulation on loop tracking performance was negligible. Making this 
same assumption here, and further assuming that cp = 0 (high loop 
signal-to-noise ratio), the equivalent noise of (13) reduces to 


Hg(t) = - jQ^^(e^) P sin 2 6^ d(t) + V2(t,0)- 

= -Jq^( 3^)P sin 2 e^d(t) - 2 cos 6^ ^s^ 


+ 2^ sin + 2N^(t)N,(t) 


09) 


The autocorrelation function R„(t) of n.(t) is easily shown to be 

6 0 






= Jg^(3g) P^ sin^ 2 R^(t) + 4P cos^ 8^ Jq^(b^) R^(t) 

+ 4P sin^ BjjQ^(B3)Rj(T) R^(t) + 4R^^(t) 


( 20 ) 


Q 


where 


Rj{t) = <d(t)d(t + T)) = [ S^(f) |G(j df 


Rd^) = l(t)N (t + T) = N (t)N (t + T) 


O 


^ p 

2 ■ .00 


|G(J 2n f)|^ df , 


( 21 ) 


The equivalent noise power spectral density at the origin is then 


O 


N = 2 

e 


Rg(x) dx 


O 


|0 


= sin^ 2e^|Sj(0) + 4 Nj,P 005^6^,0^(85) 


+ 4N„P sin^ 6 dJo^( 65 : 




S^(f) |G(j2TTf)r df 


.00 


+ 2N, 


|G(j 2TTf)|^ df 


•'—CO 
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( 22 ) 


9 




Letting 

2»f)|'*clf 

|G(j?.f)|'’ df 
IlSj(f) |G(j2»f)|2 df 

and making use of (8) and (14), we can rewrite (22) as 


( 23 ) 


m 


N. = sin2 2s/ “ 


S.(0)' 


T 


w 


2 2 
COS 6 . + Kp, sin B . 
d m D d 


where, in addition, we have defined 

2P 


(24) 


^ = NgB, 

as the signal-to-noise ratio in the arm filter bandwidth. 

4.0 EFFECTIVE LOOP SIGNAL-TO-NOISE RATIO AND TRACKING PHASE JITTER 

From the loop equation of operation in (13), we can see that the 
effective loop signal-to-noise ratio is given by 

(“ni 

- N^Bl 

where B|^ is the single- sided loop tracking bandwidth (typically 200 Hz 
for the S-band network transponder) . Substituting (24) into (26) and 
simplifying results in 

_ P 


(25) 


(26) 


P 9 ' 

4 ^L ’ 




(27) 


« 


where is the loop "squaring loss" defined by 


# 
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. 2 


K 


4 ,. ^ ^^d /^d(°^' 


Jo^(6^)[cos^S^+ D^Kp sin B^l + Rd^o'^s^ 


(28) 


m 


0 


c* 




and R(j^R^/^0 data, signal-to-noise ratio. The comparable result 

for operation in the PSK (data only) mode can be obtained from (28) 
by letting 8^ = 0 and B^ = Tt/2, namely [1], 

2 


D_ 


m 


\ ■ 

D Kn + — 
m D p. 


Againj for Manchester coded data and single-pole Butterworth 
arm filters, we have that 5^(0) = 0 and Kj^=l/2. Thus, (28) simplifies 
to 


Jo^(65)[cos2sj + Dn,KD Sin^6d]^- ^ 


(29) 


(30) 




What remains is to characterize the tracking phase jitter per- 
formance. In the linear region (high loop signal-to-noise ratio), 

2 

the mean-squared phase jitter of the loop phase error 2cp is 
given by 


1 


2<p 


p *5i 


(31) 


Since the demodulation reference for the data and ranging signals is 
at u)q rather than 2 u)q, then the jitter which degrades the power in 
these components has variance 




± 2 
4 '’Zcp 


1 


pS, 


( 32 ) 


m 


5.0 RANGING CHANNEL PERFORMANCE CHARACTERISTICS 

5.T Data Modulation On 

The ranging subcarrier with its tone modulation is extracted 
from the input signal plus noise by demodulating this input with the 


n 


# 


# 




t 


# 




! m 


m 






in-phase reference (r (t)) generated by the loop in a wideband phase 
detector and then bandpass filtering the output of this detector (see 
Figure 1). The signal and noise components at the filter output are 
obtained from (4), (5) and (6), and are given by 

y(t) = cos (o)Qt+ e)[2 sin (w^^t + 9^(t))] 

+ *^[N^(t) cos (copt + e) - Ng(t) sin (ujQt + e)]| ^ K-j cos (wQt + e ) 
COS Bfj J-|(Bg) cos (p sin (u^^t + o^(t)) 

+ N^,(t) cos cp - Ng(t) sin <p^ , (33) 

where the "hat" now denotes filtering by the ranging channel bandpass 
filter (bandwidth equal to B^). Thus, the, signal power into the 
ranging channel is 

= 2P cos^ J^^(3,) cos^<^ . (34) 
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Note from (34) that, as the value of 3^^ is increased beyond its 
critical value, the carrier becomes more suppressed and thus the Costas 
loop tends to track better. However, as 3^ is increased, cos 3^ 
decreases and thus the power in the ranging channel tends to 
decrease because of this. Since the power in the ranging channel is 
proportional to both cos 3^ and the accuracy of Costas loop tracking 
(through the factor cos^ cp ) , a tradeoff exists with regard to the 
selection of the data modulation index. 

5.2 Data Modulation Off 

When the data modulation is removed, i.e., = 0, then as pre- 

viously mentioned, the loop will now lock up around (/? = tt/2 (as opposed 
to (p = 0) . This is equivalent to saying that the in-phase and quadrature 
demodulation reference signals switch roles. Thus, since the ranging 
subcarrier component in the input signal s(t,e) is now demodulated by 
the quadrature carrier reference signal, then from (34), we have that 
(for large loop signal-to-noise ratio) P^=0. The conclusion to be 
reached then is that, in order to extract range information as in 
Figure 1, the data modulation cannot be turned off at the transmitter. 

6.0 PM VERSUS PSK THRESHOLD TRACKING PERFORMANCE 

Tracking threshold is typically specified as the minimum input 
signal level for which the mean time to loss of lock is greater than 
or equal to 10 seconds. This condition is uniquely specified by the 
loop signal-to-noise ratio or, equivalently, Thus, two modes 

of operation, e.g., PM and PSK, will have identical tracking thresholds 
if the total power- to-noise ratio P/Ng is adjusted in one relative. to 
the other to produce the same in both. 

Using (29), (30) and (31), one can determine In accordance with 
the above equivalence the additional P/Ng required in the PM mode to 
produce the sarne tracking threshold as for PSK operation. Letting 
P'/Hg (> P/Ng) denote the power-to-noise ratio required in the PM 
mode, then, 


0 , , Vl 

“m 4 P 


4NqBl 


^o(6s) 

N B . 

Oo(Sj) [cos^ Sj + D^KpSin^Bj]+ 


o 
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Solving for P'/Ng in terms of P/Ng gives the quadratic equation 


where 


A(P7Nq)^ + B(PVNg) + C = 0, 




0 V 7 ?~1/P\ ^m 

6 = -00(63) ^os 6 j + D„,KoSin 6 jJ ^ ^ 


”0 / P 


\\) „ „ /o“i • 

“m '^D TT' 


The solution to (41) is well known to be 


'_B_\2 £ 

.2A/ " A 


Note that, if arm filter degradations are ignored, i.e., = Kg = 1 

(corresponding to ->■“), then 


/ \ r 

p \ I 2 /’ '' 

0203 ^( 63 ) cos^ 26 , 


m 


which is bounded below by 


^0 Jg^(3g) [cos 26^ 


Npr 1 


and bounded above by 


0 ^^d 


for -?■ P— >> 1 

NqB. ' 
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7.0 APPLICATION TO THE S-BAND NETWORK TRANSPONDER DESIGN 

The baseline design of the S-band network transponder under 
development by TRW calls for S5=1-0 rad and B^ = l.l rad. Performance 
tests on the engineering model, as reported in the network transponder 
CDR package [4], reveal the following tracking threshold signal levels 
in the PSK mode (through TDRS): 


Table 2 


Data Rate 

Temperature 

Signal Level 

Low Data Rate 

70°F 

-104.2 dBm 

(32 kbps uncoded or 

120°F 

-102.4 dBm 

96 ksps coded) 

-20'^F 

-104.2 dBm 

High Data Rate 

70°F 

-103.7 dBm 

(72 kbps uncoded or 

120°F 

-101 .4 dBm 

216 ksps coded) 

-20°F 

(SSP will not hold 
lock -100. dBm) 


The input noise power spectral density is -151 dBm/Hz. Thus, 
the maximum measured P/Nq is 48.6 for the low rate and 49.6 for the 
high rate modes. Assuming these threshold values of P/Nq for PSK 
operation through the TDRS, Table 3 tabulates the values of various 
system parameters leading up to the calculation of P'/Ng from (43). 

Also indicated in the table is the dB increase in power- to-noi se ratio 
required to go from the PSK to the PM mode through the TDRS. 

Table 4 illustrates the power-to-noise ratio P^/Ng in the ranging 
channel [as computed from (34)] for the same parameters as in Table 3 
and operation in the PM mode. 

8.0 PERFORMANCE AS A FUNCTION OF DATA AND RANGING MODULATION INDICES 

The tracking threshold of the Costas loop when operating in the 
PM mode can be reduced (relative to that in the PSK mode) by reducing 
the ranging modulation index 3^ and/or increasing the data modulation 
index (i.e., suppressing the carrier more) . Unfortunately, however, 
both of these changes also reduce the power-to-noise ratio in the rang- 
ing channel. To see these effects quantitatively. Figure 2 illustrates 


o 
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Table 3 



P/NQ=48.6 d3, Bg = 1.0rad, 3^1= 1.1 rad 



Data Rate Pf/N^ (dB) p. (dB) P'/Nq (dB) AP/Nq (dB) 


Data Rate Pf/N^ (dB) p. (dB) P'/Nq (dB) AP/Nq (dB) 

32 kbps 0.829 3.549 -4.642 54.23 5.63 
96 ksps 0.5078 -1.223 -4.642 55.01 6.41 
72 kbps 0.8326 0.0267 -8.257 53.83 5.23 
216 ksps 0.5166 -4.745 -8.257 54.50 5.90 


Table 4 

P/NQ=48.6dB, 3g=1.0rad, rad , B^ = 200 Hz 


Data Rate 

P^/No (dB) 

32 kbps 

43.21 

96 ksps 

43.97 

72 kbps 

42.78 

216 ksps 

43.41 



the tracking threshold power-to-noise ratio for PM operation as a func- 
tion of ranging modulation index for fixed data modulation index, and 
Figure 3 illustrates the ranging channel power-to-noise ratio as a 
function of these same parameters. The value of P/Ng for PSK opera- 
tion is again chosen equal to 48.6 dB and the results are given for 
both uncoded and coded low data rates. We observe from Figure 2 that, 
as approaches tt/2 (fully suppressed carrier), the tracking threshold 
becomes virtually insensitive to data rate. A similar phenomenon is 
observed in Figure 3 for the ranging channel power-to-noise ratio. 

Using the information in these curves, along with the ranging channel 
bandwidth and the required signal-to-noise ratio in this bandwidth 
turned around to the ground, one can select values of 3^ and 3^ to 
meet the additional requirements on increase in tracking threshold for 
PM versus PSK operation. 

9.0 CONCLUSIONS 

The operation of a Costas loop in a PM mode produces performance 
degradation (relative to PSK operation) because of the presence of 
residual carrier. This appendix has numerically evaluated such degra- 
dation for parameters of interest on the Shuttle S-band uplink. The 
corresponding performance of the ranging channel is also given as a 
function of these same system parameters, thus allowing the necessary 
tradeoffs (choice of modulation indices) to be made between such per- 
formance and that of the Costas loop when tracking under threshold 
conditions. 
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